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Many trial wavefunctions for fractional quantum Hall states in a single Landau level are given 
by functions called conformal blocks, taken from some conformal field theory. Also, wavefunctions 
for certain paired states of fermions in two dimensions, such as p^, + ipy states, reduce to such 
a form at long distances. Here we investigate the adiabatic transport of such many-particle trial 
wavefunctions using methods from two-dimensional field theory. One context for this is to calculate 
the statistics of widely-separated quasiholes, which has been predicted to be non- Abelian in a variety 
of cases. The Berry phase or matrix (holonomy) resulting from adiabatic transport around a closed 
loop in parameter space is the same as the effect of analytic continuation around the same loop 
with the particle coordinates held fixed (monodromy) , provided the trial functions are orthonormal 
and holomorphic in the parameters so that the Berry vector potential (or connection) vanishes. We 
show that this is the case (up to a simple area term) for paired states (including the Moore-Read 
quantum Hall state) , and present general conditions for it to hold for other trial states (such as the 
Read-Rezayi series). We argue that trial states based on a non- unitary conformal field theory do 
not describe a gapped topological phase, at least in many cases. By considering adiabatic variation 
of the aspect ratio of the torus, we calculate the Hall viscosity, a non-dissipative viscosity coefficient 
analogous to Hall conductivity, for paired states, Laughlin states, and more general quantum Hall 
states. Hall viscosity is an invariant within a topological phase, and is generally proportional to the 
"conformal spin density" in the ground state. 



I. INTRODUCTION 

A. Background and motivation 

There has been renewed interest in the past few years 
in non-Abelian quantum states of matter, both in the 
original setting of quantum Hall states [l|, and also in 
other systems including ones in which the symmetries un- 
der time reversal and parity are unbroken (see e.g. Ref. 
0). Briefly, a non-Abelian phase of matter is a gapped 
(topological) phase in which there are quasiparticle exci- 
tations over the ground state, the adiabatic exchange of 
which produces a matrix effect on the state of the system, 
with the matrices corresponding to distinct exchanges 
not all commuting (thus the term "non-Abelian"). This 
requires that there be a degenerate space of states when 
there are quasiparticles at well-separated positions in the 
system. Given a reasonably local Hamiltonian, such be- 
havior can only occur when the dimensionality of space is 
two, for topological reasons. The basic example proposed 
by Moore and Read (MR) [l| is a paired state, which can 
be viewed as a p-wave, or more precisely a Px + ipy- 
wave, Bardeen-Cooper-Schrieffer (BCS) [3] paired state 
of spinless or spin-polarized composite fermions in zero 
or almost zero net magnetic field. It turns out that much 
of the physics of this state is also found in such paired 
states of ordinary fermions Q ■ In this paper we address 
both of these situations together, not only quantum Hall 
systems. Much of the current interest in these systems 
is driven by their potential use for topological quantum 
computing 

The major issue that we address in this paper is the 
derivation of the non-Abelian statistics of the quasi- 



holes (or vortices) when they are exchanged adiabatically, 
where adiabatic transport is calculated using trial forms 
of the wavefunctions. The trial wavefunctions are taken 
to be "conformal blocks" obtained from some conformal 
field theory (CFT), as in MR. 

A central idea of MR is that the adiabatic effect of an 
exchange of quasiholes in the trial wavefunctions given by 
conformal blocks is the same as the effect inferred from 
simple analytic continuation of the wavefunctions viewed 
as functions of the quasihole coordinates, with the par- 
ticle coordinates fixed. In contrast to analytic contin- 
uation, adiabatic transport involves integration over the 
particle coordinates for each infinitesimal time step, since 
this defines the inner product in the Hilbert space. In 
several examples with Abelian statistics, this was shown 
in MR to give the same result, but the result was not 
demonstrated for the trial states with non-Abelian statis- 
tics. Here we provide a detailed explicit calculation for 
the cases of two or four quasiholes in the MR paired state 
(with the charge part of the wavefunctions removed) ; the 
four quasihole case is the first to exhibit non-Abelian 
statistics. Finally, we also give a general criterion for 
other states, and discuss when this may hold, with ex- 
amples. 

An additional issue that we address is adiabatic varia- 
tion of the geometry in a closed finite system, for example 
a parallelogram with periodic boundary conditions (i.e. 
a torus). There is a discrete group of transformations 
that map the geometry to one equivalent to the original 
(called the modular group). It is of interest to perform 
these adiabatically also. It turns out that in these cases 
the relevant Berry connection (vector potential on the 
space of inequivalent geometries) has non-zero curvature 
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(field strength). Because varying the aspect ratio varies 
the strain on the system, this response can be related to 
a non-dissipative viscosity coefficient [ij, which we call 
Hall viscosity (earlier it was termed "odd" , or "antisym- 
metric" viscosity [6j, or "Lorentz shear modulus" [31 )■ 
Hall viscosity is so-named because it is the natural ana- 
log in viscosity of the Hall conductivity. We show that 
this viscosity arises in quantum Hall fluids and in paired 
superfluids, and its magnitude is always proportional to 
a spin density. This result should be of wide interest. For 
trial wavefunctions given by conformal blocks, the spin 
per particle is given by the conformal weight. 

In the remainder of this Introduction, we explain the 
central issues in more depth, introduce some basic ar- 
guments that will be crucial later, and also review some 
earlier work relevant to the problems of adiabatic statis- 
tics and Hall viscosity. 



B. Adiabatic transport and quasiparticle statistics 

Here for the time being we are concerned primar- 
ily with quantum Hall states. The basic examples are 
Laughlin's trial wavefunctions for a ground state with n 
quasiholes (each of the lowest possible non-zero fractional 
charge), which in their original form are Q 

Y[{zk - wi) l[{z, - z,)Qe'i (1.1) 

k,l ij 

Here the exponent Q determines the filling factor of the 
state in the thermodynamic limit as the particle number 
N goes to infinity, Zi {i — 1 , . . . , N) are the complex 
coordinates of particles, and wi are the complex coordi- 
nates of the quasiholes. When Q is even, the particles 
are bosons, and when Q is odd they are fermions, such 
as electrons. In either case, the functions are for particles 
in the lowest Landau level (LLL). 

Now we consider the adiabatic calculation of statistics. 
The adiabatic calculation is meaningful if the excitation 
spectrum of the system above the possibly degenerate 
quasiparticle state subspace is gapped (with the gap go- 
ing to a constant as the system size goes to infinity). It 
can also be valid if any gapless excitations are sufficiently 
weakly coupled to the quasiparticles. For example, in a 
quantum Hall system with an edge, there are gapless edge 
excitations. When quasiparticles are present within the 
region occupied by the particles, we will assume through- 
out this paper that the system (i.e. N) is large enough so 
that during the exchange the distance of the quasiparti- 
cles from the edge is large compared with their separa- 
tion, which is itself large compared with the microscopic 
correlation length scale (which frequently is of order the 
particle spacing). In this case the edge excitations effec- 
tively decouple from the calculation. 

The pioneering calculation of adiabatic statistics in 
a QH system was performed by Arovas, SchriefFer and 
Wilczek [l3|. (An interesting alternative approach 



has been formulated recently [ll[, and applied to non- 
Abelian states [l^ as we were finishing this paper.) 
We follow the framework of their method, but gener- 
alize to allow for non-Abelian statistics. Thus suppose 
that we have a space of degenerate states with basis 
|5'(j(u')), where a runs over some finite set, and we will 
use w to stand for the set of n quasihole coordinates w/, 
I = 1, . . . , n. We assume that the basis is orthonor- 
mal, a{w)\^ b{w)) = Sab for each w. As w varies, these 
vectors sweep out a subspace of a common multiparticle 
Hilbert space in which they all live [H,[l3|- The states for 
any fixed w are supposed to be degenerate in energy, and 
we will compute the adiabatic statistics within this sub- 
space. We neglect the "dynamical phase" which depends 
on the energy and the total time taken for an exchange. 
The desired adiabatic phase or matrix is the holonomy, 
which can be written [l^ 0, [l5| 

V exp i (j) {Aiu ■ dw + Ayj ■ dw) (1.2) 

where C is a closed directed path that begins and ends 
at a base point W(o) i in f configuration space of the 
quasiholes. (This configuration space can be thought of 
as C" where C is the complex plane or Riemann sphere, 
minus the "diagonal" on which Wk = wi for some k, I, 
k ^ I, modulo the action of the permutation group Sn for 
the identical quasiholes; if the quasiparticles are not all 
of the same type there are obvious modifications of the 
permutation group.) V is the path ordering operator, in 
which the matrix for a later point on the path is to the 
left of earlier ones, ■ dw — J2i ^w,idwi, and A^.i is a 
matrix. It is possible to change to a different basis set 
= \'^a{w))Mab{w) as a function of position on 
the path C, where M{w) is a unitary matrix for each 
point w on the path C, in order to preserve orthonor- 
mality. This is referred to as a change of gauge. The 
expression is correct as written when the vector poten- 
tial A refers to a single gauge choice for the whole path 
C [as in the case of the states in eq. (jl.ll) above]. In 
general such a choice may not be convenient, and then 
one must use patches with a gauge choice on each patch. 
The patches overlap, and in the overlap region there is 
a transition function (gauge transformation), which is a 
unitary matrix M. The holonomy is then the matrix 
composite of the path-ordered exponentials for the part 
of the path within each patch, with gauge transforma- 
tions by the transition matrices inserted in between. It 
is frequently the case that a single transition function 
M is sufficient, and the transition can be located at the 
basepoint W(o) • Then the adiabatic transport maps the 
basis states |^'6(w(o))) to J2a\'^a{w(^o)))Bab, where the 
holonomy is given by 

B ^ MVexpi (j) {A^ ■ dw + A-:uj- dw). (1.3) 

Defined this way, the holonomy is gauge invariant in all 
cases, up to conjugation by a unitary matrix that cor- 
responds to a change of orthonormal basis at W(o). We 
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note also that the use of a different base-point on the 
same path C conjugates the holonomy by some unitary 
matrix, and has no effect on the structure of the braid- 
group representation. 

The Berry connection is given in components by 



dwi 



dwi 



(1.4) 
(1.5) 



Then A:^ i{w) = Aw.i{w)\ which ensures that the holon- 
omy is unitary {A^ is not holomorphic in w in general). 
For functions in the LLL, the inner product is the usual 
(*a|*fc) = /n'^^^»*a*b- In the calculation in Ref. ^9,], 
the trial quasihole states used were the wavefunctions 
(jl.ip which are not orthonormal with respect to this inner 
product (for Laughlin quasiholes, there is only a single 
state for n quasiholes, so the label a is dropped). Arovas 
et al. appeared to neglect this point p^], but arrived at 
the correct answer nonetheless. The result is that if two 
quasiholes are exchanged around a counterclockwise path 
that does not enclose any other quasiholes (the others 
stay at fixed positions throughout), the adiabatic phase 
change is e^'"!'^ . 

The central point of the calculation is that the neces- 
sary inner products can be evaluated by using the plasma 
mapping of Laughlin [s^ , plus the fact that screening oc- 
curs in the Coulomb plasma provided Q is not too large. 
We reformulate this as follows (following a line of ar- 
gument begun by Halperin [lB|)- Consider instead the 
following wavefunctions, which represent the same quan- 
tum states as because they differ only by functions 
of the parameters w: 

*(wi,...,-u;„;2;i,...,ZAr) = 



h<X i.k 



Kj 



X e 



(1.6) 



(here and below we take the notational liberty of using 
some indices more than once in distinct factors, which 
have been separated here by the dots •). The modulus- 
squared is the Boltzmann weight for a classical plasma of 
charges Q at all Zi, 1 at all Wk, with a neutralizing back- 
ground of density l/27r, impurities at fixed positions Wk, 
with the two-dimensional Coulomb interaction potential 
between unit charges taken to be minus the natural loga- 
rithm of the distance-squared, and with temperature Q. 
Thus the overlap integral is the corresponding partition 
function: 

Z{wi, . . . ,w„) = |||^'(it;i, ■ • ■ , w„))||^ = 



jY[(fzi exp^ 



^ln|wfc - + Q^ln\zi - 



i<j 



Wk\ 



(1.7) 



Given the partition function, we can define the free en- 
ergy F as = — In Z (or as this times the temperature Q, 
but this makes no difference for us). Then by screening, 
for sufhciently large separation of the quasiholes, this free 
energy goes to a constant at sufficiently large separation 
of the quasiholes. This constant has the form 



= A/o(g) + n/qh(g), 



(1.8) 



k<l 



consisting of an extensive background term, which is the 
area A occupied by the particles plus quasiholes times 
a constant foiQ), plus a "defect free energy" /qh(Q) for 
each quasihole. (In a more general situation with dif- 
ferent types of quasiparticles, these latter terms would 
be different for each type of quasiparticle, as well as de- 
pending on the underlying ground state.) Both constants 
foiQ) and /qh(Q) are well-defined but Q-dependent, 
given the definition above, but have no universal signifi- 
cance. 

Now we can reformulate the statistics calculation. Sup- 
pose that, in addition to being (ortho-)normalized, the 
wavefunctions of the states |^'a(w)) are holomorphic in 
w, as the Halperin functions above are (except for the 
Gaussian factors, which we ignore here, but comment on 
afterwards, and except on the diagonal Wk = wi, which 
is to be avoided). Then it follows that 

A^,i.ab{w) = {^a{w) |*t(w) ) = 0. (1.9) 

OWi 

Then the holonomy is given entirely by the transition ma- 
trix or matrices: B — M. For the wavefunctions ()1.6|1 . 
the transition matrix is just a number M of modulus 1 
that is required to transform back to the original sheet 
(or gauge) after making an exchange, due to the wave- 
functions not being single valued in the WfcS. For the ex- 
change of two quasiholes along a path not enclosing any 
others (which stay at fixed positions throughout), this 
gives the phase M = e™^'^. It is important to empha- 
size that when applying this argument to wavefunctions, 
they should be functions in a common Hilbert space for 
all w, and with the integration measure independent of 
the parameter w being varied. 

The fact that the normalizing Gaussian factors are not 
holomorphic in w means that the calculation also pro- 
duces a phase factor e^^'^'-^'^^'^ , where A{C) is the area 
enclosed by the loop C Q . This can be interpreted as the 
fractionally-charged quasiparticles detecting the back- 
ground magnetic field Q through an Aharonov-Bohm 
phase, though in fact it is Q times the particle density 
that they detect, which happens to give the same re- 
sult because the particle density is uniform, with filling 
factor 1/Q, outside the quasiholes. (If the particle den- 
sity were not uniform, the normalizing Gaussian factors 
in ws would be modified.) This effect is ubiquitous for 
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QH states. We will comment on it further in connec- 
tion with systems that are not QH systems. Other than 
this effect, the holonomy depends only on the homotopy 
class of the path C in the configuration space, not on the 
precise path; that is, it is invariant under small defor- 
mations of the path such that C docs not pass through 
the diagonals wu = wi during the deformation. Then 
the holonomy gives a unitary representation of the braid 
group, acting in the space of degenerate states at W(o), 
times the path-dependent factor e'"^*^'^^/'^ which we will 
frequently just ignore. 

The transition matrices M result from the behavior 
of the wavefunctions under analytic continuation of the 
quasiholes, with the particle coordinates held fixed (or 
even if they are not fixed). In the context of solutions 
of differential equations, this is called monodromy, and 
we follow this terminology here. Halperin [l^ noted that 
the monodromy of his functions (|1.6p suggested fractional 
statistics, but did not perform the adiabatic calculation. 
The use of the present approach in the adiabatic calcula- 
tion gives the statistics without the further calculations 
performed by Arovas et al. Q, in particular, confirm- 
ing the sign of the phase (alternatively, their calculation 
Pi can also be used to produce the normalization of the 
Laughlin quasiholc functions [13] )■ 

The statistics calculation for the Laughlin quasiholes 
generalizes easily to wavefunctions that correspond to 
multicomponent Coulomb plasmas Here we wish 

to generalize it to non-Abelian cases, specifically those 
in which the trial wavefunctions are conformal blocks 
from CFT. These are holomorphic in the quasiholc co- 
ordinates, and have non-trivial, sometimes non-Abelian, 
monodromy; the Halperin form (jl.6p of Abelian quasiholc 
wavefunctions can also be viewed as conformal blocks [l[ . 
The goal will be to see whether the effect noted above for 
the latter states, that holonomy equals monodromy^ holds 
for these trial functions, as conjectured by MR. The re- 
sult will hinge on whether the screening property in the 
plasma mapping for the Laughlin states, which makes the 
Halperin functions (conformal blocks) orthonormal, also 
generalizes when other conformal blocks are used as trial 
wavefunctions. The answer will be yes in some cases. 

In addition, we apply the same formalism to consider 
adiabatic variation of the geometry of the system. We 
will show later that the curvature of the Berry connection 
is non-zero, and this determines the Hall viscosity. We 
also consider the holonomy around a loop corresponding 
to a modular transformation. In some cases we can ap- 
ply a similar argument, that the orthonormalized wave- 
functions are almost holomorphic in the relevant complex 
parameter r, except for a very simple non-holomorphic 
part. 



C. Conformal blocks as trial wavefunctions 

Next we discuss how conformal blocks from CFT can 
be used as trial wavefunctions, following MR. For back- 



ground on CFT, see Ref. [TJ. 

The trial wavefunctions for QH systems that we will 
study in this paper take the form 

^'a(wi, . . . ,w„;zi, . . . ,ZAr) = (1-10) 

^' charge • {i'izi) ■ ■ ■ -iJj{zn)t{wi) ■ ■ ■ T{Wn))a, CFT- 

Here again coordinates Zi are those of particles (either 
bosons or fermions), so '^a is single valued and either 
symmetric or antisymmetric in these variables, while co- 
ordinates wi are those of quasiholes. The label a again 
runs over a basis for a space of functions. On the right 
hand side, ^'charge is independent of a, and is a func- 
tion of the same coordinates of similar form to eq. (|1.6p . 
though the exponents (or charges of the particles and 
impurities in the corresponding plasma) may take other 
rational-fraction values. We note that the exponent in 
the particle-particle factors is the inverse of the filling 
factor, I' = P/Q, and the Gaussian in the particle coor- 
dinates is always as in ()1.6|) . The values of the exponents 
in the charge part are determined by the requirement 
that the whole function '^a be single valued and (anti-) 
symmetric in the particle coordinates; this may always be 
done consistently with a plasma form, thanks to consis- 
tency properties of conformal blocks. The exponents are 
determined only up to addition of integers, and we usu- 
ally consider the smallest possible positive values, which 
gives the highest value of the filling factor, and the lowest 
value of the quasiholc charge; further the field r is usually 
also chosen, given the CFT, to obtain the lowest possible 
quasiholc charge for a given possible filling factor. 

The expectation value J-a{wi, . . zi, . . .) = {■ ■ ■)a, cft 
stands for a conformal block in some CFT, in which tp 
and T stand for fields. The notation is somewhat formal, 
because the function is generally not single- valued, and a 
sheet should be specified. The function is holomorphic in 
the zs and ws off the diagonal on which some zs and/or 
ws coincide. The field ip must have the property that its 
monodromy is Abelian, which means that there is a sin- 
gle primary field (of some chiral algebra) in its operator 
product expansion (ope) with itself (we omit non-zero 
ope coefficients, which play no role here): 

V'(z)V'(0) V*(0) + . . . , (1.11) 

as z — > 0, where ^p* is another field, and the . . . stands 
as usual in CFT for terms smaller by positive integer 
powers of z as z ^ 0, which are descendants of ip* under 
the chiral algebra. Moreover, ijj must generate in this way 
only fields that are also Abelian. Further, there must be 
exactly one term in the ope of V' with r: 

V'(z)r(0) - z-''*-'^.+'»x-^*(o) + . . . ^ (1.12) 

where again t* is another field, is the conformal 
weight of ip (and similarly for other fields), and the same 
is true for r* and so on. These requirements mean that 
■0 is what is called a "simple current": the operation of 
taking the ope of ip with the primary fields in the the- 
ory (including ip itself) just permutes the primary fields. 
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This has the effect of guaranteeing that the full function 
has the stated properties. The ope of r with itself 
may be nontrivial, containing terms that do not differ in 
conformal weight only by integers: 

t(z)t(0) - z-2'»x + ^(i)^(l)_^^ . + Z-2'^-+^(^)T(2)+. .. + ..., 

(1.13) 

where, depending on which CFT is used, any number of 
distinct primaries could appear on the right. Those ap- 
pearing with nonzero coefficients may be summarized by 
the fusion rules (analogous to the Clebsch-Gordan for- 
mulas for SU(2) tensor products), which generally have 
the form 

7 

in which the fields (jja run over the full set of primary 
fields in the CFT used {iIj, ^p* , . . . , as well as r, r*, t'-^', 
. . . , will be among these), and the Af^p = J^Ja ^""^ "^O"^" 
negative integers (which may be larger than 1 in some 
cases). The product here is formal and simply refers to 
terms in an actual ope. We define 4>o = 1, the identity 
operator. For an Abelian field (one obeying Abelian fu- 
sion rules), A/'^„ is equal to 1 for one value of /3 (and 
zero otherwise), and similarly for iterated products; oth- 
erwise the field is non- Abelian. (As explained in Ref. [l| , 
^' charge Can also be viewed as a kind of conformal block 
in the CFT of a single scalar field with the role of the 
fields V', T played by charged fields (exponentials of ip), 
but which also includes a chiral version of the neutral- 
izing background charge density |8j].) A rational CFT is 
one in which there is a chiral algebra^ which is either the 
Virasoro algebra or an algebra extension of it, obtained 
from ope's of a finite set of Abelian conformal fields, and 
a finite set of primary fields, defined generally as confor- 
mal fields that generate irreducible highest-weight rep- 
resentations of the chiral algebra (by operator products) 
[TtI [18!| . Examples of non-trivial chiral algebras include 
the affine Lie algebras. At many places we will need to 
assume that the CFT is rational, as is the combined one 
that includes the charge part. Irrational examples will 
be discussed at the end of the paper. 

The conformal blocks are multisheeted functions when 
the exponents such as 2ft,^ — h^* are not integers, and 
because of non- Abelian fields r. The multisheetedness 
due to fractional exponents from ip is no worse than in 
the case of the Halperin functions (|1.6p . and produces 
only a phase factor when these are exchanged or encircle 
a r, so not a linearly-independent function of the par- 
ticle coordinates (as the latter vary over some open set 
that lies on a single sheet). Finally, a labels a basis for 
the space of linearly-independent conformal blocks asso- 
ciated with a given correlation function; the range of a is 
the dimension of the space of blocks, which can be cal- 
culated for the case of CFT on the sphere by repeated 
use of the fusion rules. By definition, this is the number 
of linearly-independent functions of all the variables z 
and w as they vary over some open set. When viewed as 



many-particle wavefunctions, we would instead count the 
number of linearly-independent functions of the particle 
coordinates z for fixed w. In general, the latter number 
might be less than the number of blocks (for example, 
when iV = 0). However, in examples these numbers do 
coincide when the particle number N is large enough, 
and we will assume this from here on, as large TV is the 
case of interest anyway. More generally, of course, one 
could have more than two types of primary field in the 
correlator, which could represent the particles and more 
than one type of quasihole; such could result from use of 
the ope starting from one type of "basic" quasihole. 

In CFT, correlation functions are constructed from 
combinations of blocks and their conjugates, for exam- 
ple 

{ip{zi,zi) ■ ■ ■i/j{zn,zn)t{wi,wi) ■ ■ •t(wi,u;i))cft 
= Y,\^a{wu...;z,,...)\^, (1.15) 

a 

where the sum is over the basis for the space of blocks. 
Such an expression is supposed to represent a single- 
valued correlation function of local operators in a CFT; 
the local operators tp{z, z) and so on that appear here 
differ from the chiral versions ipiz) that appeared before, 
and are related roughly by 1/1(2,2) = ip{z)'ip{z), but also 
require the sum over the blocks. The "diagonal" form 
given is single-valued if the monodromy of the space is 
given by unitary matrices Bab h^ the basis used. Such a 
diagonal form is always available in a unitary theory, and 
even in non-unitary rational theories, such as the Vira- 
soro minimal models of Belavin, Polyakov, and Zamolod- 
chikov (BPZ) [13, dl. More generally, when the CFT is 
formulated on a surface of higher genus, the number of 
blocks is in most cases larger than one even when no 
fields are inserted in the correlator. For N, n > 0, the 
single label a labels the full space of blocks. A diagonal 
theory, in which all correlators are given by diagonal ex- 
pressions like that above exists in most cases, and is valid 
for surfaces of any genus; it is referred to as the diago- 
nal modular invariant theory, where modular invariance 
refers to the case of the torus (genus 1), with no fields in- 
serted. At this point we should emphasize that what we 
call a correlation function (or correlator) here is really a 
partition function for the CFT on the given surface, with 
fields ip and r inserted at specified points. This differs 
from the more usual use of normalized correlation func- 
tions, in which one divides by the partition function for 
the case of no field insertions. For example, when -0 is a 
basic field of the theory, described by an action S[ip], the 
theory can be described by a functional integral, 

Z = J V[ij]e-^'^^\ (1.16) 

and the unnormalized correlation function is 

{ip{zi,zi) ■ ■ ■^I){zn,zn)t{wi,wi) ■ ■ ■t{wi,wi))c-ft 
= Z{wi,wi, ...■,zi,zi,.. .) 
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= J V[t/j]^{zuzi) ■ ■ ■T{wi,wi) ■ ■ ■ e-'^^'^\ (1.17) 

whereas usually one would define the correlation func- 
tion to be Z{wi, wi, . . . ; zi, zi, . . .)/Z. Even leaving aside 
a factor of the form exponential of minus a free energy 
proportional to the area of the surface, this makes a dif- 
ference for genus > 1 , because the denominator is in most 
CFTs a sum of more than one mod-square conformal 
blocks. The exceptions to the latter are "holomorphic" 
CFTs, in which the only primary field of the chiral alge- 
bra is the identity operator 1; an example is the current 
algebra or Wess-Zumino-Witten theory for the Lie group 
Eg at level 1. 

As an example of the conformal blocks and their use 
as trial wavefunctions, we give the MR state, for the case 
of no quasiholes, in which ip is a Majorana fermion field, 
with conformal weight — 1/2. The ground state trial 
wavefunction for the sphere or infinite plane contains the 
conformal block [l[ 

(V'(zi)V(22)---)lsi„g=Pf^^, (1.18) 

Zj 

where the PfafFian is defined for any even-by-even anti- 
symmetric matrix with matrix elements Mij by 

Pf M,, = ^(Mi2A/34 • • • Mn-im). (1.19) 

where the antisymmetrizer A sums over all permutations 
that produce distinct pairings (i, j), times the sign of the 
permutation. For quasiholes, one uses [H t = cr, the spin 
field of the critical Ising model; 1, ip and a are the only 
primary fields in the Ising (or Majorana fermion) CFT. 
The scaling dimensions are = 1/2, and ~ 1/16. 
Explicit conformal blocks for this case will be quoted later 
in the paper. The fusion rules for this CFT (other than 
for products with the identity 1) are 

X V = 1, (1-20) 
ip X a = a, (1.21) 
a X a = l + -)p. (1.22) 

The fusion rules imply that the number of conformal 
blocks for n quasiholes on the sphere is 2"/^^^, and n 
must be even. Other examples include the RR states 
[20l | , in which the field ijj = ijji, one of the parafermion 
currents in the Zk parafermion CFT; these fields are good 
examples of "simple currents" . One can also consider QH 
systems in which the particles carry SU(2) spin greater 
than 0, and then inner products involve the inner product 
of spin states as well as spatial integrals. 

We may now define the overlap integrals, or inner prod- 
ucts of the trial states with wavefunctions ^'a, as 

(*q(wi, . . . , w„)|*b(wi, . . . ,u;„)) 

= Zabiwi, . . . , W„) 

r N 

= / ]^(i^Zi^'a(wi, . . . ; zi, . . .) 

i—l 

x^b{wi,...;zi,...). (1.23) 



(The coordinates, and the integration measure are writ- 
ten for the plane, though the definition applies to other 
geometries with modifications that are hopefully ob- 
vious. We do not imply that the overlap matrices 
Zabiwi, . . . , Wn) are holomorphic in the ws.) The differ- 
ence between Zab and Z above should be noted: both are 
sesquilinear in conformal blocks, but Z{wi, . . . , zi, . . .) 
depends on z's as well as w's, and is a diagonal sum 
of mod-square conformal blocks, while Zab is integrated 
over z's, but not summed over the indices a, h. 

From the discussion we now see that when conformal 
blocks are used as (or in) trial wavefunctions as described 
above, then the holonomy under adiabatic transport will 
equal the monodromy provided that the overlap matrix 
Zabiwi, . . . ,Wn) is proportional to 5ab with proportion- 
ality constant independent of the positions w, asymptot- 
ically for large separations of the w's, in a basis for the 
conformal blocks in which the braiding matrices B are 
unitary. This orthonormality of the conformal blocks is 
then the desired statement generalizing screening in the 
Coulomb plasma. (MR noted in Abelian examples such 
as the Laughlin states that when conformal blocks are 
used as trial wavefunctions, the gauge is such that the 
Berry connection vanishes, and the holonomy is given 
entirely by the monodromy.) 

Let us also point out here that the so-called shift for 
the ground state on the sphere or disk geometries can 
be obtained from the CFT as well. The shift S in the 
number of flux piercing the sphere is defined by 

N^ = v-^N -S. (1.24) 

The flux can be obtained from the degree of the 
wavefunction in each coordinate z^, which itself can be 
obtained by letting z^ — > c» and extracting the leading 
power of Zi (neglecting the Gaussian factor). The CFT 
contributes —2h^ to this, as particle i and the TV — 1 re- 
maining particles must fuse to give the identity [20| . The 
charge sector contributes i'~^{N — 1). Then the shift is 

S = 2{iy-^/2 + h^). (1.25) 

We write it in this form because when the charge sector is 
interpreted as a CFT also [i| , the conformal weight of the 
fleld contributing to the particle is /2, so that 5/2 is 
simply the total conformal weight of the field representing 
the particle, including the charge sector. 



D. Bose-Einstein condensates and paired states 

Some similar trial wavefunctions also have applications 
outside of the QH effect. Let us start again with the 
simplest case, that of a Bose-Einstein condensate (BEC), 
with n > vortices included in infinite space with no 
background potential. A trial function for this is similar 
to the Laughlin function divided by its modulus (but with 
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the particle-particle factors completely removed) 

Unlike the Laughlin states, in such a condensate, while 
the average charge (or number) density is uniform, there 
are large (Poissonian) fluctuations in the density, or in 
the number in a subregion. This is connected with the in- 
finite compressibility of this BEC in non-interacting par- 
ticles. One may wonder if the vortices possess fractional 
statistics. It is simple to perform the adiabatic calcula- 
tion, using expressions similar to those in Arovas et al. 
0. There is a part of the holonomy phase factor re- 
lated to the expectation of the charge density times the 
area of the loop enclosed by the path, which is related to 
the Magnus force on the vortex. When vortices are ex- 
changed, there is a correction due to the charge deficiency 
at the vortex. However, for this trial function, the latter 
is clearly zero, as the vortices disappear from the den- 
sity calculation when the mod-square is taken. Thus, the 
vortices are bosons. The screening effect of the Coulomb 
plasma that arose on taking the mod-square, which pro- 
duced a net deficit of charge around each quasihole, which 
was so important in the discussion of the QH functions, 
is simply absent here. 

However, for a BEC, the trial wavefunction above is 
not very physical. For one, it has an unpleasant sin- 
gularity at the locations of the vortices. Even for non- 
interacting particles, it does not solve the Schrodinger 
equation (a function that does is nii(-^* ~ wi), which 
is not normalizable; the situation is better for trapped 
atoms in e.g. a harmonic potential, when the LLL can 
again be used). Further, in an interacting Bose super- 
fluid, the circulation of the fluid around the vortex pro- 
duces a "centrifugal" force effect, and due to the finite 
compressibility of the fluid there is a long-range tail in 
the deficiency of density compared with the background, 
going as ~ f^/''^, where r is the distance from a vor- 
tex, and ^ is the healing length (this result may also be 
obtained from a Gross-Pitaevskii equation analysis — see 
Ref. [2l| — and shows that the simple form above is not 
even valid asymptically outside the vortices) . As pointed 
out by Haldane and Wu this leads to a charge 

deficit in a circle of radius R centered at the vortex that 
increases logarithmically with R. Thus, the fractional 
statistics phase is not path independent; it depends on 
the separation of the vortices. (Compared with the QH 
case, it also depends differently on the signs of the vortic- 
ities, since the deficiency of particle number at a vortex is 
independent of the sign of its vorticity; the above wave- 
function represents vortices all with positive vorticity.) 
Clearly, this is connected with being in a phase of matter 
that is not "topological" , due to the existence of gapless 
Goldstone (density) modes of the superfluid, as required 
by the broken symmetry. Ref. also pointed out that if 
the interaction between the particles falls off slower than 
(in which case a neutralizing background potential 
will be required), then the fluid exhibits screening, and 



as for the trial function above, there is no net charge ac- 
cumulated at the vortex, and the net phase for exchange 
of two vortices is zero (in the borderline case of 1 /r^ in- 
teraction, a non-zero result is possible). 

In this paper we will consider BCS paired states of 
fermions in addition to QH states. In these, the vor- 
tices carry vorticity in multiples of half the usual unit, 
due to the pairing. If the fermions are not coupled to 
any gauge field (either the electromagnetic field, or the 
Chern-Simons field that arises in composite particle the- 
ory [23] ) , then the fermion wavefunctions must be single- 
valued even in the presence of vortices, while the local gap 
function (or condensate wavefunction, or pairing func- 
tion) must wind in phase by 2tt on making circuit around 
a minimum- vorticity vortex. The charge sector contribu- 
tion to the adiabatic statistics is expected to come from 
viewing the system as made of composite pairs of par- 
ticles, and the pairs behave as bosons, similar to the 
BEC. The point we wish to emphasize Q is that, like 
the BEC wavefunction above, the "nice" trial wavefunc- 
tions that will be considered here do not include the re- 
sult of the self-consistent calculation of the gap. That 
is, the gap function for the pairing should be calculated 
by solving self-consistency conditions from BCS mean- 
field theory that incorporates the presence of the vortices. 
This is difficult, and in general these details should not 
be relevant to the topological properties of a topological 
phase. The result of such a calculation should be simi- 
lar to the hydrodynamic or Gross-Pitaevskii calculation 
for the Bose superfluid, and (for short range interaction 
of the fermions) the density deficiency at the vortices, 
and consequently the charge sector contribution to the 
adiabatic statistics, will not be well defined. Again, the 
neutral paired fermion superfluid is a gapless phase with 
a Goldstone mode, if there is no long-range interaction. 
Thus, while there may be a non-Abelian contribution to 
statistics from the paired wavefunctions, as discussed in 
Ref. Q and the present paper, the Abelian contribution 
is not well-defined; in this sense, these systems are not 
in a topological phase. The same applies to other gap- 
less degrees of freedom, for example when spin-rotation 
symmetry is broken. This should be kept in mind when 
considering the use of non-QH paired systems for topo- 
logical quantum computation, such as the half-flux vor- 
tices in He'^ 24] . These will most likely only be successful 
if the Abelian phases drop out of computations. (Even 
in QH systems, the Aharonov-Bohm phase is an incon- 
venience.) 

For interactions falling slower than with a neu- 

tralizing background added, the contribution to adiabatic 
statistics from the charge sector is well-defined, but is 
zero. Thus for interacting electrons with 1/r Coulomb 
interactions, only the effects other than the charge sec- 
tor will be left. These can be calculated using the trial 
wavefunctions, if there are no other Goldstone modes in 
the system apart from the charge mode. 

Returning to trial paired states in the presence of vor- 
tices but without the self-consistent calculation of the 
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gap function, one can make a singular gauge transforma- 
tion as a function of the fermion coordinates that turns 
the trial state into one with a gap function that does 
not wind in phase, with the many-particle wavefunction 
for the fcrmions changing in sign on making a circuit 
around the vortex. If all the vortices have positive vor- 
ticity, the transformation is multiplication by the inverse 
square root of the above BEC trial function. Then for 
this second gauge choice, it was shown in Ref. Q that the 
long-distance behavior oi a p + ip paired state of spinless 
or spin-polarized fermions is the same as that of the MR 
trial state, as given by the Ising conformal block (or Pfaf- 
fian) above. More generally, we might view conformal 
blocks, now without the charge sector factor 'J charge, as 
trial wavefunctions for Abelian anyons in zero magnetic 
field. These represent superfluid states of anyons, and we 
will include this possibility in the following discussion. 
Wc should point out that this relation of QH wavefunc- 
tions to those for particles (possibly of different, though 
still Abelian, statistics) in zero magnetic field is precisely 
the idea of "composite particles" [23| , here specialized to 
composite particles in zero net magnetic field. Later we 
will argue when performing the calculations of overlaps 
for the QH trial functions that the charge part may as 
well be removed, implying that the results also have ap- 
plication to particle systems in zero magnetic field, con- 
sistent with the composite particle point of view. How- 
ever, as we have just seen for the boson case (which is 
related to composite bosons W§\ ) , in practice there can be 
important differences in the behavior in the charge sector 
that distinguishes these two types of physical systems. 

For such trial functions (conformal blocks) that are not 
single-valued in the particle coordinates, some technical 
issues must be dealt with in order to discuss the adia- 
batic statistics of the vortices. The monodromy of the 
functions is well-defined if one keeps track of all the zs 
and ws. But for our purposes we wish to move w only, 
with z's fixed, and the result depends on the precise path 
taken relative to the z's, due to the square roots in the 
trial function (this generalizes directly to particles that 
are Abelian anyons) . As we wish to compare with holon- 
omy calculated by integrating out the particles in each 
infinitesimal time step, this dependence on particle po- 
sitions is unacceptable, even though it leads only to an 
ambiguity in sign in the present case (more generally, 
to some root of unity). This effect was not present for 
the QH trial wavefunctions above, which included the 
charge sector and were single valued in the particle co- 
ordinates (and the particles were bosons or fermions). 
[We note that the Berry connection is well-defined even 
for the non-single valued wavefunctions, because the de- 
pendence on the particle coordinates is just a square (or 
other) root, and the phase change cancels in the overlap 
(thus relies on the particles obeying Abehan statistics).] 
One solution when the charge sector is removed from the 
QH wavefunctions to obtain the blocks considered here is 
to retain the gauge-transformation phase factor so that 
the functions are single valued in the particle coordinates. 



which amounts to using the original gauge. Then we ex- 
pect still to obtain the charge sector contribution that 
is the area enclosed by the loop (because the particle 
density is uniform), in addition to the contribution of 
the conformal blocks. Another simple solution, which we 
adopt here, is that when the exchange of vortices is made 
along a path C that is contractible to a limit point on the 
intersection of some diagonals, then we may define it as 
made along a different path, homotopic to the original in 
the vortex (or quasihole) configuration space, that does 
not enclose any particles. This is acceptable for the mon- 
odromy as it does not require well-separated vortices (or 
quasiholes). This approach is not available in the case 
of exchange by non-contractible paths, for example on a 
surface of non-trivial topology, but we will not enter into 
this in this paper. 

The preceding applies to fermions that either are not 
coupled to a gauge field, or are but the penetration depth 
for the gauge field is large. If the penetration depth is in- 
stead small, and we consider exchange of vortices at sep- 
arations larger than the penetration depth, then there is 
a circulating pure-gauge vector potential outside a pen- 
etration depth from the vortices, and the gap function 
is covariantly constant. This corresponds to the use of 
the second gauge choice above. The calculation of adi- 
abatic statistics may be made well defined by the gauge 
transformation technique as described above. 



E. Earlier work 

There have been various earlier steps towards demon- 
strating that non-Abelian adiabatic statistics occurs in 
trial QH wavefunctions based on conformal blocks, and 
in certain BCS paired states. The idea that in the MR 
wavefunction, the holonomy equals the monodromy was 
re-emphasized (though not using this terminology) by 
Nayak and Wilczek (NW) \23\ , who also emphasized that 
this generalizes screening in the plasma mapping for the 
Laughlin states. They also found explicitly the two con- 
formal blocks corresponding to any even number N of 
particles and n = 4 quasiholes. Even though they did 
not find these in the form of Majorana fermion zero mode 
states on the quasiholes (which was found at around the 
same time in Ref. [11]), they guessed that this interpre- 
tation was correct for any number of quasiholes. This 
led them to conjecture the form of the braid group repre- 
sentation in the monodromy, which apart from Abelian 
factors (i.e. tensor product with an Abelian representa- 
tion of the braid group) can be viewed as an image of the 
braid group in the spinor representation of the rotation 
group in n dimensions [25*]. This representation of the 
braid group was known [27] , and was also known to oc- 
cur in the Ising_^CFT (its structure is described in more 
detail in Ref. j28|). A similar argument was spelled out in 
greater detail in Ref. jlJ] , after the work of Ref. [4] on the 
Majorana zero modes in p + ip paired states. However, 
the argument gives only the monodromy of the states 
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(modulo Abelian factors) , and it is not clear if adiabatic 
transport is actually considered in Ref. [11] (no expres- 
sion for the Berry connection appears there). Another 
argument of NW is somewhat similar to the one we will 
give in Sec. fVl below. 

Adiabatic transport of quasiholes or vortices in the 
paired state was considered further more recently in Ref. 
[29| |. and especially clearly in Ref. [13] (Appendix A), 
where it is shown that the Berry connection is propor- 
tional to the identity matrix, thus proving that the holon- 
omy is given by the monodromy found by NW, up to 
some Abelian factor. 

Other approaches to the problem for the MR QH state 
should also be mentioned. Gurarie and Nayak [3l|] used 
another Coulomb gas method from CFT to represent the 
overlap integrals. For the case of only two quasiholes, 
they succeeded in obtaining the vanishing of the Berry 
connection, and hence that the holonomy equals the mon- 
odromy in this case. For four quasiholes their result de- 
pended on some assumptions, the validity of which does 
not appear to be obvious. Other groups [H, [s^ for- 
mulated field-theoretic arguments, but seem to assume 
that the ed ge t heory is the expected CFT. Tserkovnyak 
and Simon [34| evaluated the holonomy numerically for 
two and four quasiholes by Monte Carlo methods, find- 
ing agreement with the expected result, at some degree 
of accuracy. 

For most other states, such as those of Read-Rezayi 
[20| . much less has been shown. But there is a series of 
spin-singlet trial states due to Blok and Wen f35| for par- 
ticles of SU(2) spin fc/2 (fc = 1, 2, . . . ; fc = 1 is Halperin's 
Abelian spin-singlet state, see e.g. Ref. [l|), in which the 
CFT is SU(2) level k, which have many nice properties. 
These authors were able to show that the Berry con- 
nection vanishes, and so holonomy equals monodromy, 
for these states by u sing the Knizhnik-Zamolodchikov 
equation from CFT [l7[, and making an assumption 
that screening holds in an SU(2) generalization of the 
Coulomb plasma, in a manner closely parallel to the work 
of Arovas et al. @. The screening assumption implies 
that the trial ground state has short-range spin correla- 
tions. We will comment on this further in Sec. fVl below. 
For recent further results on monodromy of blocks in the 
RR series, see Refs. ^36i.i37i]. 

We conclude from this survey that with few exceptions 
existing results in the literature are either only partial 
ones for the MR state (not calculating the Abelian fac- 
tors), or else depend on unproven assumptions, or apply 
only to particular states. By contrast, the results pre- 
sented below for the MR state are complete in that they 
yield the full holonomy for up to four quasiholes, or for 
the ground states on a surface of any genus. An argument 
that we invoke frequently, which is supported by renor- 
malization group arguments herein, is that the charge 
sector factor of the QH wavefunctions can be dropped 
without jeopardizing the results for holonomy in the CFT 
sector; this is related to conventional lore about compos- 
ite particle methods, and has of course also been used by 



others, for examples, Refs. [l^,!!^]. The results for quasi- 
holes rest on an assumption, that screening occurs in a 
certain very conventional two-component plasma, which 
will be accepted by most physicists. The arguments given 
in Sec.lVlapplv to any trial state given in terms of confor- 
mal blocks as explained here, and show that holonomy 
equals monodromy under some general conditions that 
can be checked for each particular trial state as a well- 
posed physical question in two-dimensional field theory. 
Further, there is a simple easily-checked criterion (rele- 
vance or irrelevance of a perturbation) that may provide 
important clues as to whether or not the general condi- 
tions hold. 

Now we turn to earlier work on adiabatic variation 
of the aspect ratio of a QH system on a torus. In the 
complex plane, the torus is defined by identifying points 
under z — > z -|- L and z — > 2; -I- Lt, where r is in the 
upper half plane Imr > 0; thus the upper half plane is 
the parameter space on which we may study adiabatic 
transport. This was considered for the filled LLL in an 
elegant paper by Avron et al. [^. They showed that 
there is a contribution that is not holonomy, but curva- 
ture (anholonomy or field strength) of the Berry connec- 
tion (vector potential) on the upper half plane. (This is 
somewhat analogous to the Aharonov-Bohm phase that 
is proportional to the area enclosed by a loop, in the 
Arovas et al. calculation for moving one quasihole.) 
For the case of the filled LLL, the result is proportional 
to the total flux through the torus [3]. Physically, for 
a homogeneous many-particle fluid state, this adiabatic 
curvature divided by the system area represents the Hall 
viscosity, here denoted 77^^^ a non-dissipative transport 
coefficient (38i . i39,] that is known in plasma physics but 
often overlooked elsewhere in fluid dynamics; 77*^ '^^ is the 
only coefficient of viscosity that can be non-zero in a 
two-dimensional isotropic incompressible fluid (but must 
vanish if time-reversal symmetry is present) [^. The re- 
sult for the filled LLL can also be extracted from the 
detailed single-particle results of Levay [i^]; this paper 
clarifies many aspects of this problem. More generally, 
for the state in which the lowest v Landau levels are filled 
{v integer), the result r/"^-* = fivn/A is quoted for inte- 
ger V, \v\ > 1 in Ref. [33, again using results from Ref. 
(40] {n is the particle density n = ]z^]/(27ri!^); h and the 
magnetic length £b are set to 1 elsewhere in this paper) . 

Note that following comments from the authors of Ref. 

3, we have corrected the coefficient to 1/4 to account 
for apparent typos in Ref. [38j.] 

A recent paper ^ has argued that the result of Ref. 
generalizes to an arbitrary QH state in the LLL to 
give 77^"*' = Tiv / {^TTpg) for arbitrary (possibly negative) 
values of v, independent of the state (this paper [7[ , and 
references therein, use the Hall viscosity in an interest- 
ing hydrodynamic approach to collective modes in QH 
fluids). Unfortunately, the end of the calculation uses 
the incorrect argument that all particles are close to the 
X-axis. The single-particle results of Levay [40] can be 
applied to obtain a seemingly similar result. He finds. 
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for the Berry calculation in the iV^-fold degenerate space 
of single-particle states in a single Landau level {N^ is the 
number of flux piercing the torus), that the Berry con- 
nection and its curvature are proportional to the identity 
matrix in this space, with a coefficient related to the ki- 
netic energy of the Landau level. This implies (e.g. by 
using the Slater determinant basis) that for the space of 
all many-particle states with all the particles in a sin- 
gle Landau level, the Berry connection and curvature are 
again proportional to the identity matrix, and so in this 
sense are independent of the state. (This generalizes fur- 
ther if we consider a space consisting of all the many- 
particle states with a given number of particles in each 
Landau level.) However, this adiabatic transport is not 
what should be considered for a quantum fluid state, in 
which because of the presence of interactions we take a 
family consisting of a single ground state (as considered 
in Ref. Q) for each value of t, or more generally a space of 
"degenerate" states for each r, but generally not the full 
space of many-particle single Landau level states. This is 
an example of the general set-up for adiabatic transport, 
in which a family of subspaces within a common Hilbert 
space is considered, as discussed above in Sec. II Bl (The 
T-dependent Landau level states themselves arise in this 
way, as all are subspaces of the space of square- integrable 
hmctions on the torus with given boundary conditions.) 
The curvature of the Berry connection definitely does 
in general depend on the choice of a subspace for each 
value of T (though not on the gauge, that is the basis 
for the subspace), and cannot be reduced to a calcula- 
tion that ignores restriction to the subspace. Hence the 
result that is independent of the state chosen, as claimed 
in Ref. 01 , cannot be as general as was stated. Later, we 
will show that the correct result depends on the form of 
the ground state, not only on the density, though it is 
universal within a topological phase. (We will also find 
related results for paired supcrfluids. We are not aware 
of any earlier results for Hall viscosity of paired states.) 

From Levay [40| , the result of Ref. [7] can be viewed as 
the correct one for non-interacting particles in the LLL. 
Also, using Levay's results for non-interacting particles 
but at non-zero temperature T, and going to high tem- 
perature using classical equipartition of e nerg y, we find 
that the result agrees with that in Ref. [39|, which is 
Ty^'^) — kBTfi/ {2uoc), where ujc is the cyclotron frequency. 



F. Structure of paper 

In Sec. ini we first review the essentials of paired states 
at the BCS mean-field level, then specialize to ground 
states on the torus. We calculate the normalization 
factors, then examine the monodromy under modular 
transformations, and then the curvature and holonomy 
of the Berry connection for changes in the aspect ratio 
r; this determines the Hall viscosity of BCS paired sys- 
tems, for which we also give a simple direct calculation. 
We also discuss higher genus surfaces, and the strong- 



pairing phases. In Sec. IIIIl we examine similar questions 
for the Laughlin QH states, relate the Hall viscosity for 
trial states given by conformal blocks to the conformal 
weight of the field for the particle, and conclude the gen- 
eral discussion for Hall viscosity. In Sec. IIVI we present 
direct arguments for the non-Abelian adiabatic statistics 
of two and four quasiholes in the MR state on the sphere 
or plane in the thermodynamic limit. The calculations 
work by "doubling" (taking two copies) of the system 
with the charge part removed, and using an argument 
that a plasma is in a screening phase. In Sec. |Vl we 
present general arguments that amount to necessary and 
sufficient conditions for trial wavefunctions given by con- 
formal blocks to describe a topological phase with adi- 
abatic statistics given by the monodromy of the blocks. 
The condition is that related correlation functions in two 
dimensions should go to a constant, as in a ordered phase. 
We discuss numerous examples in this light. Cases not 
obeying the condition are argued to be gapless phases or 
critical points. We argue that use of non-unitary ratio- 
nal CFTs in our way cannot produce such a topological 
phase if there are any negative quantum dimensions in 
the theory (which follows if there are any negative con- 
formal weights) . The argument assumes that the twist in 
the theory, defined adiabatically, is also the same as in 
the CFT, which has not been shown. In an Appendix we 
discuss this, which is the last step in deriving a modular 
tensor category from the construction, and show that the 
consistent twists are very limited, so that the argument 
does go through in at least one family of examples. 



II. PAIRED STATES IN A CLOSED FINITE 
SYSTEM 



In this section, we consider the p + ip paired states on 
compact two-dimensional surfaces (no boundary). The 
set-up was already described in Ref. ^] (see also Ref. 
[sot), but we review some essential steps, and add a few 
points. The basic case, other than the sphere, is the 
torus, or equivalently periodic boundary conditions on 
a parallelogram in the plane. The other cases, surfaces 
of genus higher than one, require more work to set up, 
and we will be more brief. For all these problems, we 
can show starting from the general pairing problem, that 
at long wavelengths the orthonormalized wavefunctions 
(within the BCS mean-field formulation) are the confor- 
mal blocks of the Majorana fermion (Ising) CFT. This 
allows us to calculate explicitly the adiabatic transport 
of the states as the aspect ratio of the torus is varied; this 
leads to the Hall viscosity and the modular transforma- 
tion group. We also consider the generalization to other 
paired phases of fermions, including the strong-pairing 
phases. 
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A. General equations for pairing 

We begin with BCS mean field theory [4lj. The effec- 
tive grand-canonical Hamiltonian for the fermions is, in 
the most general case (we assume for the present that the 
fermions are in a finite set of orbitals labeled by Greek 
indices like a, taking values a = I, . . . , M) 



(2.1) 

Here the creation and annihilation operators, c\ and c^, 
which are adjoints of one another, obey the canonical 
anticommutation relations {ca,cjj} = 5aj3, {ca,C/3} = 0. 
The M X M matrix h must be Hermitian, hap = h/Sa 
(the "is complex conjugation), and represents the usual 
kinetic energy and any one-body potential terms, while 
A is an antisymmetric M x M matrix, A^^ = — A^gct and 
corresponds to the gap function. The former condition 
ensures that Kcs is Hermitian, while the latter means 
that there are no redundant components of A. 

It will be convenient to write everything in terms of 
a 2M-component column vector of field operators C, in 
which the first M components are Ca, the remaining AI 
are cj,. Then, up to an additive constant. 



I , h A 
2^ I -A -h 



C 



-c^nc. 



(2.2) 



The 2M-dimensional vector space can be viewed as 
a tensor product of A/-dimensional space with a two- 
dimensional "particle-hole" space. Then we may use a 
tensor product notation for matrices, with Pauli matri- 
ces cr^, fi = X, z, acting in the two-dimensional space. 
Then the 2M x 2M matrix H obeys 



(2.3) 



where is transpose, Yi^ = <Jx ® Ia/, and \m is the 
M X M identity matrix. 

We will now diagonalize and find the ground state, 
by performing a Bogoliubov transformation. First we 
note the anticommutation relations for C, in terms of a, 
/3=1, ...,2M, 



{Cq,C^} — 



(2.4) 



To diagonalize the Hamiltonian, we require creation and 
annihilation operators for quasiparticle modes, Xr, xl 
{r — 1, . . . , M), which we combine into a 2M-component 
column vector x, similar to C. These operators x must 
obey the canonical anticommutation relations of the same 
form as those for C above. If we define explicitly 



Ux 



u V 
V u 



X 



(2.5) 



where u, i', are complex M x M matrices, then requir- 
ing C and X to satisfy the canonical anticommutation 
relations, we must have 



0, 

1m. 



(2.6) 



Given the above form for J7, these equations are equiv- 
alent to UU^ — I2M, and it follows that J7~^ = JJ^ = 
'Y^xU^'^x- These imply that U is an element of the Lie 
group 0(2M). We note that the condition U^U = I2M 
leads to additional relations 



u^v + u — 0, 



t„ - 



LA/, 



(2.7) 



which may be useful, but are equivalent to those above 
as we are dealing here with finite matrices for which the 
right and left inverses of a matrix are equal. 
To make K^s diagonal, we require that 



off 



E 








X = ^X^Ex, 



(2., 



where E is an M x M diagonal real matrix, i? is a 2AI x 
2M matrix, and Kq is the ground state energy, where by 
"ground state" we mean the state annihilated by all Xr- 
The conditions for E to be diagonal, which can be found 
by commuting Xr with K^g, are the matrix equations 
(a complete eigenvalue problem) HU = UE, or more 
explicitly 



hu + Av = uE, 
Au + hv = ^vE, 



(2.9) 



which are the BdG equations for this problem. (The 
BdG equations may look more familiar if one takes the 
rth column of each of these, which gives the eigenvalue 
problem for the rth eigenvalue Er, ior r = 1, . . . , M .) We 
notice that the equations for u, TJ, which are the complex 
conjugates of these, are in fact the same, by replacing 
u V, V ^ u, E —E. This can also be applied to 
an individual column of u, v, say the rth, to effectively 
exchange the corresponding Er with —Er- 

The BdG equations are not easy to solve in closed form 
in this general case. Fortunately we do not require these 
general solutions. For the moment, all we need is the 
form of the ground state wavefunction. We will concen- 
trate on the case in which the ground state contains only 
states of even particle number. The ground state can be 
written in the unnormalized BCS form 



= exp 



a,/3 



|0), 



(2.10) 



where |0) is the vacuum, which is annihilated by all c^s, 
and g is an antisymmetric M x M matrix. The condition 
that all Xr annihilate yields the relation 



T 

u g 



(2.11) 
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as a matrix equation, and so g = vu~^ by taking the 
transpose, or by using the first of eqs. (|2.7p . The wave- 
function of the ground state is then a Pfafhan, up to an 
overall sign: 

^'(ai, . . . , QfAf) = {0\Cai ■ ■ ■ Ca„ \^) = ±Pf (2.12) 

Note that for any state jfi), with wavefunctions 'J defined 
in this manner for each iV-particle subspace, the norm 
square is 

("l^>=E E I*(«i,---,«A^)IV^!- (2.13) 

Ql,...,QiV 

We now obtain the normalization of the ground state 
wavefunction. To calculate we expand the expo- 

nentials and use Wick's theorem. With the help of some 
simple combinatorics, we obtain 

= cxp ^Tr ln(lA/ + g) - ^Jdct{lM + g^ g)- 

(2.14) 

Using the second of eqs. (|2.7p . we obtain 

uu^ ^{Ui + g^g)-'. (2.15) 

Hence {det uu^Y^^\fl) is normalized, a generalization of 
the usual normalization for BCS theory in momentum 
space (see e.g. Refs. [3, !4l!]). Also the expectation value 
N of the particle number N — ^^cl^Ca is 

N = Ttv^v ^ Trg^g{lM + g^gy^. (2.16) 

By writing the BdG equations in terms oi E = uEu^^ 
in place of one finds that 

g = -K'\E-h), (2.17) 

and then by eliminating E, g obeys a matrix quadratic 
equation, 

A + gh + hg + gAg = 0, (2.18) 

so that g can be found without first solving for E or E. 
(This is similar to a diagonal equation in Ref. p. 
68.) This equation may also be obtained more quickly 
by applying Kcs to jfi), whence one also finds that the 
ground state energy is Kq = —TTAg/2. 

In some situations of interest, h is a constant multiple 
of the identity matrix, say h = In this case the BdG 
equations (j2.9p can be analyzed further. By multiplying 
the first of equations (|2.9p on the right by {E — ^Im), 
and the second on the left by A, then eliminating v, we 
obtain 

uE^ ^ {^J.'^lM + A^ A)u. (2.19) 

Here we used also —A = A'^, and we note that A^^A is 
Hermitian and positive. Thus we see that the columns of 
u are the eigenvectors of A^^ A. Similarly, the columns of v 
are the eigenvectors of A A^^ . It follows that the diagonal 



entries of E"^ are greater than or equal to /i^. In fact, 
a standard result from linear algebra allows a complex 
antisymmetric matrix A to be written in the form 

A = w^Dw, (2.20) 

where w is unitary, and D is an antisymmetric AI x M 
block-diagonal matrix of [M/2] {[x] means the largest 
integer < a;) 2 x 2 blocks of the form Xa(Jy, where Aq 
(a = 1, . . . , [M/2]) obey Aa > 0, and for M odd also 
a single 1x1 block containing 0. The transformation 
can be considered as a basis change among the orbitals 
C only, or as a Bogoliubov transformation with u ~ w 
and w = 0, that brings 7i to a form with A replaced by 
D, so that the gap function links orbitals in pairs only 
(similar to the case of BCS theory for a translationally- 
invariant system in k space, where the gap functions links 
k and — k only). It follows that the eigenvalues in E^ are 
ji^ + with multiplicity 2 for each a, and an additional 
^2 if M is odd. 

Taking the positive solution for E^ we have 

E^w^/n^M + A'^Aw-^, (2.21) 

and the positive square root of a positive Hermitian ma- 
trix can be defined through spectral theory. At the same 
time, we also have 

E = u-^Eu = Vm^Im + At Am, (2.22) 
and consequently g is 

g = A~' (^Im + Vm'Im + At a) , (2.23) 

which solves the quadratic equation for the present case, 
which is 

A - 2fig + gAg = 0. (2.24) 
This solution leads to 

5^3= (m+ TM^TAtA) (AtA)-i (m+ V^^^TAtA) 

= (m+ V/^2 + AtA)'(AtA)-\ (2.25) 

The two expressions for E are consistent, because we 
can write the polar decomposition of u, u — \u\w'~^, 

where |u| is defined as \u\ = VuvJ = {1m + 9^g)^^^'^ 
and w' is unitary. Then we see that |u| commutes with 

II^Im + At A, and we can take w' = w. As E is posi- 
tive, it follows that the solution for g describes the lowest- 
energy ground state. 

Next we will pass to a continuum limit, in which 
the finite matrices are now operators as in an infinite- 
dimensional single particle Hilbert space. This applies 
to systems in continuous space, and also for the large 
distance limit of a lattice. In all cases, our interest will 
be in scales larger than the coherence length, which can 
be thought of as the characteristic size of a Cooper pair 
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(and which we will assume does not depend significantly 
on position). Then the terms in Ti. will be assumed, first 
for the translationally-invariant case, to be of the form 
h = k"^ /{2m*) — and A — A{kx — iky) when written 
in momentum space; at present m* (the effective mass 
of the fermions) , /i (the chemical potential) , and A (the 
magnitude of the gap function) are constant. We have 
assumed for A the characteristic form for — ipy pair- 
ing, as in Ref. [4] (the coefficients of and py are the 
same, but a different non-zero ratio could be removed 
by a rescaling once and for all). Then in position space 
these become differential operators; to represent these as 
functions of two positions ri, r2 as in the above formal- 
ism, with a replaced by r = ix,y) everywhere, one would 
consider them as integral kernels, by including suitable 
(5- functions. They are 



h = -VV(2to*) - A* 
A = -2iAd 



(2.26) 



where d = d/dz = \{dx — idy). In fact, the k^ term in 
h can also be dropped for our purposes, even at = 
which is the transition point between weak and strong 
pairing . The behavior of both operators is implicitly 
assumed to change at large wavevectors, so as to include 
suitable physical cutoffs, including the coherence length. 
For long-wavelength purposes, the above results for h = 
—fi a constant are now applicable. 



B. Periodic boundary condition ground states 

The case of periodic boundary conditions was already 
considered in Ref. [3], so we will be brief, and just make 
some additional points. For h = —fi, A both constant, 
we can apply the above results in position space. From 
here on, we use units in which |A| = 1 in this case. Then 
as g acts from the left on , use of the standard Fourier 
transform definitions, and by comparison with 

Keff = -/icj^Ck + ^ (A^C_kCk + Ak^cli^ 



(2.27) 

we obtain Ak ^ k^ — iky, and 

5k = (a* + + |k|2) /{k^ + iky), (2.28) 

which agrees with Ref. up to an unimportant minus 
sign [1]. In position space, for > 0, we find g(ri, r2) = 
g{ri - r2) 2n/{zi - Z2) for |ri - r2| > 1/1^^1, where 
zi = xi +iyi, etc. For distances less than of order 
the behavior changes 

This behavior of g{r), which is what appears in the 
many-particle wavefunction (the Pfaffian), might sur- 
prise us as we usually think of the fermions as forming 
bound Cooper pairs, whereas the integral of our |.9(r)p 
over infinite space diverges. However, the usual physi- 
cal way to think of these pairs is to examine the order- 
parameter function (c(ri)c(r2)), which is a function of 



ri — r2 only. This is the Fourier transform of (ckC_k) = 
'"k'^^k = -Ak/i?k (where Ey, = y/ fi'^ + |kp). Using the 
same approximations, this is analytic near k = 0, and so 
has no long-range tail at distances > 0(l/|/i|), though 
it does at shorter distances. If = 0, one finds the 
same behavior ^ l/(z|z|) as \z\ — > 00 (as for g{r) in this 
case). This function is normalizable, and so the size of the 
Cooper pair is still the coherence length ^0 l/l/^l that 
must be defined using shorter-distance or larger wave- 
length behavior of the pairing (note, however, that these 
arguments show that if one attempts to use the expecta- 
tion of |rp in this function to define the size of a pair, one 
encounters a logarithmic divergence at a* = 0). Then in 
the topological (weak-pairing, fJ. > 0) phase of interest, 
the interesting effects will be transmitted by the long- 
range behavior of g{r), yet this does not show up if one 
looks at the natural order parameter as above. This is to 
be expected for a topological property of a phase. 

We now consider the ground states on the torus, as 
this is a warm-up problem for the issues of most interest 
in this paper. Here we can take the torus to be the r 
plane, modulo translations by two linearly-independent 
vectors lji, u)2, with u>i along the positive a;-axis and 
of length Lx- If we identify these vectors with complex 
numbers, then the first becomes Lx, and we define the 
second to be L^t. We assume that Imr > 0, and define 
Ly = Xj^Imr. When we vary the aspect ratio r we will 
assume that L^Ly is fixed. First we assume that the 
fermions obey a periodic boundary condition c(r-|-a;i) = 
c(r + u)2) = c(r), and that there is no vector potential. 
Then from the underlying physics that produces p + ip 
pairing, we can assume that there is a ground state with 
A constant. In the weak-pairing phase /i > 0, this has 
odd particle numbers N in all components of the state 
By boosting all the fermions by a momentum of 
order 2'k/Lx, we can imagine that there is also a low 
energy state with a net circulation around the x cycle 
parallel to Wi, or similarly for the y cycle parallel to 0)2, 
or both. In terms of the small wavevector behavior of 
the gap function A, these involve A(r) that winds by 2tt 
along the x cycle, or the y cycle, or both, respectively. 
Because the change in the state is small, locally, these 
should exist as low energy states. 

In the case when A(r) varies with position, some gen- 
eralization of the above expressions is necessary. The 
most important possible variation is in the phase of A(r), 
though the magnitude also varies in some important sit- 
uations. Then in either case it becomes 



A 



-i{A(r),9} 



(2.29) 



which involves the anticommutator of one-body opera- 
tors. Thus the effective Hamiltonian becomes 

K^s = j d^r c^(r)/ic(r) + i (c^(r)Ac^(r) -t- c(r)Ac(r)) 

(2.30) 

In the present section we can see that A is antisymmetric 
by integrating by parts. We can also see that Kcs is 
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gauge invariant under 

c(r) ^ e*^('-)c(r), 

A(r) ^ e2'^("")A(r), (2.31) 

with no need to include a vector potential. The physi- 
cal motivation for the form of A, at least for A constant 
in magnitude, is to ensure this gauge invariance and an- 
tisymmetry. Thus these forms are valid at long wave- 
lengths whether the superfluid is neutral (no U(l) vector 
potential), or charged (coupled to the scalar and vector 
potentials of electromagnetism, or a Chern-Simons (CS) 
gauge potential in the CS formalism for the FQH effect). 
In the charged or CS cases, the vector potential would 
enter in the higher-derivative terms that we neglect here. 
(We correct here some slight mis-statements about the 
role of the vector potential in Ref. Q-) 

From the form of A = — i{A(r), d}, we can see that in 
the states in which A winds in phase around the torus, 
instead of k pairing with — k, k is paired with — k shifted 
by a small amount. These states correspond to the mini- 
mum quantized amounts of vorticity or circulation of the 
superflow around the two cycles (we do not call these 
flux as there is no vector potential present at the mo- 
ment). They can be generalized further to states with 
TOi, 7712 units of quantized circulation around the two cy- 
cles. These states do not all have the same energies in the 
present situation, because there will be contributions to 
the total energy which can be expressed as a Ginzburg- 
Landau (GL) type functional of A, and this energy den- 
sity will presumably include a term like |VAp (we as- 
sume that the gap function, which after all is actually a 
function in k space, can be adequately characterized by 
a single complex number for this purpose, as in the usual 
GL theory for s-wave paired states). This term arises 
from the extra kinetic energy of the fermions in these 
states, for example. For example, for the case 7712 — 0, 
the total energy cost will go as 7n\Ly/Lx, multiplied by 
constants independent of system size. When rrii, 7712 are 
both even, these states occur at odd particle numbers in 
the weak-pairing phase, while in the other cases (and in 
all cases in the strong-pairing phase /i < 0) they occur at 
even particle numbers. 

An alternative point of view is also useful. As men- 
tioned in Sec. II Dl for each of the winding states, we 
may perform a gauge transformation to bring A back 
to a constant, which will presumably be close in magni- 
tude to that in the (0, 0) ground state. Acting on the 
fermions, this is a phase change by — 7rmi^2 along the 
two cycles, and produces a corresponding change in the 
boundary condition. Thus the boundary conditions on 
the fermions can be described by whether or not there 
is a sign change around each of the two cycles. We will 
write the boundary condition as -I- if the corresponding 
m is even, — if it is odd. Thus there are four possible 
boundary conditions corresponding to the four values of 
(mi, 7712) modulo 2. For each of the four choices, there is 
a tower of states in which the winding numbers relative 



to the lowest energy state in the tower are (mod 2). 
We note that the gauge transformations used here are 
"large" transformations (though not singular anywhere), 
as for mi, m2 not both zero they are not continuously 
connected to the trivial transformation by a path in the 
space of gauge transformations. 

Returning to the original gauge (boundary condition), 
for superfluids in which a U(l) vector potential is present, 
the same winding states still exist, but their energies are 
now shifted because of the possible holonomy of the vec- 
tor potential along the two cycles (again, this effect is 
due to the entrance of the vector potential into the larger 
wavevector aspects of the physics). This holonomy is the 
usual gauge-invariant notion of the flux of the vector po- 
tential threading the two cycles on the torus. Although 
the holonomies for the two cycles can take any real val- 
ues, when they are equal to 7r(mi,m2) for some inte- 
gers mi, m2 (in units of hc/e, in ordinary units) then 
the corresponding winding state has the lowest energy, 
and this energy is independent of (mi,m2) by a (large) 
gauge transformation argument. This is referred to as 
flux quantization [4l[. When the gauge field is viewed 
as an additional quantum degree of freedom, all of these 
states are distinct in the combined Hilbert space of par- 
ticle and electromagnetic field states. Note that physical 
states can be viewed as invariant under gauge transfor- 
mations that are continuously connected to the identity, 
but not under "large" gauge transformations, because the 
holonomies represent real physical flux. 

In the QH applications, the physics is slightly different. 
In this case the vector potential is not a truly independent 
degree of freedom, but instead is tied to the total mo- 
mentum of the fermions. Further this momentum space 
is compactified; it is topologically a torus [i^ . While the 
different winding states corresponding to the four choices 
± for each cycle are still distinct valid ground states [3| 
(when the signs are not both -I- they lie on the bound- 
ary of the Brillouin zone in pseudomomentum [43j). it 
appears that the towers of winding states over each of 
these are not distinct states, but can all be identified 
with the lowest one. 

In the gauge-transformed view of the {mi, 7712) wind- 
ing states, whether a vector potential is present or not, 
the long-distance part of the pairing function g is inde- 
pendent of the vector potential (if present), as we have 
shown, and so the long-distance wavefunction only de- 
pends on mi, m2 (mod 2), that is on the signs ±. The 
pairing functions g in the four cases have been determined 
previously, and in the weak-pairing phase at scales larger 
than 1 /jj, are essentially Jacobi elliptic functions (except 
for some slight complications in the +-|- case) At 
large lengths (i.e. separation |ri — r2| of the particles), 
ga in boundary condition a is 2^ times the inverse of A, 
and thus of d, up to a constant factor independent of the 
system size, aspect ratio, and boundary condition. These 
functions may be determined using complex analysis ar- 
guments, and for mi, m2 not both zero (mod 2) are (up 
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to such a constant) 

3mi,m2(r|T) CX 

'^(mi + l)/2,(m2 + l)/2(^^/-^x|-r)5^1/2,l/2(0|T) 
Lx-dl/2.1/2(zlLx\T)d(„i^+l)/2,(ni2 + l)/2{Q\T) 



(2.32) 



where mi, TO2 = 0, 1, and 'da,b{z\''') ^i'^ elhptic theta 
functions with characteristics, defined by 



'&a,b{z\'i') = 



i7rr(n4-(i) +27r't(n+a)(2:+6) 



(2.33) 



We note that 'di/2.i/2{z\T) is the odd theta function (odd 
under z — )■ —z), and so has a zero at z = and points 
related by translation by 1 or r, while the others we 
use are even about z = (see e.g. [ill). We define 
5i?i/2,i/2(0|r) = 9i?i/2,i/2(^|t)U=o- Thus the elliptic 
functions on the right-hand side of eq. (|2.32p have a sim- 
ple pole with residue 1 at z = 0, and the required (anti- 
)periodicities. For mi, m2 both zero, we have instead 

, I X 5^^1/2,1/2 (^/i^rk) , 2mhxYz/Lx 
Vi/2,i/2{z/Lx\t) Imr 

(with the same proportionality constant as before) which 
is periodic, but not analytic. The non-analytic term can- 
cels in the ground states with odd particle number [3|, 
in which the k = state is occupied by an unpaired 
fermion. That is, the wavefunction for N fermions is the 
product of gs for the pairs, times for the unpaired 

fermion [A = L^Imr is the area), antisymmetrized over 
all fermions. For each pair i, j, the non- analytic term is 
proportional to Im Zi — Im Zj . It breaks into terms Im Zi 
(independent of Zj) and the same with i, j interchanged. 
The former of these two puts j in the constant mode, 
which is already occupied by at least one other parti- 
cle, so this contribution vanishes on antisymmetrization. 
(For a long proof, see Ref. [13] ■) Hence all such terms 
give zero, and the non-analytic term can be dropped for 
all pairs; the wavefunction of course is still periodic [2^ . 

One further step will be useful for later purposes. The 
factor 1/Lx in each pairing function can be rewritten as 
(ImT)^/^/A^/^, where A = L-^lmr is the area. 

Now we turn to the norms of these paired states, in 
the gauge transformed point of view. From the general 
theory above, the factor that is needed in the state to 
normalize it is 



[det(H-5^5)] 
(up to an overall phase) , which equals 

[detAtA]^/' 

det (At A + {fi+ Va^^ + AtA)^^ ^ 



(2.35) 



(2.36) 



In our case, A = —2iD, and A^ = —A = —2iD [here 
D is d, with the understanding that it acts on functions 



of the chosen (anti-)periodicity]. On passing to momen- 
tum space, we find that these expressions are identical to 
those in Ref. [3]. For long wavelengths, the denominator 
is non-singular and we can ignore it; it depends exponen- 
tially on the area of the system, but not on its aspect 
ratio T because effects of finite size are exponentially de- 
caying corrections (formally this is shown by applying 
Euler-Maclaurin methods to the logarithm of the prod- 
uct of factors, and noting the analytic behavior in k near 
k = 0). But the numerator is detD times its conjugate 
(times a constant), and is nonanalytic near k = 0. With 
a momentum cutoff of order /i, this product can be eval- 
uated (for ^ 1/m) ■ The discreteness of k is important 
here, and hence the boundary conditions enter. For the 
three cases other than -|--|-, k = is absent from the 
product, and we consider these first. The logarithm of 
the product is 



\j2^n\k\' = L.Ly I ^lnlklV2 + ..., (2, 



37) 



where the subleading terms represent the corrections due 
to the discreteness of the sum. The first (integral) term 
is again extensive (proportional to the area LxLy) and 
independent of r, and will be discarded; this corresponds 
to renormalizing the determinant. 

At this stage it is helpful to realize that for each of the 
three boundary conditions stated, the product of interest, 
det (£)£)), can be interpreted as the partition function of a 
massless Dirac fermion on a Euclidean torus. (The oper- 
ators D and D are parts of the massless Dirac operator in 
two Euclidean dimensions, see Ref. [3].) The square root 
is the partition function of a massless Majorana fermion, 
which has half as many degrees of freedom. This the- 
ory arises in connection with the two-dimensional Ising 
model on a torus, at criticality. The finite-size correc- 
tions have been computed in that context (Ref. [4^, or 
see e.g. Refs. [13, H^)- To describe these, we introduce 



some functions ¥„ 



of q 



J2'KiT 



that correspond to 



the boundary conditions, m,i.2 = 0, I (mod 2) 

oo 

Fii = g-i/^«n(l + '?"^'^')' 

oo 

Fio = 9"^/^«n(l-'?"^'^')' 

n=0 

Foi = V29i/i6-i/48-Q(i + ^n)_ (2.38) 

These products are convergent for Imr sufficiently large, 
and it is convenient to take the limit of large Lx^i at fixed 

T. 

For the ++ boundary condition, the ground state of 
the paired supcrfluid is at odd fermion number, and there 
is an unpaired fermion occupying the k = state [3] . The 
normalization of the ground state is given by a similar 
product, but with k = omitted. The interesting part is 
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then related to the function 

Foo = V2(Imr) 1/4^1/16-1/48 -Q(i _ ^ny (2.39) 

(The here could be dropped, though not in Fqi.) 

In each of the four cases, the partition function of the 
Majorana ferniion field theory with fixed boundary con- 
dition is \Fm-1.m2 The normalizing factor for our paired 
wavefunctions is \Fm1.m2 1, up to an overall phase and fac- 
tors independent of r and the boundary conditions. We 
now notice that the phase can be chosen so that the nor- 
malizing factors are simply the respective Fmi,m2 them- 
selves, which are holomorphic in t [except for the factor 
(Imr)i/4 in Fqq]. We view these as the chiral partition 
functions, proportional to VdetD (with the zero mode 
deleted for +-I-). As D is antisymmetric, this can be 
viewed as a Pfaffian, Pi D. Note that we are saying that 
V det D is holomorphic in t [again, except for the factor 
(Imr)i/4 in the case]. 

These partition functions can also be conveniently cal- 
culated from a Hamiltonian point of view of the (non- 
chiral) massless Majorana field theory in imaginary time 
[l3]- Then [Fn^.m^l^ are given by, respectively, 

T,qLo-c/2i-Lo-c/2i ^^_y 

Tr(-l)^g^°-'=/24^^"-^/24 (_+)^ (2.40) 

where Lq ± Lq are proportional to the Hamiltonian (with 
the energy of the ground state in finite size for — spa- 
tial boundary condition subtracted off) and momentum 
of the system, F is the fermion number, c = 1/2 is the 
central charge, and in the first two cases the trace is over 
states in the sector with either the — or -I- boundary con- 
dition in the space direction, respectively. The factor of 
V2 in Foi arises because for -I- spatial boundary condi- 
tion there is a zero energy mode shared between right and 
left-moving fermions that can be occupied at most once, 
giving two states. These functions arose in the same way 
as partition functions of edge Majorana fermion modes 
of the MR state on the cylinder [47| , which emphasizes 
again the relation between bulk and edge of the QH 
states. For the final function |Foo|^ for ++, one must 
consider the trace weighted with (—1)^ as for [-Fioj^, but 
with fermion operators ipip inserted also, so that the zero 
mode does not give a zero answer. This treatment of the 
zero mode produces the factor (Imr)!/*. 

The wavefunction of the A'^-particle component of the 
normalized BCS state is then, in the three even cases, 
TOi, 7712 not both zero, 

Fmi ,m2 ('^)-P^ Qttii ,m2 

(r,-r,|r). (2.41) 

We emphasize that the normalization is correct (up to a 
T-independent constant at least) even on passing to the 
fixed-particle number components, provided the particle 
number is chosen close to (say, within ^/N of) the ex- 
pected value in the BCS ground state (which we assume 



is large), by standard arguments about the equivalence 
of canonical and grand canonical approaches. For these 
even cases, these functions are precisely the three confor- 
mal blocks for TV Majorana fermions on the torus, for all 
even N including = 0, whereas the function without 
the factor Fmi.m.2 is the "normalized" chiral correlator, 
in which the block is here divided by the correspond- 
ing one for = (see the discussion in Sec. II C|) . For 
TOi = m2 = 0, there are similar statements for odd A^, 
but the "normalized" chiral correlator involves dividing 
by the A^ = 1 function Fqq, which seems less natural. 
Note that the function Fqo itself is the conformal block 
for A^ = 1, a single Majorana fermion on the torus, which 
is independent of the coordinate zi . 

Now we introduce the familiar fact that the 
parametrization of torii of given area by t is redundant: 
infinitely many different ts describe the same torus, up 
to an isometry. The transformations from one to an- 
other of these are of the form r {ar + b)/{cT + d), 



where 




is a matrix of integers of determinant 1. 



These matrices form a group, SL(2,Z), called the mod- 
ular group. (The transformations arise by changing to 
a different basis for the lattice generated by uji, 0)2.) 

SL(2, Z) is generated by the elements ^ = ^ q \ ^ '' 

S" = ^ ^ 0^ ) ' correspond to the transformations 

T : r — !■ r + 1 and S : t ^ ~^/t] it can be defined 
abstractly as the group generated by elements S, T with 
relations (ST)^ ^ S'^, = 1 (where 1 is the identity 
element). Note that the action on t of any matrix in 
SL(2,Z) is unchanged by multiplying it by — 12, so that 
strictly speaking the group of modular transformations 
on T is the quotient PSL(2, Z) = SL(2,Z)/Z2, where 
Z2 ~ {±12}, defined by generators S, T with relations 
{ST)^ — — 1 [we will not distinguish between S and 
T in SL(2, Z) and their images in PSL(2, Z)]. 

Given our states |toi , 7712, t) on the torus for each r and 
for each boimdary condition (mi, 7712) (the system size 
will not play a role here, provided ^iL^, ^iLy are large, 
and will be dropped from the notation), we can consider 
the effect on them of modular transformations. This is a 
warm-up exercise for the statistics calculation later. First 
we consider the monodromy of the states viewed as func- 
tions of T. Starting from any state [ttji, 7712, t), we can 
vary t along a path to reach either t-|- 1 or — 1/t, keeping 
(7771, 7^2) fixed. Note that (tmi, 7^2) are defined using the 
parametrization of the torus determined by r. The final 
result is a torus equivalent to the original one, but the 
boundary conditions in terms of the original r are now 
different. The effect is T : (701,7)7.2) — > (7771,7)7.2 — tii), 
and S : (7771,7772) — ^ (jn2, —mi). When reduced modulo 
2, these also describe the eff'ect on the boundary condi- 
tions around the two cycles that are important in the 
long-distance part of the state. The three even cases are 
permuted by these transformations, while the odd (-1 — h) 
case maps to itself (it was clear this must be so, because 
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fermion number is conserved by these operations) . 

If we consider the winding numbers only modulo 2, 
then the three even states form some representation 
of the modular group, while the odd state is a one- 
dimensional representation. First we consider the nor- 
malizing factors Fmi,m2- The transformations of these 
under monodromy, simply comparing -Fym .moij + 1) or 
Fmumii-'^l'r) with i^„i;,„i^(r), are known [l3]; the ac- 
tion of T is easily found to be 

^^ii(t + 1) = e~^^''^^Fw{T), 
i^io(r + l) = e-2^'/4«i^ii(T), 
^^oi(r + l) = e'^^I^^Foiir), (2.42) 

while the action of S is 

Fi^{-1/t) - i^^nM, 
^^io(-1/t) = ^oi(t), 
Fmi-l/r) = Fw{t). (2.43) 

Thus the coefficients in these transformations form uni- 
tary matrices, which generate a unitary representation 
of the modular group. The sum X^mi.ms l-^™i,m2p (over 
TOi^2 = 0, 1, not both zero) is a modular invariant, and 
is the partition function of the critical Ising model. 

The function i^oo is V2(Inir)^/'* times the Dedekind 
eta function, 

oo 

ri{r)=q''^^X{{l-q^), (2.44) 

n=l 

which is a modular form of weight 1/2 (with multi- 
plier system — see e.g. Ref. (3) for the modular group: 
,y(T -h 1) = e2W24^(^)^ r/(-l/T) = V^r]{T). As 
Im(— 1/t) = Imr/lrp, {lmTY^^ri(T) transforms purely 
by phases under the modular group, like the other i^mi,m2 
above, though in this case the phase depends on r. The 
significance of these facts will be explained below. 

Next we turn to the transformation of the pairing func- 
tions g under modular transformations implemented sim- 
ply by monodromy. The transformation of the various 
functions can be easily obtained [11], and for z = is 
similar to that of 77(t) (i.e. the weight is 1/2), but also 
involves the appropriate permutation of boundary condi- 
tions. For general z, the t?s are said to be automorphic 
forms under the larger group of translations of z and 
modular transformations that act on both z and r [4^ : 
for example, one has 

^1/2,1/2{z\t+1) ^ e^""/"^?l/2,l/2(^|T), 

A/2A/2(z/t\ - 1/r) - -(-zr)i/2e-^V- 

X ^i/2a/2{z\t). (2.45) 

In the ratio of i!)s giving the pairing function, most such 
factors cancel, the exception being due to the derivative 
of t?i/2,i/2 in the numerator. Consequently, one finds for 
all TOi, 7712 that the pairing functions transform under the 



generators S, T, simply by the permutation of boundary 
conditions, already described above, and by a factor of 
t/\t\ in the case of S. The factor of t/\t\ arises from the 
interplay of the elliptic function and the factor (Imr)"'^/^, 
similarly to the case of ri{T) discussed above. Let us 
point out that the net number of factors (Imr)^/^ in the 
A^-particle component of the wavefunction is just N, for 
all cases, including the boundary condition in which 
N is odd. These factors are essential in order that the 
monodromy matrices be unitary, as they clearly must be 
when orthonormalized states are analytically continued. 

The interpretation for the phases in the transforma- 
tion under S is then that, apart from the N- and t- 
independent factors arising from the -F'mi,m2S, the phase 
(t/|t|)^/^ is connected with the rotation of the geome- 
try that brings the torus back to an equivalent one un- 
der S; the rotation angle is — argr, and the factor for 
each particle reflects the conformal weight 1/2 for each 
Majorana fermion in the conformal block. Put another 
way, each Cooper pair carries angular momentum — 1, 
and thus each particle effectively carries (orbital) angu- 
lar momentum —1/2. This is true even for the unpaired 
particle in the case of N odd, though there the factor 
arises in the BCS point of view through the normaliza- 
tion factor. No such phase arises from T as no rotation 
of the system is involved. 

Thus we have determined the monodromy representa- 
tion of the modular group on our ground states. Read- 
ers are cautioned that because the full states depend on 
(7771,7712) as integers, and not just modulo 2, we can- 
not say that the representation on the states is three- 
dimensional in the case of paired superfluids. In the QH 
case, this would be true, but there one must also con- 
sider the charge sector or Laughlin factors in the wave- 
functions, which we will do in the following section. 



C. Hall viscosity and adiabatic modular 
transformations 

Finally we are ready to consider the more physical op- 
eration of adiabatically dragging the states to perform 
a modular transformation of r. Because of the intrinsic 
interest of the result, we will present two arguments that 
both give the same result. In the first, we will take the 
long-distance forms derived above, and assume that they 
apply at all length scales (down to a short distance cut- 
off, so that the state can be normalized). The functions 
are then conformal blocks, and an elegant calculation is 
possible. In the second argument, we instead consider a 
general p — ip paired state, and show that the result can 
be obtained quite directly in the language of the pairing 
function. Thus the second argument is very general, and 
applies even for the strong-pairing phase, and for other 
paired states, and even in higher dimensional space. 

Now we begin the first argument. We have already 
constructed for each r a set of four ground states that 
are orthonormal, and depend holomorphically on r, ex- 
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cept for a factor (ImT)^/"*. The transition functions are 
determined by the S and T transformations given above. 
The adiabatic transport can be considered separately for 
the one odd and three even cases, as these two types 
clearly cannot mix with one another; here we will denote 
the four as ^'a(T) for brevity. 

We can now apply the reasoning of Sec. II Bl varying 
r instead of quasihole coordinates. If the wavefunctions 
were completely holomorphic in r, then the Berry con- 
nection would vanish. As they are not holomorphic, this 
is not the case, but as in the quasiparticle transport cal- 
culation, the non-holomorphic dependence takes a simple 
form, and one can easily take it into account, while still 
exploiting the holomorphy of the remaining more com- 
plicated factors in the wavefunction. (We should also 
point out here that the integration measure should be 
independent of t. If we write z — Lx{pi + TP2), where 
pi and p2 are in the unit interval [0, 1], then the integra- 
tion measure for each becomes (Pvi = Adpidp2, and 
the variation of r is always performed with A fixed. The 
crucial point is that in this parametrization, the wave- 
functions are still holomorphic in t up to the same factor 
as before.) Then for the Berry connection we have 

= ldA^a\^b)^ '^fj^'^^ {^a\^b), (2.46) 

4Imr 

where we used the holomorphy of 5* up to a known factor. 
Because the states are orthonormal, this reduces to 



Ar,ab A-i-f-i^i) 



Slmr' 



(2.47) 



(In the paired states with indefinite particle number, N 
can clearly be replaced by its quantum average N .) The 
result happens to be the same as in Ref. [6|, in which 
the filled LLL state has the same form, (Imr)^/'* times 
a function holomorphic in t (however, we will see that it 
is something of a coincidence that the power of Im r is 
exactly the same) . Following the reasoning of that refer- 
ence, the curvature (or field strength) of this connection, 
which is — 7V/(2ImT)^, corresponds to a Hall viscosity of 



hn 



(2.48) 



Here we restored h to obtain physical units for viscosity, 
and n is the average particle density. We emphasize that 
the viscosity came out independent of the aspect ratio 
T and system size L^, as is appropriate for an intrinsic 
local property of a fluid. 

We may now also consider modular transformations, 
applied adiabatically. Because of the curvature of the 
Berry connection that we just calculated, the adiabatic 
effect (holonomy) depends on the path taken in the 
upper-half plane). To remove this effect, we can consider 
special paths that enclose vanishing area. That is, the 
"moduli space" of distinct torii is the upper half plane, 
Imr > 0, modulo the action of PSL(2, Z), which is gen- 
erated by S and T. As S sends r to — l/r, we see that 



T = i is a fixed point of S. We may consider adiabatic 
transport along a path lying in the fundamental domain 
llmrj < 1/2, \t\ > 1, that connects two points on the 
unit circle that are mapped to one another by S. If the 
path is a semicircle centered at i with radius shrinking 
to zero, then it implements the S transformation but the 
line integral of the Berry connection along the path van- 
ishes (we note that the Berry connection is smooth at 
i). Consequently, for this limiting path, the holonomy 
is once again given solely by the monodromy of the con- 
formal blocks. We may make a similar argument for T, 
using a straight path from —1/2-1- ilmr to 1/2-1- ilmr, 
and let Imr — > 00. Then once again, the line integral 
of the Berry connection gives zero, and in the limit the 
action of T by holonomy is the same as the monodromy. 

These results then fully characterize arbitrary modular 
transformations implemented adiabatically. To describe 
the results succinctly, we can check the defining relations 
of the modular group in our representation. This requires 
that we compose S and T operations. To do so, it is 
convenient to view all group elements as implemented by 
paths in the upper half plane that start and end at the 
same point, which we will take to be t = e^^^^. To reach 
the points i and 1/2 + ioo at which we defined S and 
T above, we follow a path along the unit circle (resp., 
Rer = 1/2), apply S (resp., T), and retrace the path 
afterwards. Then the Berry phases related to the Hall 
viscosity cancel. [In similar calculations for QH states, a 
gauge transformation is sometimes required as part of the 
adiabatic transport (holonomy), because of a net rotation 
of space, to return to the original basis set. But for the 
operations (ST)^ and S'^ we wish to check, the rotation 
is by TT, under which commonly- used gauge choices are 
invariant, and this simplifies the calculation.] We find 
then that as holonomy equals monodromy, we can read 
off the transformations, and verify that (ST)^ = S"^ = 1 
in the present case. Thus we do have a representation of 
the modular group. 

In the argument so far, we treated the long-distance 
wavefunction as if it were valid for all particles, and so 
wavefunctions are just conformal blocks. One may won- 
der if the result is different when the form of the pairing 
functions is different at short distances, or just how gen- 
eral (or universal) the result is. For this reason we now 
include a second argument that shows that the result for 
the Hall viscosity is, in fact, completely general. The adi- 
abatic modular transformations are then also invariant, 
within a given topological phase. 

We will calculate the Berry connection in the upper- 
half plane directly from the normalized BCS ground 
states. As the different boundary-condition states are 
orthonormal and the final result is diagonal, we will 
simply drop the index a for simplicity. Using the ex- 
pressions for the normalized BCS ground state jfi) — 
[det(l + g'^ g)]~^/^\D,) (note that this may differ from the 
T-dependent family of states used above by a r-dependent 
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phase), we obtain 



A, 



-Tr [5^9,5 - (a,5t)g](l+gt^)-i. (2.49) 



We evaluate this in k space, in which case g is a diagonal 
function of k, and the values of k are determined by 
kx + iky — 27ri(n2 — tiit) / {LxhuT). The values of ni^2 
are ni = rrixjl (modulo integers), and similarly for ^2 
(for 00, i.e. H — h boundary conditions, ni, 71,2 both zero is 
omitted from the sum). Now for a, p ~ ip state, the gap 
function Ak is of the form kx ~ iky times a function of 
|kp, while the kinetic energy is a function of |kp only. 
Hence the pairing function is of the form [k^ + iky)^^ 
times a function of |kp. Further, the only dependence 
on T is through the discrete values of k. In A^., we see 
that the dr acting on the function of |kp cancels between 
g and contributions, and the only remaining term is 
from the prefactor {k^ + iky)~^ in g. The result then 
reduces to 



A, 



ni 



9^9 



4 ^-^ no — niT 1 + 

ni,n2 



9'9 



(2.50) 



As Lx 00 with r fixed, the leading part of this is 



Ar 



A 



8Imr 

TV 



dPk 



9'^9 



(2^)2 1 + g^g 



ST ' (2.51) 
81m r 

which agrees with the previous result. We should men- 
tion that the result here for -f-|- boundary condition has 
A'^ — 1 in place of A'^, as the unpaired fermion is not 
counted by Tr g'^ g{l + g'^ g)~'^ = N —1 for ++, and this is 
a discrepancy with the calculation in terms of conformal 
blocks. This is because we were somewhat sloppy in the 
last step, replacing the sum by the integral, and a more 
careful treatment of the region near k = would lead to 
similar calculations as the preceding argument. But in 
any case, in the thermodynamic limit this discrepancy is 
negligible. 

We can see from this calculation that the Hall viscosity 
is unchanged throughout the weak-pairing phase, and is 
the same in the strong-pairing phase. For different angu- 
lar momentum values in the pairing, and even for mix- 
tures of different values, it generalizes easily to give Ti/2 
times the average angular momentum per pair, times half 
the particle density. A similar form also holds for pairing 
in higher space dimensions, in which the direction of the 
orbital angular momentum of the Cooper pairs is a vec- 
tor. Thus the Hall viscosity is a consequence of the net 
orbital angular momentum of each pair about its center 
of mass. The effect appears to be overlooked in standard 
references 42|] on He^. 



D. Higher genus surfaces and mapping class group 

The discussion for the torus generalizes to oriented 
surfaces (without boundary) of higher genus Q (genus 



1 being the torus). The set-up was discussed in Ref. [1]. 
The fermions should experience a net magnetic field of 
Q — 1 flux, in order that the ground state contain no vor- 
tices. That is. A, which can be defined within coordinate 
patches, is non-vanishing everywhere. The gap function 
A is then a part of the Dirac operator on the surface. It 
depends on a choice of a boundary condition, which can 
be either -f or — for each of the 2Q cycles on the sur- 
face; we will label these by a, and write Aq. Of these 2^^ 
boundary conditions or "spin structures", 2^~^(2^ -I- 1) 
are "even", and the remaining 2^"-'^ (2^ — 1) are "odd". 
In the odd cases, the Dirac operator A^, has one zero 
mode (a holomorphic function, really a section of the 
bundle), and in the even case it has none. Consequently, 
there is a paired ground state with even particle number 
in the even cases, and with odd particle number in the 
odd cases The pairing function is proportional to 
1/Aq, as for the plane, sphere, and torus. The general 
theory for h = —fi pairing still applies, and so the impor- 
tant long-distance part of the normalizing factor in the 
ground state wavefunctions can be taken to be \/ det Aq 
(with the zero mode deleted in the odd cases). This again 
is the (chiral) partition of a massless Majorana fermion 
field, and the normalized N-partic\e wavefunctions are 
chiral partition functions with N insertions of the Majo- 
rana field ip, which are essentially the conformal blocks 
of the Majorana CFT. 

The last statement contains the word "essentially" be- 
cause one further point was neglected. We saw for the 
case of the torus (t/ = 1) that the chiral partition function 

(without insertions), \/ det Aq , tends to a limit which is 
a conformal block, apart from a factor of the form of the 
exponential of the area times a constant independent of 
the parameter r; this factor is all that remains of the 
cutoff. For 5 7^ 1, the corresponding statement does not 
quite hold. There is an additional residual dependence 
on the scale size A^/^ of the system (relative to the cutoff 
scale). It is given for a general conformal field theory by 
a factor A'^^^^'^, where c is the central charge (c = 1/2 
for the Majorana fermion CFT), and x is the Euler char- 
acteristic of the surface. This originates from the trace 
anomaly, which relates the expectation value of the trace 
of the stress tensor, T^j, of the CFT, which generates a 
change of scale, to c times the local curvature of the man- 
ifold evaluated [l7[ . When integrated, the latter gives the 
Euler characteristic x = 2 — 2fJ of the surface (see e.g. 
Ref. [i^l, which includes also the more general case of 
a surface with a boundary). Notice that this result ap- 
plies to the sphere, with Q = 0. More generally, there 
is a dependence on the local cutoff scale give n by the 
exponential of c times the Liouville action [50l . [sii ] ; this 
is independent of the boundary condition choice labeled 
by a. Such factors are not expected to contribute to the 
holonomy, when we vary the geometry holding the area 
A constant, so we ignore them. 

A genus G > I surface endowed with a complex struc- 
ture is described by 3G — 3 complex parameters called 
moduli, in place of the single r for Q = \ [5l|. The 
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moduli space A^g of inequivalent surfaces of genus G is 
described by a domain denoted Tg (TeichmuUer space) in 
Q3e-3 jjiodulo the action of a group Tg of equivalences 
analogous to the modular group, called the mapping class 
gro up (or TeichmuUer modular group): Mg = Tg/Tg 
|52 . [53] ■ The structure of Tg is difficult to describe for 
Q > 1, especially for Q > 2. The conformal blocks are 
holomorphic in the moduli, up to a factor analogous to 
(Imr)^/^ which is raised to the power N (when consid- 
ered at fixed A). While the explicit functions are diffi- 
cult to obtain, we expect that the curvature of the Berry 
connection is again given by the Hall viscosity, which is 
independent of the genus of the surface, and that the rep- 
resentation (which is actually projective, and unitary) of 
the mapping class group defined by the holonomy of the 
adiabatic transport is the same as the monodromy rep- 
resentation, up to path-dependent phases determined by 
the Berry curvature. 



E. Generalizations and strong-pairing phase 

In this subsection we briefly consider other paired 
states of fermions from a similar point of view, temporar- 
ily leaving behind trial wavefunctions that are conformal 
blocks. The paired states we consider are states that 
(except for /i = 0) are fully-gapped, and result from a 
gap function that is a rotational eigenstate in momen- 
tum space, like the p + ip state above. These include 
pairing of any odd angular momentum i for spinless or 
spin-polarized fermions. In the weak-pairing phase, the 
ground state on the torus can be normalized by a fac- 
tor \/det"A^, where again a labels the possible bound- 
ary conditions. By the same argument as above, this 
equals (PfDy. This is clearly equivalent to £ Majorana 
fermions. For i odd, as required here, these phases are 
non-Abelian, as pointed out in Ref. j^. However, since all 
the Majorana fermions see the same boundary condition, 
the theories are not simply direct products of £ copies of 
the Majorana fermion theory above; rather, they possess 
a hidden SO(^) symmetry. Similarly, fermions with spin, 
or other multi-component fermions have other possibil- 
ities, such as the 331 p-wave state related to the Dirac 
fermion, which is Abelian. Another case is spin-singlet 
pairing for spin-1/2 fermions, where the angular momen- 
tum £ of the pairing must be even. In the weak-pairing 
phase, the long-wavelength normalization factor for the 
ground states on the torus is (detD)^. This may be 
viewed as 2£ Majorana fermions, and there is a hidden 
S0{2£) symmetry [thus S0(4) for £ = 2,the d + id-wave 
case; the argument given here is more direct than that 
given in Ref. [3]]. These cases are all Abelian. The hid- 
den symmetries under SO(M ) for some M are present in 
the topological properties of the bulk of the paired states 
(e.g. in counting numbers of conformal blocks), but may 
be fully or partially broken by the dynamics at the edge. 

By contrast, in the corresponding strong-pairing 
phases, the pairing function g is non-singular at k — > 



(and short ranged in coordinate space), and the normal- 
ization factor (|2.36p has no singular behavior at k = 0. 
Moreover, all four ground states occur at even particle 
number. Thus the ground state wavefunctions are con- 
centrated on configurations in which all fermions have 
positions equal in pairs (within a coherence length), and 
the long-wavelength part of the ground state wavefunc- 
tions is trivial. These correspond to conformal blocks in 
a CFT that is trivial (except for the charge sector). The 
orthonormal ground states are not holomorphic in r, but 
have T-dependence related to the Hall viscosity already 
discussed above. They also have non-trivial monodromy 
(and hence also holonomy) in so far as they are sim- 
ply permuted by modular transformations, without the 
phase factors that originated from the non-trivial nor- 
malization factors Fa (under T). This is not the same 
as possessing no holonomy whatsoever. We also point 
out here that the normalization factor is non-singular at 
k = at the transition point fj, — 0, even though g has 
weakly singular behavior. (The leading singular behavior 
in (7 at k ^ cancels in g^g; possibly there is some sub- 
leading effect.) Thus this point appears to have the same 
holonomy as the corresponding strong-pairing phase. Of 
course, adiabatic transport is presumably irrelevant here, 
because the excitation spectrum of the fermions is gapless 
at the transition point. 

The general framework for describing topological 
phases is that of modular tensor categories (MTCs) [H, 
0, HE] (strictly, for QH systems, this is for bosons; some 
modifications are required for QH systems of fermions, 
which are not fully modular whenever the chiral alge- 
bra of the CFT contains the fermion, becoming a chiral 
superalgebra |l|). For paired s-wave superfluids (when 
viewed as a fully-gapped topological phase, such as when 
a Coulomb interaction is present) , Kitaev has argued [s^ 
that the correct MTC is the "toric code" or Z2 gauge 
theory, with Abelian fusion rules corresponding to the 
group Z2 X Z2, and we find that this applies here for 
strong-pairing phases (without the charge sector) also; 
it accounts for the four ground states. Other than the 
identity, or trivial quasiparticle, this theory contains two 
quasiparticle types that are bosons, call them "electric" 
and "magnetic" charges (these charges do not refer to 
particle number), and one that is a fermion, under ex- 
change of a type with itself. Adiabatically making a 
circuit of a magnetic around an electric charge gives a 
phase of —1, and hence the composite of an electric with 
a magnetic charge makes the fermion, which picks up 
a —1 factor under a circuit around either an electric or 
a magnetic charge. (This behavior was calculated from 
trial wavefunctions in a lattice model that resembles the 
strong-pairing phase [st}.) 

When the charge (particle number) sector is included 
to make a QH system, from say the p+ip state, the filling 
factor is V = l/Q, where Q is a positive integer, and is 
even when the particles are fermions, and odd for bosons. 
It was argued in Ref. [4] that the strong-pairing phase is 
equivalent to a Laughlin phase of charge- 2 bosons at fill- 
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ing factor 1/(4(5). This relation predicts the statistics of 
the family of fractionally-charged excitations generated 
by the quasihole of charge 1/(2(5) under fusion; those of 
particle number m/(2(5) have fractional statistics phase 
r7i^7r/(4(5) under an exchange. The fusion rules for this 
theory are Z4Q. However, at first sight it does not appear 
to account for the excitation that consists of an unpaired 
composite (neutral) fermion, which exist in this phase as 
well as in the weak-pairing phase. Under adiabatic trans- 
port in a circuit around the quasihole of charge 1/(2(5) 
(corresponding to the vortex in the paired state), the 
neutral fermion should pick up a phase — 1 (in the weak- 
pairing p + ip phase, this follows from the calculation in 
Section |lVl we will assume this here in general, and also 
in the strong-pairing phase) , and this phase clearly does 
not come from the charge sector. 

These statements can, however, be reconciled. The 
neutral fermion times the quasihole of charge 1 gives us 
back the underlying particle (a boson if Q is odd, and a 
fermion if Q is even). Further, using the assertions just 
made, this particle is local with respect to all the quasi- 
holes (that is, a phase -t-1 is produced by a circuit of 
one around the other). Therefore, quasiparticles that are 
related by fusion with the particle should be identified. 
This leaves a total of 4(5 distinct types, with fusion rules, 
charges and statistics as in the charge sector part alone 
(however, the formal charge which was defined modulo 2 
equals the actual particle number charge only modulo 1). 
Usually we would have said that the particle is in the chi- 
ral algebra, which leads to such identifications. However, 
in the CFT of the edge excitations, the bulk particle does 
not appear (and is not part of the chiral algebra), and 
the theory is simply that for the charge sector as in Ref. 
[3, [2^ . We believe this is the correct description of the 
strong-pairing phases that involve pairing of composite 
fermions. (For spin-singlet phases, the fermion carries 
spin-1/2, but this does not affect the MTC properties.) 



III. HALL VISCOSITY OF THE LAUGHLIN 
STATES 

In this section we calculate the Hall viscosity of the 
Laughlin states. Using the ideas already presented above, 
the main point is to evaluate the normalization of the 
Laughlin states on the torus, up to constants indepen- 
dent of the size and aspect ratio r. We will find that 
the result is again given by a spin density, where in this 
case the spin of each particle is (5/2, the conformal weight 
of the field representing the particles in the CFT point of 
view [l[ (these states include the filled LLL, with Q = 1, 
as a special case providing a check on the calcula- 
tion). We discuss the general situation for QH states, 
and adiabatic modular transformations. 

The key to obtaining the normalization of the Laughlin 
states in this sense is the plasma mapping from Laugh- 
lin's original paper Q , which here will be generalized to 
the torus. Thus we will rely on the physical argument 



that screening holds in the plasma when Q is not too 
large. Most of the calculation is aimed at establishing 
that the wavefunctions on the torus have modulus- 
squared that is the Boltzmann weight in a plasma with 
a neutralizing background, with control over L^;- and r- 
dependent factors that are commonly omitted. Then the 
normalization integral equals e~''^^°^'^\ independent of 
T. A particular issue to be dealt with is the neutraliz- 
ing background charge (indeed, if the plasma were neu- 
tral when only the mobile charges were considered, the 
results would be almost immediate, as we will describe 
later). We will first obtain the Boltzmann weight for N 
point charges Q and a~^N = (M integer) smaller 
charges of —aQ on the torus. The latter will be fixed 
on the points of a grid. As a ^ 0, these become the 
uniform background. By taking this limit carefully, we 
can determine the necessary normalizing factors. Even 
though the potential due to the neutralizing background 
appears to be a constant because of translational sym- 
metry, its dependence on t has to be evaluated, with the 
area A fixed at the correct value for N particles. 

The Boltzmann weight for a Coulomb plasma in two di- 
mensions can be obtained from a 2D conformal field the- 
ory, the (non-chiral) free massless scalar field (p [13, 13| ■ 
Its partition function is given by the functional integral 



Z 



where the action is 
S 



(3.1) 



(3.2) 



The normalization is such that the two-point function in 
the plane is 



(^(r)^(O)) =-ln|r|- 



(3.3) 



Some of the important operators in this theory are the 
electric charges, given by Oe{Y) = e"^''''^'"^/^ for some 
charge e (i? is a parameter that will be used later). These 
have conformal weights h = h = e^/(2i?^). 

We now consider this theory on the torus of side Lx 
and aspect ratio r as before. The field (p will be com- 
pactified, which means that the configurations ip to be 
integrated over must be allowed to obey generalized pe- 
riodic boundary conditions (for convenience, we write z 
for r, even though the fields are not holomorphic) 

p{z + Lx) = p[z) + 2TTRn, (3.4) 
p{z + Lxt) = p{z) + 2TrRn'. (3.5) 

Here i? > is a parameter, the compactification radius. 
The meaning of this is that, if V0 is an "electric" field in 
two dimensions, these boundary conditions allow for elec- 
tric fields that circulate the system, without any source. 
This is important in obtaining periodic correlation func- 
tions when charges are present. A general field configu- 
ration for these boundary conditions can be written as 



ip : 



fnn' 



(3.6) 
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where 'p'^ni is a fixed solution to the boundary conditions, 
say 



„, = 27ri?Ini 



nr 



while ^ is a periodic function. For any n, n' , we have 



(3.7) 



G(r) = (^(r)^(0)> = -ln 
The function £{z/ Lx\t), 



£{z/L,\t) 



£{z\r) 



9.^?i(0|r) 



-7r(Im z) /Im r 



(3.8) 



(3.9) 



where we have used the more familiar notation 'di{z\T) = 
—'&i/2,i/2{A''')i is periodic up to phase factors, and is 
called the prime form for the torus (or else this name 
is used for the ratio of ds without the Gaussian). The 
function G(r) is the Green's function for the Laplacian 
on the torus, defined as the inverse of — in the space 
orthogonal to its zero mode. 



Sir) - J 



(3.10) 



At short distances, it is asymptotically the same as in 
the plane. 

Now we evaluate the (unnormalizcd) correlation func- 
tion of N charges e = Q at positions Zj, and N/a 
charges —aQ at positions Wk, where it is convenient to 
choose N/a — M^, with M > an integer. We take 



(The terminology for the charges thus 



agrees with that used earlier, with as the tempera- 
ture.) We will assume that Q is even, and we will recover 
the Laughlin states for bosons. The states for fermions 
require a slightly different sum over boundary conditions, 
and one eventually obtains a unitary representation of 
a subgroup of the modular group [47|, whereas in the 
bosonic case one obtains a representation of the full mod- 
ular group (hence the correlation functions we are calcu- 
lating are invariant under these respective groups). The 
unnormalized correlation function can be evaluated by 
standard techniques El], including a Poisson sum- 
mation which replaces the n' sum by a sum over integers 
771. Omitting some algebra, we obtain 

i k e 

(3.11) 

where e runs over the Q integer values —Q/2 + 1, —Q/2 + 
2, . . . , g/2, and 



^e/QfiiQ^/L^lQ'^) ^-^Q (Im Zf/A 

7y(r) 



k<l 



where Z = YjiZi -oY^f. Wk , z^j = Zi- Zj , and Wij = Wi- 
Wj (compare Ref. |59l|). The evaluation of the sums over 
77, 77' has led to a sum of modulus- squares of functions 
^'e that we will see are, when a — > 0, essentially the Q 
Laughlin ground states on the torus. The sums over 77, 
77' render the correlation function doubly periodic in the 
Zi variables, though not in the Wk variables, because of 
the fractional values of the corresponding charges which 
enter through Z. Hence we should be careful to specify 
precisely the positions Wk- We will assume that Wk = 
0. 

Now we relate the functions ^'e to the standard forms 
for the Laughlin states on the torus [H, and take 
the a ^ limit. First we deal with the Gaussian factors. 
As they stand, all of them actually cancel. This was 
expected for a conformal block of a collection of charges 
with zero total charge. But we wish to separate out the 
i^feS, and take a limit. So let us assume that J^k — 0, 
and omit the ik and kl Gaussian factors. The others then 
simplify to 



(3.13) 



j^-NQ/2-aNQ/2 



(3.12) 



which is the usual Gaussian in the LLL in the Landau 
gauge (we used N/A — v/{2tt; note we implicitly selected 
a gauge when extracting from I^Pep)). In addition, 
the factors diizij / Lx\t)^ and the center of mass fac- 
tor -d^QiQZ/ Lx\Qt) also occur in the Laughlin states 
[2^ . [Sq . The periodic Laughlin states can be defined for 
arbitrary phases in the boundary condition in each direc- 
tion [sll , and we have obtained particular values of these 
phases, which depend on J2k '^k, so that other values are 
also available if desired. Thus the Zi dependence of these 
functions is exactly that of the Laughlin states. 

We can also see that the powers of become 
A-^Q/4-"^Q/4(Imr)^Q/4+"^Q/4. As a -> (with N 
fixed), this yields a power of A times (Imr)^'^/'*^ 
which the exponent is related to the conformal weight 
(5/2 of each Oq that represents the particles in the same 
way as in the paired states. The correlation function 
is modular invariant, and this entails that the vector of 
functions ^'e transforms by a unitary matrix. As we saw 
in connection with the paired states, the factors Imr 
do not transform by phases under S, and there must 
be modular forms of non-zero net weight to compensate 
them, in order that the transformations involve t/\t\ (to 
the power of conformal weights). Noting that the prime 
forms contain 92i9i(0|r) = 2nr]{T)^ [i^l, we can extract 
77(t)^ from each £, that is a net factor 7y(T-)^Q+"^Q^ 
and this combines with the power of Im r to give the 
expected transformation behavior. (Another way to see 
this is to realize that each factor \£{z/ L.j;\t)/ L^l'^ must 
be modular invariant when the area is held fixed, so that 
£{z / Lx\t) / Lx must be invariant up to a phase, and the 
exponent of t/\t\ is related to conformal weight.) The 
remaining factors, including the e-dependent one in the 
center of mass, are all ratios of the form 7?/?7(r), which 
transform with weight zero. Now we want to show that, 
of these remaining factors, the product of those that de- 
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pend on Wk tends to 1 as a —s- (up to r- and L^- 
independent factors). 

As we already stated, the uniform neutralizing back- 
ground will be simulated by placing the charges on a grid, 
for which we choose the lattice sites Lx{ni +n2r)/Af , for 
integers ni and ri2. The ratio 'di/rj of interest can be 
written as 



r]{T) 

7r(Im , 
oo 



(3.14) 



— 7r(lmz) /Imr 27riT/12/ iTTZ —inz\ 



27ri7iT-\~2'iTiz 



)(l-e 



27ri7iT~27Tiz\ 



n=l 



using the product formula for [M 
that 



Then we can show 



M-l M-1 

n n^" 



z + 



M 



m— 712—0 

-i7rRoT(M-l)(2M-l)/6-i7rM(M-l)Ro z-i7r(M-l)(3M-2)/2 



X £{Mz\t). 



(3.15) 



Similar results, but differing by a phase, hold if any of 
the set of points on the grid in the complex plane are 
shifted by integers or by integer multiples of r; this is ac- 
tually needed if we use the above choice of a set of Wk such 
that Wk = 0. However, these phases are not impor- 
tant. We now apply the above formula, with z replaced 
by Zi/Lx for each i, in the above expression for ^'e. In 
that case, it must be raised to the power aQ — NQ/AP. 
Ignoring the overall phase, we see that f (Mzi/ij^jr)^/*^ 
vanishes when Zi is at any of the points on the grid, but 
(roughly speaking) tends to 1 as M — > oo away from 
these points. This is a little crude, as the points of the 
grid themselves depend on M. If we rescale Zi so that the 
grid points in the Zi plane are independent of M, then the 
statement is clear. Moreover, for large M, the fraction of 
the plane (or of a unit cell) of the grid at which the func- 
tion is smaller than 1 — e [for any e in (0, 1)] is of order 



(1- 



for large M, which is exponentially small in M. 



Without the rescaling, it follows that SiMzi/L^lr)'^/'^ 
tends to 1 at almost every Zi, that is except on a set of 
measure zero; the latter set is essentially the set of points 
Lx{o.i + (12t) where ai and 02 are rational numbers. 

The factors containing can be handled similarly, 
with both Wk and wi on the grid. There is the slight 
complication that k = I must be excluded. This can be 
handled by taking the derivative of the identity (|3.15p at 
z = 0, and extracting a factor 9zi?i(0|r) from both sides. 

Putting all these results together, and dropping phase 
factors and r-independent factors (including A), we find 
that the wavefunctions are 



'i!e{Zi,Wk\T) = 



n 

i<j 



[(Ir 



Nl/4 



(3.16) 



which essentially agrees with Refs. [26|, |58| up to z- 
independent factors, and we emphasize again that differ- 
ent boundary condition phases on the particles can be in- 
corporated by changing the value of — limjv/^oo '^Tlk'^k 
in Z. This was derived for bosons {Q even), but the 
result is similar for fermions (Q odd). 

Let us now summarize, and complete the calculation of 
the Hall viscosity. The form (|3.12p for the wavefunctions 
(before the M ^ 00 limit) shows that for particles at 
separations much less than L^, the interactions in the 
plasma are independent of the system geometry, as they 
should be. Hence the sum of the norm-squares reduces to 
g-^/o(Q)^ independent of the geometry. The other form 
p.l6p exhibits the effect of a modular transformation 
(in particular S) most clearly. Also, it shows that the 
dependence on r is holomorphic, apart from the power 
of Imr and the Gaussian. Using Zi = L^ipu + Tp2i), we 
see that 



-i(Imz,)^ 



(3.17) 



which differs only by a phase from e'^'^^''''^'''^' . Hence 
by using an appropriate gauge choice, the wavefunc- 
tions are holomorphic in r except for the explicit factor 
(Imr)^*^/"*. We have actually calculated only the sum of 
the modulus-squares of the wavefunctions, while we want 
to show that the set is orthonormal. Because the set 
transforms by unitary matrices under modular transfor- 
mations (by monodromy), and these form an irreducible 
representation, we can borrow an argument from Sec. 
IVl below, and conclude that they are in fact orthonor- 
mal (up to a constant). Alternatively, we can obtain 
orthonormality from the known fact that the functions 
are mutually orthogonal, which follows as they have dis- 
tinct quantum numbers under certain translations under 
which they are eigenfunctions ^43,]. Then it follows from 
use of center-of-mass translations, which map them to 
each other, that they have the same normalization con- 
stant, and so are orthonormal up to a common constant. 

By a similar calculation as the first one in Sec. Ill C[ 
it now follows that the Berry connection on the space 
of Q states as a function of r, and its curvature, are 
proportional to the matrix Ig. The value of the Hall 
viscosity of the Laughlin states is finally 



V 



{A) _ 



hQn 



(3.18) 



where the density here is n = l/(27rQ£^) if we restore £b 
which is 1 elsewhere. The result agrees with the special 
case Q = 1 (the filled LLL) which was treated in Ref. Q 
by using a Slater determinant of single-particle wavefunc- 
tions; there should be a version of the Vandermonde de- 
terminant identity for the torus, which will show that the 
wavefunctions themselves are the same, up to an L^- and 
T- independent constant. For Q 7^ 1, the result disagrees 
with a more general result that was claimed recently in 
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Ref. [7[. Further, adiabatic modular transformations can 
be performed using the ideas of Sec. Ill C[ and are given 
by the monodromy, up to the path-dependent phases. 

Finally, we address the more general QH trial states 
that include a conformal block from a CFT. On the torus, 
the charge part is similar to that in the case of the Laugh- 
lin states above, the differences being in the center of 
mass factor, which plays only a minor role here anyway, 
and in t he exponent Q which more generally is replaced 
by (while the denominator Q of the filling factor still 
appears in the center of mass factor). Under modular 
transformations, the conformal blocks behave as auto- 
morphic forms, and at fixed area a factor Imr must be 
included to an appropriate power, so that the S transfor- 
mation acts by unitary matrices that include factors t /\t\ 
to the total conformal weight of the fields. This is rather 
more clear for these blocks than for the Laughlin states, 
where we had to deal with the neutralizing background. 
For general blocks, the factor Lx must appear for each 
particle, by scaling considerations. This produces the de- 
sired power of Imr when written in terms of A, which 
is fixed. (The paired states in Sec. Ill Cl are an example, 
in which — 1/2 for Major ana fermions.) Hence it 
becomes clear that the full functions always include the 
factor 

(ImT)^(''"/2+''-^'/2 (3^^9) 

times a function holomorphic in r. More generally, the 
trial functions with quasiholes included contain similar 
factors depending on their total conformal weight (in- 
cluding the charge sector) . Under monodromy, S is given 
by T-dependent phases because of these factors, and r- 
independent coefficients from the center of mass factors, 
which are also A^-independent and occur even for 0, 
in which case the wavefunctions are the conformal blocks 
for the CFT (including charge sector) on the torus with 
no particle insertions. We will argue in Sec.|V]that these 
trial wavefunctions are orthonormal at large N , at least 
in some cases. The arguments presented here then show 
that the Hall viscosity is generally given by 

This will be a universal result within a topological phase, 
when the trial functions represent such a phase. We em- 
phasize that the result here, and that for the paired 
states, is always a certain constant times the average 
"conformal spin density" . For the trial QH states given 
by conformal blocks, the conformal spin of each parti- 
cle is just its total conformal weight, /2 -\- h^. This 
explains why the p + ip paired state, and the filled LLL, 
both gave the same Hall viscosity (when written in terms 
of density): both involve the same conformal weight, 1/2 
(the value for Dirac or Majorana fermions; for the filled 
LLL, the associated CFT is the Dirac fermion theory, 
which appears e.g. as the theory on the edge). 

We recall that the conformal spin also enters in the 
shift S. More generally, in non-chiral CFTs, the confor- 



mal spin of a field is the holomorphic minus the anti- 
holomorphic conformal weight, which equals the confor- 
mal weight in a chiral theory. For v filled Landau lev- 
els, the result [H, can also be viewed as the aver- 
age spin per particle, where the spin arises from the cy- 
clotron motion [in the JVth Landau level, Af ^ 0, 1, . . . , 
it is {Af + 1/2) and this appears also to agree with 
the classical non-zero temperature result (3^ . For other 
topological QH phases of matter, such as the hierarchy 
or composite fermion phases, the conformal spin density 
in the ground state can be readily calculated by utilizing 
the preceding results, and the average spin per particle 
always equals half the shift. For the gapless Fermi- liquid- 
like phases (gO], for which, in the simplest cases {ly = 1/Q 
with Q even), trial wavefunctions are related to Laughlin 
wavefunctions, we expect that the Hall viscosity takes the 
same form as in the Laughlin states as calculated here; 
however, as these phases are gapless, other components 
of viscosity may be non-zero also. Note that "spin" here 
always refers to orbital effects, as we ignore the intrinsic 
spin of the electrons. We expect similar effects involving 
the "real" or intrinsic spin of spinful particles. 

It is worth emphasizing that the value h times density 
is the natural quantum of viscosity, in any dimension, 
and so is analogous to e^//i for electrical conductance, or 
TrkgT/h for thermal conductance. Values of components 
of the viscosity tensor in a quantum fiuid (or respectively, 
of a conductivity tensor in two dimensions) at some tem- 
perature might be either larger or smaller than this value, 
but the non-dissipative part is quantized at zero temper- 
ature (in the sense of taking a universal value) through- 
out a topological phase, and helps characterize the phase. 
However, viscosity relates to momentum transport, and 
so the Hall viscosity probably does not mean much if 
translational symmetry is violated at short length scales, 
say by impurities, unlike electrical conductivity, which 
relates to charge transport, and charge conservation is 
not violated by impurities (similarly for thermal conduc- 
tivity) . 



IV. ADIABATIC STATISTICS CALCULATION 

In this section we present direct arguments for the adi- 
abatic statistics of two and four quasiholes in the MR 
state on the sphere or plane in the thermodynamic limit. 
The calculations work by "doubling" (taking two copies) 
of the system, which turns the problem (the Ising, or 
Majorana fermion, CFT) into an Abelian one, the Dirac 
fermion CFT, which can be bosonized, meaning turned 
into a neutral Coulomb plasma without a uniform back- 
ground charge density. (The Abelian problem can be 
viewed in the QH context as the 331 state which applies 
for example to bilayer systems, but with the charge part 
removed.) It is straightforward in the Abelian problem 
to argue that various blocks are orthonormal (up to a w- 
independent constant) for well-separated quasiholes, as 
well as holomorphic in w, and that will conclude the ar- 
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gument. Readers who prefer a more conceptual approach 
that is also more general (and avoids doubling) may pre- 
fer to skip this section and go to the following one. 

We will consider here only the trial wavefunctions with 
the charge sector removed. These presumably correspond 
to the long-distance behavior of the p + ip paired state in 
the presence of vortices, though this has not been shown 
directly from the point of view of BCS theory. Arguments 
that the charge part can be dropped from the QH trial 
wavefunctions for the present purposes are considered in 
more detail in Sec. |V] below. The trial wavefunctions 
will be written in terms of conformal blocks from confor- 
mal field theories. We will first recall the relation of two 
copies of the Majorana fermion field theory to a Dirac 
theory, beginning with zero quasiholes, then progressing 
to two, and finally to four quasiholes. 

The relation of the MR Pfaffian wavefunction (with 
charge sector removed) to the Ising (or Majorana) CFT 
was reviewed in the introduction. Here we will also need 
the similar relation for the Dirac CFT. The Dirac fermion 
CFT is defined using the fields -0, -0^ with the opes 

0t(z)0(O) ~ - + ..., (4.1) 

z 

0(z)0(O) ~ + ..., (4.2) 

(and similarly for tp^ip^), in which as usual omitted terms 
tned to zero as z ^ 0. These expressions then also give 
the chiral correlators for the fermion fields, using Wick's 
theorem (and note that the fields anticommute) . The 
theory can be related to two Majorana fermion field the- 
ories ■01 J V'a by writing 

V'^(z) = [Mz) + tM^)]/V2 (4.3) 
0(z) = [0i(z)-z^2(2)]/\/2. (4.4) 

The relation also works in the presence of spin fields a± ; 
by definition, the Dirac fields acquire a change of sign on 
continuation around the location of a spin field. More 
precisely, these are defined by the opes 

0t(z)a_(O) ~ -i^a+(0) + ..., (4.5) 
0(z)(7_(O) - 0+..., (4.6) 

and similarly with 0^ ^ 0, c+ ^ ct-- The spin field 
tT+ (fJ-) carries charge 1/2 (—1/2), in units where 
has charge 1, and (t± both have conformal weight 1/8. 
(The U(l) charge is also referred to as pseudospin in the 
bilayer context.) In terms of the two Majorana theories, 
the Dirac spin fields become the product tTiO'2 of the Ma- 
jorana (Ising) spin fields, which were defined earlier. 

The basic property we wish to exploit is the rela- 
tion of the Dirac fermion theory to a Coulomb plasma. 
In field theory language, this is bosonization, in which 
by introducing a chiral scalar field f'{z), we represent 
0t(^) ^ e*'^'(^), 0(z) = e-'"^'''^\ We note that this 
is the chiral version of what we had for Q = 1 (the 
filled LLL) in Sec. IIIII (but now with ip' as this is not 



the charge sector). Then the spin fields are represented 
by (T±(z) = e^*''' '^^^/^ (all these fields are present in the 
modular- invariant Q = A bosonic theory in Sec. IIIip . 
These relations imply that conformal blocks in the Dirac 
theory can be written as Coulomb gas functions, that is 
as products. This is exemplified by the Cauchy determi- 
nant identity. The conformal block for Dirac fermions on 
the plane with N 0's, N ip^^s, at Zi, i = 1, ... , 2N, with 
i even for 0, i odd for 0^, is 

(0^^ (Zl)0(z2) ■ • ■ 0(^2A'))Dirac = dct ^ 

Z-i Zj 

= / -,NjV(Ar-l)/2 ni<i'(^» ~ ^»')nj<j'(^3 ~ ^j'l , 

^ ' TkM^) ^ ' 

in which in the matrix with entries l/{zi — Zj), and in 
the products, all i's run over odd values, all j's over even 
values; we are not concerned about the overall sign. The 
last equality is the determinant identity, and exhibits the 
block as a Coulomb gas Boltzmann weight, similar to 
those discussed earlier. In terms of the Majorana fermion 
theories, the determinant becomes a sum of products of 
Pfaffians, 

det = 2-^ Pf^^i ^^Pfs. (4.8) 

Zj^ Zj ^ Zi Zj Zi Zj 

Here the sum is over the distinct subsets Si of 
{1,2,..., N} with an even number of elements \ Si\ = iVi, 
and 5*2 is the complement, 5*1 U S'2 = {1,2..., N}] 6*1,2 
refer to the type 1 and 2 Majorana fermions. Pfg^ 2 de- 
notes the Pfaffian with the indices i, j , ranging over the 
set 5*1,2. Evaluating a Pfaffian requires an ordering on 
the index set to fix the overall sign; 5*1 2 are viewed as 
inheriting the lexicographic ordering from {1,2,..., N}). 

The sum over different numbers of Majorana insertions 
01 and 02 in the last expression is clearly inconvenient. 
To remove the technical difficulty, we will treat all the 
states grand-canonically, which is natural anyway for the 
paired states, as in Sec. [Hi To do this, let us introduce 
fermion creation operators cj , cj , for particles with pseu- 
dospin in two space dimensions, with canonical antico- 
muutation relations. These are related to two copies of 
the spinless fermion creation operators (as in Sec. |TT]) by 

c\ = [4(r) +i4(r)]/\/2 (4.9) 
c\ = [4(r)-z4(r)]/V2, (4.10) 

similar to those for the 2D fermion operators 0ji,2 (and 
again, local in r) . (If we think of the Dirac theory as part 
of the 331 QH states, then C| and represent fermions 
in the two layers with definite U(l) "pseudospin" values, 
while ci, C2 are similar to the even and odd combina- 
tions of the pseudospins that are also sometimes used 
[2^, except for the factor of i which will be convenient 
here.) Now using the vacuum |0) annihilated by the C|,| 
operators, an unnormalized BCS state with pairing only 
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between t and | has the form 



n 



) = el '^'"■d=r'M(r,r')c|(r)c|(r')|Q^^ 



(4.11) 



and the 2A^-particle wavcfunction is 



^'(ri, . . . ,r2Ar) = (0|c|(ri) • • •C|(r2Ar)|f7||). 

= detjvM(r.i,rj), (4.12) 

where we have used the same convention that i is odd and 
J is even as in the wavefunctions above. In this case, the 
norm square of any state (not only a paired state) 
with equal numbers of "f and | particles is computed from 
the 2iV-particle functions ^' as 



(%i|f^Ti) 



N 



|^(ri,...,r2Ar)p 



(4.13) 



(which also generalizes to ^ Ni , in which one divides 
by -/V|!7Vj^!). For the paired states above, this norm square 
reduces to det(l + MM^). 

For these fermions, M(r, r') need not have any par- 
ticular symmetry in general. If we choose M(z, z') = 
fi/{z — z') (switching to z's corresponding to r's), where 
/i is a fugacity similar to the chemical potential in Sec.lTTl 
then we obtain the above wavefunctions. If we now use 
the relation to fermions of types 1,2, then the two types 
decouple in irij^), which equals the (tensor) product of 
two BCS states for spinless fermions, as in eq. (j2.10p . 
with g = M: 



(4.14) 



Clearly, this is consistent with the normalizing factors 
for these states, also. In carrying out integrations over 
space in these states, we will disregard boundary condi- 
tions and the integration domain, but one can think of 
them as performed on the sphere, with suitable changes 
in the details of the functions to include rotational sym- 
metry and so on. The vortices will always be separated 
by much less than the system size, and we consider the 
thermodynamic limit at fixed density. In addition, short 
distance divergences as fermions come together must be 
cut off; this will not be shown explicitly. Because the 
leading term in the ope is independent of the positions 
of other fields (including the vortices) , the cutoff effects 
should have no effect on the results. 

With these few remarks, we have completed the expla- 
nation of the relation of the Dirac and doubled Majorana 
theories and conformal blocks (wavefunctions) without 
spin fields (vortices). We should point out that the re- 
moval of the charge sector was technically necessary to 
obtain the tensor product form of paired states. If the 
charge sector were included in the wavcfunction for | 
and I particles (as in the 331 states), it would couple the 
1 and 2 particle types together, and the system would 
not factor into decoupled MR states (however, one might 
argue that this does not matter, along similar lines to 
arguments for dropping the charge sector). 



Now we turn to the introduction of spin fields. The 
simplest functions to write down are those in the 
Coulomb gas form. For t particles (with labels 
i odd), Ni I particles (with labels j even), with 
+1/2 charges (or cr+'s) at positions Wk, k odd, n_ 
— 1/2 charges (or (t_'s) at positions wi, I even, and with 
2N^+n+ = 2Ni+n- for charge neutrality, the conformal 
block (times /x to the power half the number of fermions) 



.(Nr+N, )/2 Ilk<k' K - Wk' Y/^ YIkI' ( 



Wl - Wl 



,)l/4 



Wl) 



1/2 



Wk 



,1/2 



(4.15) 



The 2-dependent factors can also be obtained from a 
generalization of the Cauchy determinant identity. For 



Ni , the same block is equal to (up to a sign) 



nk<k'i^k-Wk')'/'Ui<Am-w,)^/^ 



UkM-wiV/' 

where M{zi, Zj) is given by 



■detAf(zi,Zj), 
(4.16) 



M{zi, Zj) = /i 



(4.17) 



This is easily proved as the w-dependent factors in M can 
be taken outside the determinant. Hence, apart from a 
function of w's only, the grand-canonical state takes the 
BCS form \n^i{M)) for this M. This special case of 
iV| = Ni will be sufficient for our purposes, but can be 
generalized if required. 

Now we can address the normalization (or overlap) 
integrals for the grand-canonical states with A''-particle 
components (I4.15P . For the self-overlap, we can note that 
the sum (over N) of the integrals is simply the partition 
function of a Coulomb plasma, with some "impurities" 
inserted at the Wfc's (note the factors l/{N^\Ni\) are es- 
sential for this). For the values of the exponent used here, 
the plasma is in a screening phase. Then if the impurities 
at Wl are far apart compared with the screening length, 
the overlap integral takes the form 



(4.18) 



similar to that for the case with a neutralizing back- 
ground. Here /o and are constants, and A is the area 
of the system. Note that this is independent of the posi- 
tions and of the charges of the spin fields. 

Though the Dirac theory is Abelian, we can nonethe- 
less obtain not just one but a set of functions from it, by 



choosing a set of positions wi, 



. , and varying 



the assignment of which are +1/2 and which are —1/2 
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charges. This gives us ("^J^"") distinct states. We can 
form the overlap matrix for these states, and so far we dis- 
cussed the diagonal elements. For the off-diagonal over- 
laps, let us suppose we take the overlap of two states 
in which the charge assignments are the same except for 
one +1/2 and one —1/2 which are switched in one state 
compared with the other. Then the phase of the inte- 
grand winds by ±27r as one fcrmion is carried around 
just one of these two locations (by monodromy), instead 
of the integrand being purely real in all previous cases. 
We want to argue that when these locations are far apart 
compared with the screening length, the overlap goes to 
zero, because of the oscillations of phase. 

One way to see this is to consider the sine- Gordon for- 
mulation of this Coulomb plasma, which is represented 
by a field theory with action 



J on 



(4.19) 



in which the screening is due to a term — 2 fj, cos tp' added 
to the action of the (non-chiral) scalar field ip' . In this 
phase one can consider the latter term as simply produc- 
ing a mass term (by expansion of the cosine around a 
maximum). In this language, the phase winding corre- 
sponds to the two ws in question being locations of "mag- 
netic" charges that produce a "vortex" in tf', instead of 
"electric" charges. That is, (p' changes by ±27r on go- 
ing around a magnetic charge. This should be viewed as 
allowed in the functional integral over ip' , by identifying 
values (^'-f 27r with ip' . But as the cosine term tries to pin 
ip>' to a multiple of 2n, all of this winding has to be acco- 
modated by a "domain wall" on a straight line connecting 
the two vortices (with exponentially decaying tail off the 
wall) . The system then incurs a positive domain- wall free 
energy proportional to the length of the wall (when these 
vortices are separated by more than the screening length, 
which is also the thickness of the wall) , which means that 
the overlap we were computing decays exponentially with 
separation. This argument clearly generalizes to any of 
the off-diagonal overlaps. Hence the matrix of overlaps 
Zab (in the language of Sec. II Cp is proportional to the 
identity, and as the wavefunctions are holomorphic in 
ui's we can obtain the holonomy of these states under 
adiabatic transport of the w's by reading it off from the 
monodromy of the wavefunctions, that is from the fac- 
tors {wk — Wk')^^^^, where the exponent depends on the 
signs of the charges. 

The main problem now is to relate the conformal 
blocks of the Dirac and doubled Majorana theories with 
spin field insertions. Once we do so, we can use the or- 
thonormality of the blocks in the Dirac CFT to infer the 
same for the blocks in the Majorana CFT, which will 
complete the argument. 

First we will consider the simple case of only two vor- 
tices (spin fields), at positions wi and W2- Above, we took 
-1-1/2 at wi, —1/2 at ^2- If we compare with the simi- 
lar state with these charge assignments reversed, which 
means in the wavefunctions that we exchange wi and W2, 



and take the sum of the two states (with a certain phase 
relation to be determined in a moment), then the nor- 
malization changes only by a factor of V2 when the spin 
fields are far apart, using the above argument. Now in 
the form of the functions as paired states, we can write 
for M{z, z') 



M{z,z') = M{z,z')+ L{z,z'), 



(4.20) 



where M is antisymmetric, M(z',z) = — M(z,z'), and 
L is symmetric, L{z',z) — L{z,z'). Moreover, for the 
case here of n+ = n_ = 1, M is symmetric under the 
exchange wi ^ W2, while L is antisymmetric. Explicitly, 



M{z,z') 
and 



I \l/2/ ' \l/2 / ^l/2/ / n1/2 

{Z — Wl) ' [z —W2) I {z — W2) {z —WI) ' 

1 {z-W2)'-^^{z'- Wi)'-/^ (2-«Ii)l/2(z'_^2)l/2 



Z — Z 



Wi - W2 



2Ui=iAi^- my/Hz' -wiyf^ 



(4.21) 



(4.22) 



Now when we sum the states over the exchange of wi,2, 
the phase is chosen so that if the factor {wi — W2)~^^^ 
is omitted, the functions are simply added. Then in the 
expansion in many-particle wavefunctions, all terms with 
an odd number of L's drop out. L{z, z') is non-singular 
as 2 ^ z'. Thus, in the language of Ref. 26\, in the terms 
containing L's some ti pairs are broken, and the unpaired 
"f fermions occupy a mode with constant wavefunction (as 
do the i fermions) . But as only even numbers of L factors 
can appear in our state, the mode must be occupied by 
more than one t fermion, and so these terms vanish on 
antisymmetrization (i.e. by Fermi statistics). 

The wavefunction of the state thus has thejorm of the 
above paired state, with the antisymmetric M in place of 
M. Hence, it is equal to the tensor product of two paired 
states of spinless fermions containing the same g = M. 
The form of each of these paired states is therefore 



{w,-w2r'^'MM)) 



1,2- 



(4.23) 



The form of M is simply fi times the propagator for the 
Majorana fermion in the presence of two spin fields. The 
factor {wi — W2)~^^^ has exponent equal to minus twice 
the conformal weight 1/16 of the spin field, and the N- 
particle wavefunction is exactly the conformal block for 
N Majorana fermions and two spin fields (here N is even) 
1]. Because the doubled version of this state is normal- 
ized for well-separated spin fields, and holomorphic, we 
can conclude the same for this state. Hence we can read 
off the phase for adiabatic exchange of the two vortices, 
g-i7r/8^ This is a simple phase factor, yet does not come 
from the charge sector of a QH state, which was removed 
here. This contribution was expected and emphasized 
early on [l| , and previously derived by a different method 
[3l| . We note that the same result can be obtained for N 
odd instead of even. This is obtained by taking the linear 
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combination of the two charge assignments for wi , W2 in 
the Dirac theory that is orthogonal to that used above, 
in which case ah the terms with exactly one L survive, 
and give rise to blocks in the Majorana theory that all 
have odd particle number N (these are non-vanishing in 
the presence of a positive number of spin fields). These 
correspond to occupying the zero mode with one fcrmion 

For the case of four vortices in the paired state, we use 
similar but more involved arguments. We begin by in- 
troducing notation for the six states with n+ = ~ 2. 
We will write the BCS states previously called |ri||(Af)) 
as |Oi3,24), and the remaining five are obtained by per- 
muting the indices; thus the first two subscripts are the 
positions of the +1/2 charges, and the last two are the po- 
sitions of the —1/2 charges. The grand-canonical states 
are in full 



I* 



n 



13,24/ 



k<k' 



1/4 T-r 
"kk' llKV 



„l/4 



n 



1/4 
k,l ''^kl 



13,24/ 



(4.24) 



in which again fc, fc' — 1, 3, and Z, 1' — 2, 4, and five other 
states obtained by permuting the indices; these form an 
orthonormal set when the spin fields are well separated. 
Once again, we can make use of the decomposition of 
M{z, z') into antisymmetric and symmetric pieces under 
z ^ z', called M(z, z') and L{z, z') respectively. These 
are respectively even and odd under exchange of all of 
the coordinates of +1/2 with coordinates of —1/2 charge 
spin fields, that is wi <-> W2j ''^'3 W4. Explicitly, 



M{zt, Zj)i3,: 



(4.25) 



Ul){z — W3){z' — W2){z'—W4) 
V2){z-W4,)(z' -Wi)(z' -W3) 



1/2 



+ (z ^ z') 



-(4.26) 



L is non-singular as z — > z' for symmetry reasons. When 
we write the states using the 1, 2 types of fermions in- 
stead of t, i, each M term involves types 11 or 22, while 
each L term involves types 12. Hence if we sum and form 
the combinations like 



\n 



13.24/ 



24,13/ 



(4.27) 



(and two similar ones), then these are sums of states 
that contain only even numbers Ni and N2 of the 1, 2 
fermions. They are symmetric under exchanging 1 with 
2 type fermions. There are just two Majorana conformal 
blocks for even fermion number (and also two for odd 
fermion number). We are able to show that these three 
states in the doubled theory are constructed from the 
symmetrized tensor products of the two states in each of 
the spinless theories that produce the Majorana blocks 
(similar to the construction of the triplet of spin 1 states 
from the symmetrized tensor product of two spins 1/2 in 
angular momentum theory). This is sufficient to deduce 
the orthonormality of the two Majorana blocks for N 
even. 



Explicitly, the two conformal blocks for four spin fields 
and N fermions (N even) in the Majorana theory were 
found by NW [2g|, who used a related doubling tech- 
nique. The expressions are 



N/2 



•^0,1/2 — 



1 



W13W24 



1/8 



\/l ± VT~^ VW12W32W34W14 

Pf M(zi,Zj)i3,24 

iVT^Pf M(z„z,) 14,23] , (4.28) 

where x — u'i2W34/(wi3W24) is a cross-ratio, and the up- 
per sign is for 0, the lower for 1/2. For iV = 0, in which 
case the Pfaffians are replaced by 1, these two blocks 
were found in Ref. f6l| (see also Ref. [46]). From these, 
the four-point spin correlation function of the critical 
Ising model is X]s=o 1/2 l-^^l^- same combination 

should give the correlation function in the presence of 
fermion (energy operator) insertions. Then if we are to 
obtain holonomy equal to monodromy of these blocks, 
these states should be an orthonormal pair [2^ . We will 
denote the two grand-canonical states corresponding to 
these TV-particle wavefunctions by |±). 

The relations of the Majorana and Dirac blocks, when 
there are no fermions, also dictate what the relations 
must be when fermions are present, most conveniently for 
the grand-canonical forms. After a considerable amount 
of algebra, we can verify that 



W13W24 



1/4 



\w12w32w3iWii , 
l+>i 



1) + |f^24,13)) = 



~Wl4W23 



(u;i2W42W43Wl3, 

I+>1 



1/4 



(1^^14,23) 

)|+>2-| 



i» 



(4.29) 



(4.30) 



-W12W34 



1/4 



W13W23W24W14 , 

I+>1 



(|f^l2,34) + |f^34,12)) = 



(4.31) 



The left-hand sides of these are three orthonormal states 
(in the usual sense), and it is then easy to see that |+)i 
and |— )i are an orthonormal pair (in the usual sense), 
and of course similarly for type 2 fermions. As we have 
explained, this implies that the holonomy for an adiabatic 
exchange of well-separated vortices in the paired states of 
spinless fermions equals the (non-Abelian) monodromy. 

Finally, we can consider more than four quasiholes, 
at the level of counting arguments. For n ~ 6, similar 
to n = 4, the number of symmetrized tensor products 
of Majorana conformal blocks with Ni^2 even equals the 
number of Dirac blocks with n-|_ ~ n_ = 3, when they 
are symmetrized between + and — spin fields (i.e. both 
numbers equal 10). In this case a similar argument may 



29 



go through. But for n > 8, the use of n+ = ri- = n/2 and 
iV| = N-^ and symmetrization produces a set of orthonor- 
mal states the number of which is less than the number 
of Majorana blocks (with both even) symmetrized 
between 1 and 2. We are then forced to consider all par- 
ticle numbers together. In the Majorana theory, there 
are 2"/2-i blocks for N even, and to these we can add 
the number for N odd. The square of this is 2", the total 
number of blocks in the Dirac theory if we allow all pos- 
sible assignments of the charges on the spin fields, with 
the numbers and N-^ chosen to ensure neutrality. The 
functions must agree because of the interpretation of the 
Dirac theory as the doubled Majorana theory, and hence 
we expect that a version of the argument goes through 
for all values of n. 



V. GENERAL CONDITIONS FOR HOLONOMY 
TO EQUAL MONODROMY 

In this section we discuss general conditions for the 
holonomy of trial wavefunctions given by conformal 
blocks, including those for QH systems that include the 
charge sector, to equal their monodromy. In view of the 
preliminary discussion in Sec.[Tl it is sufficient to find gen- 
eral conditions for the conformal blocks to be orthonor- 
mal in the basis in which the monodromy matrices are 
unitary. After giving the general arguments, we turn to 
examine a number of examples, and discuss non-unitary 
and irrational CFTs. 



A. Relevant and irrelevant perturbations of the 
CFT 

We begin with the functions '^a from Sec. II C[ 
which include the charge sector, or with the functions 
J-a{wi, . . . ; zi, . . .) which are the blocks of the CFT with 
the charge part removed. In either case, we can view 
these as related to a perturbation of an underlying CFT 
by some operators (viewed within two-dimensional field 
theory); here by the underlying CFT we mean including 
the charge sector contribution, which as we mentioned is 
a massless scalar field [l| . Then the particles are repre- 
sented by tpc{z) — e^^"~^ '''(^)-0(z), where v = P/Q is the 
filling factor (if we view the CFT with the charge sector 
removed as a theory of Abelian anyons, then these parti- 
cles are represented by V' alone) . Similarly, the quasiholes 
are represented by some fields (j){z) = e^9<ii'V(2)r(z) (or 
just by T, respectively; it should be clear how to general- 
ize to include several types of qhasihole). (These "fields" 
are really "chiral vertex operators" [13] ■) Here (/? has two- 
point function {{p{z)ip{Q)) = — Inz, is determined by 
single-valuedness of m the particle coordinates, and 
the charge (particle number relative to the ground state) 
of the quasihole is —qqhy/v- The functions '^a are repro- 
duced as a chiral correlator of these fields together with 



the background charge density factor 

exp[-i j cPz' ^^{z')/{2t:)]- (5.1) 

the approach to handle this factor is discussed in MR [l[ 
(for the functions J- a, the background charge density is 
omitted). 

The overlap integrals Zab represent the inner products 
of these a, viewed as quantum mechanical states of 
N particles (there are similar ones without the charge 
sector). It is convenient to introduce a grand-canonical 
point of view [6^ . as in the previous sections. Then the 
overlap integral can be written formally as an expression 

x^{w^)^{w^)---), (5.2) 

where the expectation is taken in the underlying CFT, 
and the notation for a choice of blocks is suppressed. The 
fields with bars "over them are antiholomorphic partners 
of the chiral fields. Expanding the exponential in ■(/'oV'e 
gives a series of integrals over N particles, iV = 0, 1, .... 
For the plane the background charge is correct as written, 
but for the sphere one should define the series with the 
appropriate background charge, spread uniformly on the 
sphere, for each particle number, so as to ensure total 
charge neutrality (related to the total number of mag- 
netic flux quanta through the sphere in the QH effect). 
For a given (real) value of A, the dominant terms in the 
sum cluster around some value of N , and we choose A 
so that this is large, and the fluctuations in N will be 
of order ^/N . In some cases, the correlator due to the 
non-charge part will be nonzero only for N congruent to 
some number, modulo another constant (for example, N 
even, for the case of no quasiholes in the MR state). 

Leaving aside the background charge density for a mo- 
ment, the overlap integral has a form similar to a cor- 
relator in a perturbed CFT. In fact, the diagonal sum 
X^a is thc (un-uormalizcd, in the sense of Sec. II C|) 
correlator in the diagonal version of the underlying CFT, 
perturbed by adding the term Xtpci^c to the Lagrangian. 
For example, omitting the charge sector, if ip were a Ma- 
jorana fermion, this term would be the mass term that 
drives the continuum theory of the Ising model off its 
critical point. For the charge sector, we have in addition 
the background charge density, which we will view as 
also part of the perturbation. This notion of a perturbed 
CFT will be crucial to our approach from here on. 

In general, given a CFT and a perturbation, the basic 
question is whether the perturbation is relevant, irrele- 
vant, or marginal. This analysis (and even referring to 
it as a perturbation) presupposes that the perturbation 
is weak, i.e. that the coefficient A is small. Let us once 
again ignore the charge sector entirely for a moment, or 
assume that it has been removed. Then the perturbation 
is by the operator '^p^l). It is relevant if its scaling dimen- 
sion 2h^ is less than two, irrelevant if this is greater than 
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two, and marginal if it is two. This is because the scahng 
dimension of A is the codimension d — 2h^ = 2 — 2h^, 
and the perturbation is relevant if the codimension is 
positive, because then the coefficient grows under renor- 
malization group (RG) rescaling transformations to ap- 
proach larger length scales, with the opposite for the ir- 
relevant case. For marginal perturbations, one must go 
to higher (non-linear) order in A to see whether it grows 
or shrinks under the RG on going to larger length scales. 
If the perturbation is relevant, then a crossover length 
scale ~ X-'^/i'^-^'^^) can be defined. We should note 
that the conformal blocks are singular as two ips approach 
one another, diverging as z~'^'^^~^^ as discussed in the 
introduction. This may be strong enough that a cutoff 
must be introduced, by cutting tube-shaped regions out 
of the integration along these diagonals, and also along 
those at Zi = wi, with radius say e. Then the limit 
e ^ requires renormalization of the operators. We are 
suppressing dependence on e from these expressions, by 
holding it fixed. The parameter A may be viewed as de- 
fined at the cutoff length scale, and the crossover length 
is defined in units of e. Then ^x is much larger than 
e when A is small (and relevant). For our purposes, ^x 
can be viewed as the typical spacing of the insertions of 
'0V' in the integral. Note that while the conformal blocks 
are viewed as conformally covariant, conformal, and in 
particular scale, invariance of the correlators is violated 
by the presence of the perturbation. 

It is now clear that in general we will need to exam- 
ine the relevance of the perturbation. When the charge 
sector is included, the situation is different. The per- 
turbation is certainly relevant, but the scaling dimension 
of e^^^iv+v) in the underlying CFT, which is is 
apparently of little importance. From the basic example 
of the Laughlin state [8], we know that the mod-square 
of the Laughlin wavefunction is a one-component plasma 
with a background charge density, which is in a screen- 
ing phase when the exponent q is larger than about 70. 
For larger values of , a crystalline state is formed. 
A scaling analysis does not seem to be of much use in 
finding this physics. (We note that Ref. [1^, section 8, 
began to formulate arguments about the relevance of the 
perturbation, but incorrectly states that ip^ has negative 
scaling dimension.) In contrast, for the two-component 
plasma, which is charge-neutral without a background, 
the physics is the Kosterlitz-Thouless theory fGSj, which 
at small A ("fugacity") is determined by the relevance 
of the perturbation; when it is relevant (which would be 
< 2, if we use the same conventions but include 
charges of the opposite sign instead of the neutralizing 
background), the plasma arrives at a screening phase. 
This was emphasized in Ref. [g^l in a context similar to 
the present one. 

For the more general cases when the CFT includes 
more than the charge sector, and the operator ipe con- 
tains ip ^ 1, it appears that as long as z/ is such that the 
charge sector is in the screening phase when tp is ignored, 
this will still be true in the presence of ^p. The correla- 



tors in the CFT with the charge sector removed do not 
contain high powers of z as one tp{z) is moved away to 
large distance, and so are unable to affect the screening 
properties very much (except perhaps near the transi- 
tion to the crystal). This is generally taken for granted 
in treatments of trial wavefunctions in the fractional QH 
effect. 

It is now natural to adopt the following analysis. To 
consider the effect of the perturbation of the underlying 
CFT, including the charge sector, we begin by noting 
that if z/~^ is less than about 70, the RG flow will be to 
screening behavior of the charge sector. This sets in by 
the time the screening length is reached; this scale is of 
order the particle spacing for flUing factors not too far 
(say within a factor of 10) from I. On large scales, we 
then should consider how the remaining CFT degrees of 
freedom behave. The operators ip are still attached to 
the charges in the plasma. If this was ignored, we would 
integrate over the positions of the tps as if the char_ge 
sector were absent. Then the scaling dimension of ip'4' 
would come into play. Because the insertions of ipijj are 
in fact on the charges in the plasma, fluctuations in the 
density of the ipips, which would be Poissonian if there 
were no plasma, are in fact suppressed if one looks at a 
length scale larger than the screening length. In the pres- 
ence of the charge sector, the screening length is the scale 
from which the RG analysis must be started, rather than 
the scale e which tended to zero. As this scale is simi- 
lar to the spacing of particles, and thus of ipip insertions, 
this corresponds to a fugacity for the ^p'lp insertions of 
order 1. In this situation, it is hardly surprising that the 
fluctuations in the number density of insertions are not 
Poissonian. This is somewhat like a lattice model with 
■ipip on each lattice site. In effect, the strength of the per- 
turbation in the CFT (with the charge sector removed) 
is here of order 1. 

We are then forced to consider the effect of the per- 
turbation by ipTp when its strength is of order one, rather 
than small. When the perturbation is weak, the RG flow 
will carry the system back to original CFT if the per- 
turbation is irrelevant or marginally irrelevant. If the 
perturbation is relevant, the flow will usually carry the 
system to a different RG fixed point (possibly passing 
close to other fixed points along the way) . If the pertur- 
bation is strong, the RG flow might be to a different fixed 
point (or phase of the 2D theory) than when it is weak. 
In particular, a perturbation that is irrelevant when weak 
might still lead to a flow to different phase when strong. 
An example is the KT theory again, according to Koster- 
litz's RG flow diagram for the two-component plasma 
[6^. On the other hand, a perturbation that is relevant 
when weak seems unlikely to lead to a flow back to the 
original CFT when it is strong. RG flows that travel 
in a cycle back to the original flxed point are very un- 
usual, and im pos sible in unitary theories according to 
the c-theorem |65| . Thus a perturbation that is relevant 
for small A should lead to a flow away from the original 
CFT towards another flxed point also when it is strong. 
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possibly to the same one as starting from weak A. If a 
perturbation that is irrelevant when weak does flow to a 
different fixed point when strong, it may be difficult to 
identify what fixed point that is. 

Thus the simplest assumption to make for the QH sys- 
tems would be that the behavior of the system viewed as 
a CFT without the charge sector, perturbed by 'ipip with 
coefficient A of order 1 is the same as for the same pertur- 
bation with A small. In that case, if A is irrelevant, the 
system flows back to the original CFT, and if A is relevant 
it flows to the same fixed point different from the start- 
ing CFT, even though the perturbation is not weak. We 
note that this point of view is essentially a generalization 
of what has been used successfully in studying composite 
particle approaches to the QH states, in which the charge 
sector is factored off, and results for composite particles 
in zero magnetic field (to correspond to the present case) 
are used as a "mean field" approximation [23]. But we 
will learn from simple examples later that this plausible 
assumption can fail in practice, in that even for cases in 
which A is irrelevant (when weak), there can be flows to 
some other fixed point. So we will not proceed by using 
this assumption. On the other hand, if there are sys- 
tems in which we are willing to use the conformal blocks 
without the charge sector as trial wavefunctions (with 
cutoff e much less than the particle spacing), then the 
relevance or irrelevance of A comes directly into play. (In 
real non-QH systems, there could be similar problems 
also, for example for paired states in the weak-pairing 
phases. But for the paired states as studied in Sec. [TTl 
this difficulty does not in fact seem to occur.) 



B. Long-distance fixed point 

Whether or not the simplest assumption for the effect 
of a perturbation holds when its strength is of order one, 
we can analyze the possible results for the holonomy at 
large length scales. If we refer to the starting CFT as 
the short distance or ultraviolet (UV) fixed point, then 
we will now be considering the possible long-distance or 
infrared (IR) fixed point reached at large length scales 
(possibly after the RG flow passes near other fixed points 
at intermediate scale) . There are generally two classes of 
IR fixed points to consider. One is a "massive" fixed 
point, in which all correlations go to constants with ex- 
ponentially decaying corrections. The other is a fixed 
point that contains "massless" degrees of freedom, corre- 
sponding to a critical theory, possibly a conformal field 
theory. In this there are power-law correlations, though 
we should point out that the power could be zero, and in 
all cases there will be power-law subleading corrections 
also. We emphasize that these terms are applied here to 
the behavior of two-dimensional (2D) field theories. A 
further possibility is that of exponentially growing corre- 
lations. This cannot happen in a unitary theory, but may 
be allowed in non-unitary theories; in non-unitary CFTs, 
one finds power-law correlations that increase with dis- 



tance (negative scaling dimensions), so it is not clear that 
exponentially growing behavior is ruled out in the "mas- 
sive" cases. This would probably indicate an instability 
of the vacuum of the 2D theory. We will not have much 
to say about this possibility. 

1. Massive IR fixed point 

We can consider the cases with or without the charge 
sector together. In order to obtain holonomy equal to 
monodromy for the trial functions given by conformal 
blocks (possibly including the charge sector, on both sides 
of this relation), we need the overlap integrals Zab to 
approach Sab at large separations of the ws (up to a w- 
independent factor). We now interpret this in terms of 
the RG flow of the underlying CFT. The behavior is very 
much like correlation functions approaching a constant 
in the perturbed phase, where the correlation functions 
would however be ^^^^a (as noted in Ref. [2^, where 
the role of more than one conformal block is neglected at 
this point). Thus we will suppose that the RG flow goes 
to an IR fixed point that is "massive" , that is one in which 
correlations go to constants or to zero, with exponentially 
small corrections at finite separation. (No power-law cor- 
rections are possible in the absence of massless degrees 
of freedom in the IR fixed point theory.) 

We wish to consider then the behavior of the over- 
lap matrices Zab in such a phase. Important conditions 
are placed on the possible behavior of these by the mon- 
odromy properties they possess. It is easy to see that they 
inherit the monodromy properties that they had in the 
starting CFT (in which in the conformal blocks, the fields 
Tpc or Tp were either omitted entirely, or located at fixed 
positions). That is, under a braiding operation, either 
exchange of identical quasiholes, or a circuit of one type 
of quasihole around another one, the overlap changes by 

2ab Y.^M^'^acZcdMdb, (5.3) 

where Mab is the monodromy matrix, describing the cor- 
responding effect on the conformal block. Now as noted 
above, Zaa is a correlation function of local operators, 
which goes to some constant. As the braiding matrices 
are unitary, it is invariant under monodromy. (If the 
braiding matrices in the monodromy were not unitary, 
it would be impossible for the holonomy in the quantum 
mechanical theory to equal the monodromy!) In a pure 
phase (in the sense here of 2D statistical mechanics), this 
correlation approaches ri/(^(''^')''"('^')); a product of the 
expectation values of the local operator. If the perturbed 
phase is not pure, this factorization may not hold, but 
the value of the correlation approaches the product of the 
expectations that would hold in one of the pure compo- 
nents of the phase, when the phases are related by sym- 
metry (for example, the low-temperature phase of the 
Ising model). More generally, all the components of Zab 
can be considered as (linear combinations of) correlation 
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functions of some operators. The non-trivial monodromy 
(|5.3p implies that the operators are not all mutually lo- 
cal. (Notice that Abelian factors in the braiding drop 
out of this argument). In this situation, it is well-known 
that the correlator cannot go to a constant. Thus, if the 
phase were pure, it should reduce to a product of expec- 
tation values of the various operators. As each expec- 
tation value is independent of the locations of the other 
operators, the product cannot have the required mon- 
odromy. Consequently, the overlaps must go to zero at 
large distances (exponentially), except for combinations 
that are invariant under all allowed braiding operations, 
which can survive as constants. 

The simplest example of this is the Ising model again, 
in which ipijj (for ijj the Majorana fermion) is the pertur- 
bation. The Ising order operator is cr(u;, id), and there is 
also the disorder operator fi{w,w); these two fields are 
dual to each other. Under a circuit of one around the 
other (along a path containing no others of either type) 
in a correlator containing both, the monodromy is a fac- 
tor — 1. The fields are mutually non-local, and cannot 
both have expectation values in the same phase, a has a 
non-zero expectation in the low-temperature phase, while 
fi has one in the high-temperature phase. On the sphere, 
there are 2"/^~^ conformal blocks for the chiral a, so 
2"^^ components Zab in the overlap integral matrix. In 
fact there are this many blocks both for even and for odd 
fermion numbers in the blocks (for n > 0), while n must 
be even. In correlations of as and fis in the critical theory 
(A = 0), there are 2" possible choices (either a or ^) of 
the field at each wi. However, interchanging the types a 
and iJ, leads to the same functions by duality, so there are 
only 2"~^ distinct functions. Further, in the absence of 
any fermions {N = 0), the number of as and of /zs must 
both be even (this is because fusion of the nonchiral a 
and IJ, produces either a ip or a, 'tp [Toj). More generally, 
the number of fis plus N must be even, and (hence) also 
the same for as. Hence there are only 2"~^ distinct non- 
zero functions. This essentially shows that the terms in 
Zab correspond to the various ways to choose order and 
disorder fields in the correlator. In more detail, in the 
grand-canonical version of the perturbed partition func- 
tion, both even and odd TV occur. The odd values of N 
drop out if we sum the cases A, —A. The counting of 
functions still works, but one cannot say if it is the a or 
the /i that has an expectation value. The only non-zero 
limit of large separation corresponds to the correlator of 
all as (or all /xs), and this is equal to the trace J^a^aa- 
If one does not sum over the two signs of A, but instead 
fixes this sign, then it turns out that the two functions 
(iV even or odd, but close in value) are nearly equal, and 
so for one sign of A they cancel, leaving a correlation 
function that decays with distance (exponentially, with 
the correlation scale of order the separation of tpip inser- 
tions) . For the other sign of A they add. This is how the 
distinct behavior of, for example, aa and /x/i correlations 
in either the high- or low-temperature phase arises. Du- 
ality acts by turning a to fj, and /i to a, and reversing the 



sign of A. 

Returning to the general argument, if the overlap ma- 
trix approaches a constant at large distances, then this 
must be invariant under monodromy. Now suppose that 
the unitary monodromy representation of the braid group 
on the conformal blocks is irreducible (as it is in many 
CFTs). It now follows from Schur's lemma that the only 
possible invariant form for the overlap matrix is a con- 
stant times Sab (here we view the matrix Zat as a map of 
the vector space of conformal blocks into itself, defined 
in the given components) . If the representation is not ir- 
reducible, then we can find a (unitary) change of basis of 
the space of conformal blocks so that the representation 
is block diagonal with irreducible diagonal blocks. If the 
blocks are mutually non-isomorphic as representations of 
the braid group, then the overlaps of vectors in distinct 
blocks vanish, and the same argument then applies to 
the overlap matrix of vectors from the same block. It is 
not obvious that the constant multiplying each of these 
is the same. If some of the distinct braid group irre- 
ducibles occur more than once, then the overlap matrix 
can only be reduced to a direct sum of tensor products of 
the identity matrix on each irreducible, times a Hermi- 
tian matrix of overlaps of size given by the multiplicity 
of that irreducible. By a choice of basis that preserves 
the braid group action, these matrices (and hence the 
full overlap matrix) can be made diagonal. The diago- 
nal w-independent constants in this matrix are expected 
generically to be non-zero, so that they can be absorbed 
into rescalings of the conformal blocks, though this may 
be forbidden in conformal field theory, and it is possible 
that there are deeper reasons why the constants must be 
equal. From the point of view of the blocks as trial wave- 
functions, rescaling them so that they are all normalized 
is harmless, and by hypothesis can be done independent 
of w. 

This completes the argument, because when the over- 
lap matrix Zab is proportional to the identity, so the 
basis is orthonormal, as well as holomorphic in w, the 
Berry connection vanishes, showing that holonomy equals 
monodromy. In particular, for the Ising (or MR) ex- 
ample, this shows that the expected non-Abelian statis- 
tics occurs, using only known properties of the massive 
Ising (Majorana) field theory. Generally, we have shown 
that the condition that the IR fixed point of the two- 
dimensional theory is massive (or that correlation func- 
tions of the T insertions go to non-zero constants with 
exponentially decaying corrections) is a sufficient condi- 
tion for holonomy to equal monodromy. It is also a nec- 
essary condition, as a non-vanishing Berry connection is 
unlikely to have no effect on the holonomy. 

The question of whether or not the braid group rep- 
resentation is irreducible is not crucial for the argument, 
as we have seen, but is nonetheless worth a comment. In 
many familiar CFTs, it is irreducible, including in those 
for the entire RR series. More generally, we may make 
a connection with whether a theory, or QH state, is uni- 
versal for topological quantum computation. The latter 
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reduces to the question of whether the braid group rep- 
resentations (for increasing numbers of quasiholes) are 
dense in U(M), where M is the dimension of the space 
of blocks (more accurately, as the U(l) factors cannot be 
dense, it is whether it is dense in the group PU(M) of 
unitary matrices modulo phase factors). If the represen- 
tation is dense, then it is certainly irreducible. Thus in 
all cases of interest for universal quantum computation, 
the representation is irreducible. 

We should add some further remarks on the case in 
which the monodromy matrices M are not unitary. This 
occurs when the correlation functions in the CFT are not 
the sum of modulus-squares of the blocks (i.e. a positive- 
definite sesquilincar form is used to form the norm square 
of the vector of blocks, and the basis choice brings this to 
standard form) as assumed above, but some terms appear 
with a minus sign (this may occur in non-unitary CFTs). 
In this case, the monodromy matrices are pseudo- unitary. 
If the blocks form an irreducible representation of the 
braid group, then the invariance of the large-separation 
limit of the overlap matrix Z under monodromy forces it 
to be a non-positive sesquilincar form also (proportional 
to the same one as in the correlation functions). But the 
overlap matrix is a positive-definite form by construction, 
because of the positive-definite inner product in quantum 
mechanics, and so no such long-distance form is possible. 
This contradiction means that the hypothesis of a flow 
to the massive fixed point is untenable in such a case. 
Of course, as mentioned earlier, it is also impossible for 
the necessarily unitary holonomy to be equal to any non- 
unitary monodromy. 

We also wish to make a comment regarding the map- 
ping class (or modular) group, and the so-called twist op- 
eration on quasiholes. This latter is an essential part of 
the structure of a modular tensor category (or even just a 
ribbon category), or of a topological phase or topological 
quantum field theory (TQFT) in 2 + 1 dimensions. It is 
effectively a counter-clockwise rotation of a quasiparticle 
about its center. For a quasiparticle of a definite type 
a, the twist 9a is expected to be given by e'^'^^^" , where 
ha is the conformal weight of the corresponding field in 
the CFT (including the charge sector). This cannot be 
produced using braiding only, which translates quasipar- 
ticles without rotating them (we generally consider doing 
so in the plane, which is flat). It is nonetheless a physical 
operation on quasiparticles that can be obtained adiabat- 
ically, as we will now explain. 

First, we point out that for any non-negative number 
of quasiparticles at marked points on a surface of genus 
Q > 0, there is a more general mapping class group of 
diffeomorphisms that preserve the complex structure. It 
is generated by the "Dehn twists" , which are obtained by 
cutting the surface along a closed simple curve, making a 
2tt rotation, and rejoining; the closed curves can be non- 
contractible cycles, either wrapping around a handle, or 
enclosing some number of marked points. Thus the twist 
of a single quasiparticle is one of these. In the absence 
of any marked points, this mapping class group reduces 



to that mentioned in Sec. Ill Dl with marked points, it 
includes the braid group for the plane as a subgroup. 
Further, in the presence of marked points, their coordi- 
nates (modulo complex analytic transformations) can be 
viewed as moduli. 

Now the idea is to calculate the effect of any Dehn 
twist adiabatically by deforming the metric of the sur- 
face until one reaches a metric equivalent to the one we 
started with. This involves a generalization of what was 
done when we varied the aspect ratio of the torus. While 
we will not attempt to write down explicit expressions, 
we expect that the conformal blocks depend holomorphi- 
cally on the moduli, except for simple contributions in 
the charge part and some that relate to the conformal 
weight of each insertion, when written with the inflnites- 
imal area element of the metric held fixed (these will 
produce effects involving the Hall viscosity). Moreover, 
the monodromy of the functions will generate a unitary 
"monodromy representation" of the mapping class group 
(strictly, it will be projective, so we should pass to a 
covering group). This can be used to show that the con- 
formal blocks are orthonormal for each deformed met- 
ric, by a similar argument as for the braid group. Then 
adiabatic transport around a closed loop in the moduli 
space, which corresponds to implementing an element of 
the mapping class group, should produce holonomy equal 
to the monodromy, up to a phase related to Hall vis- 
cosity that depends on the path taken in moduli space. 
Applying this to the twist, we expect to show that af- 
ter removing the Hall viscosity part (much like removing 
the phase related to the Magnus force from the statis- 
tics calculation), the twist 9 for each quasiparticle type 
will be the same as in the underlying CFT (including the 
charge sector) . This seems very likely to succeed in view 
of the analogy with the argument for the statistics given 
here, and would essentially complete the derivation of the 
TQFT from the trial wavefunctions given by conformal 
blocks, when the 2D theory flows to a massive IR fixed 
point, up to some remaining details that will be discussed 
in App. [X] Unfortunately, the explicit calculation must 
wait for another occasion. 

There is an easier alternative approach, but it yields 
incomplete information. We can follow the techniques of 
sections IIIII and III C[ and consider modular transforma- 
tions of the ground states on the torus, implemented adi- 
abatically. When the hypothesis of a massive fixed point 
holds, we then find that (as in the CFT) the modular 
transformations obey (ST)^ = = C, where C, which 
obeys = 1, is the "charge conjugation" permutation 
matrix that exchanges type a with its antiparticle a*. 
The eigenvalues of the T matrix are given by 0^6"^'^*'^/^'*, 
where c is the central charge. Unfortunately, we do not 
have an adiabatic calculation of the central charge at 
present. If we know which eigenvalue corresponds to the 
identity object, which has twist 9=1, then the problem 
is solved. In general we can still check the values of the 
ratios of the twists for distinct quasiparticles, and these 
clearly agree with the CFT (including the charge sector) . 
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2. Massless IR fixed point 

Now we turn to the case in which the IR fixed point 
is massless, and thus scale invariant. We note that such 
a theory must in many cases be conformally invariant as 
well, for example, in unitary cases [t^]. Then for a CFT 
we can again apply the machinery of conformal blocks, 
etc. We begin here by considering the case in which the 
perturbation is irrelevant (or marginally irrelevant), so 
that the IR CFT is the same as the one with which we 
began (with the charge part removed). It will be con- 
venient throughout to ignore the charge sector, which 
behaves exactly as in the other cases. 

The irrelevance of the perturbation implies that the 
overlap matrix Zab of interest approaches at large sepa- 
rations of the quasiholes the product of conformal blocks 
for the T fields {without ip insertions): 



Zab ^ J^a{wi, . . .)J^b{wi, . . .). (5.4) 

Thus it certainly does not approach a constant. 

In this situation, we need more general expressions for 
the Berry connection than those we used when we as- 
sumed that the states were orthonormal. Let us write 

Z=JV-^ (5.5) 

as matrices {Af is Hermitian, and we choose signs so that 
it is positive), so that \'^a{w))Afab are orthonormal. 
Then using the holomorphy of I^P <!(?«)), we find the ma- 
trix equations 

= -idAfAf-'^, (5.6) 
A:^ = iAf-^dJ\f, (5.7) 

in which we suppress the index I on A and on d/dwi. 
These generalize expressions for the case of a single block 
in Ref. fio'l, in which there appears to be a sign error. 

On applying this to the limiting form of Z in the 
present case, we see immediately that when there is 
more than one block, there is a problem because, while 
= Z^l"^ exists, its inverse M does not. Let us ig- 
nore this for a moment by considering the two-quasihole 
case, for which there is always only a single block, which 
has the form T\{w\, w^) — {w\ — W2)^^'*^ . Then a short 
calculation shows that the Berry connection is non-zero, 
and in the holonomy the path-ordered exponential ex- 
actly cancels the monodromy of the function. (One can 
normalize the blocks, and then the cancellation is exactly 
the same as occurs if one starts with a trivial example 
and makes a gauge transformation to obtain both the 
connection and the monodromy.) 

Hence the adiabatic exchange of any quasiparticle type 
with itself is trivial. If the trial wavefunctions repre- 
sent a topological phase of matter, then the quasiparticle 
properties are described by a modular tensor category 
(when the underlying particles are bosons; the fermion 
case should be similar). It is possible to show that the 



only (unitary) MTC in which the exchange of two iden- 
tical quasiparticles is trivial for all quasiparticle types is 
the trivial MTC containing only the identity or trivial 
quasiparticle [6^. This shows that if the wavefunctions 
in one of the present cases describe a topological phase, 
then it is trivial. This is the opposite extreme from holon- 
omy equal to monodromy. 

This however is not the whole story. As the IR fixed 
point CFT is massless, one in general expects that there 
will be irrelevant (and possibly also marginally irrelevant) 
operators that decay to zero at large length scales. In 
the situation at present, in which the perturbation of 
the UV theory is irrelevant, these can include the same 
perturbation we began with. But there will be others 
also, especially because our perturbation is not small. 
In general, there is always an irrelevant operator TT, 
where T and T are the components of the stress tensor; 
this operator has dimension 4. Thus there will in fact 
generally be non-zero corrections of the form 



Zab^Ta{w^,...)H{w^....) + , (5.8) 

as r — > oo, where r stands for \wk — wi\ for all pairs fc, 
and if there are other irrelevant operators more relevant 
than TT, then the correction factor will decay as an even 
slower power. (The correction term will not exhibit the 
factorized form of the first term.) When this is included 
in a calculation of the Berry connection, the matrix M 
may be well-defined for any number of conformal blocks, 
and the resulting Berry connection will take the form 
\A\ ^ 0(l/r)+0(l/r'^). In particular, for two quasiholes, 
there will be corrections to the trivial holonomy which 
decay as 1 jr^ . 

More generally, even if the perturbation is not irrele- 
vant, the IR fixed point might be a non-trivial massless 
theory different from the underlying CFT. If so, then the 
overlap matrices will possess power-law corrections (even 
if the leading behavior is going to a constant). If the lead- 
ing behavior is a non-constant power-law behavior, then 
the holonomy will disagree with the monodromy of the 
blocks (which is always that of the underlying CFT with 
which we started). If the massless theory is conformal, 
then there are conformal blocks. Note that when the 
leading part oiM exists, the integral of the leading part 
of the Berry connection contributes just (minus) the log- 
arithm of the monodromy of the blocks of the IR CFT 
to the total holonomy in all cases. It is not clear whether 
there will be any MTC that would contain the result- 
ing behavior (although the subtraction is reminiscent of 
a coset construction jlTj). In all these cases there will be 
power corrections to the leading results for the holonomy. 

When the overlap matrix is viewed as that of trial 
states in a quantum system in 2 -I- 1 dimensions, the 
presence of power law corrections to the holonomy signals 
that the behavior is not that of a topological phase. Such 
a phase is supposed to be fully gapped, and so the holon- 
omy should possess only exponentially-decaying correc- 
tions. The power-law corrections are a clear break-down 
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of the screening properties we were seeking. We conclude 
then that the system is gapless. It may be in a gapless 
phase, or it may be at a phase transition (a subspace 
of codimension at least one) in the space of Hamiltoni- 
ans for 2 + 1 systems (that is, if there is a Hamiltonian 
for which the trial wavefunctions are exact eigenstates). 
In such cases, adiabatic transport presumably does not 
mean much anyway. Possibly, the trial functions simply 
do not represent any respectable phase of matter in 2 + 1 
dimensions, even one at a critical point. 



3. Summary 

Thus we have found that for the two main classes of 
possible IR fixed points, massive and massless, there is 
corresponding behavior of the holonomy: either (i) that 
consistent with a topological phase of matter, and agree- 
ing with the monodromy of the underlying conformal 
field theory, or else (ii) none of these statements hold 
and the 2 + 1 system is gapless, respectively. In general, 
although this is a little crude in case (ii), the holonomy 
in the adiabatic transport can be thought of as the mon- 
odromy of the blocks, which comes from the underlying 
(or UV) CFT, divided by the monodromy of the corre- 
sponding blocks of the IR CFT (the trivial CFT in the 
case of a massive phase) . This is also reminiscent of the 
normalization factors for the paired states in the Sec. 
im which is a ratio referring to the underlying CFT and 
the large distance behavior. Notice that we do not find, 
for example, behavior like that of a strong-pairing phase 
which is topological but trivial (apart from the charge 
part). Trial wavefunctions that are conformal blocks are 
not known for all the possible QH phases, so this may 
indicate that none can exist for the strong-pairing type 
of phases. For the latter, the closest one can get may be 
to use the trivial CFT (i.e. the charge part only) in the 
UV, treating pairs or larger clusters of particles as the 
basic particles represented by V'e- This could be general- 
ized to less trivial CFTs, but still applied to clusters of 
particles. Finally, the third possibility, overlap matrices 
that grow exponentially with distance, may mean a more 
gross instability in the 2 + 1 particle system, so that the 
trial wavefunctions do not represent a true phase of mat- 
ter at all, in contrast to the preceding two possibilities 
which are well defined phases, even when gapless. 



C. Examples 

1. Abehan states 

We begin with trial wavefunctions for some Abelian 
QH phases. We focus on the simplest ones (other than 
the Laughlin states), which are the Halperin mm'n states 
for particles that come in two species or "pseudospins" , 
which we will denote t and |. The states can be ap- 
plied to bilayer QH systems, with the pseudospin being 



the layer index. The trial states exist for non-negative 
integers to, to', n, but we will concentrate on the cases 
to' = TO for which there is a symmetry under the opera- 
tion t^ij i^T for all particles. The trial wavefunctions 
for the ground state (no quasiholes) are 
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The filling factor is 2/(to + n), and the 331 case, which 
means the cases to — n — 2, has been discussed above 
in Sec. IIVI There are two-body local Hamiltonians for 
which the mmn states are exact zero-energy eigenstates 
(see e.g. Ref. [ll]). On removing the charge sector, one 

where f' is a chiral 



has the fields V'T,i = 
scalar field, distinct from, but obeying the same proper- 
ties as, that for the charge sector. The conformal weights 
for these are both (to — rt)/4. Thus the perturbation of 
the CFT is relevant when to — rt < 4, and irrelevant when 
TO — n > 4. We have assumed here that to — n > 0. If 
TO — n < 0, the plasma is unstable 67|. At the same 
time, the scaling dimension of the fields in the CFT be- 
comes negative, and this corresponds to a non-unitary 
CFT. Also, the case m = n represents a ferromagnetic 
state. We assume to — n > from here on. 

When the charge part is removed, this problem is the 
two-component Coulomb gas in two dimensions, which 
was considered by Kosterlitz and Thouless (KT) [63|. 
Then if the perturbation of the CFT was weak, the KT 
transition occurs at to — n = 4. We would then expect 
flow to the original CFT in the IR if to — n > 4, result- 
ing in no screening, and gapless behavior, but flow to a 
screening phase if m — n < 4. However, for larger values 
of the perturbation (in this context, called fugacity of 
the Coulomb charges) the system can flow to the screen- 
ing phase even for to — n > 4, according to the RG flow 
diagram of Kosterlitz, illustrating a point made above. 

On the other hand, if we consider the cases m > 0, 
n ~ 0, including the charge part, then we see that the 
two components of pseudospin decouple, and we have a 
product of Laughlin states. Hence in this case, screening 
does occur, even when to — n = to is large (for to < 70 
anyway). So in these cases it appears that the pertur- 
bation is effectively so large that the perturbed CFT is 
in the screening phase. For n > the plasma is again 
stable [67] , and the widely- held belief that this plasma is 
in a screening phase for both charge and pseudospin is 
consistent, at least. Consequently, these examples form a 
cautionary note that the irrelevance of the perturbation 
does not necessarily mean a flow to a massless phase. 

We note that if the flow had gone back to the CFT, 
it would have been easy (using standard Coulomb gas 
arguments) to see that screening fails in the ip' sector, 
and that the pseudospin associated with quasiholes is not 
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localized near the locations w, but spread over the region 
in between, consistent with a gapless (possibly critical) 
phase. 



2. Read-Rezayi series 

We have already dealt at length with the Laughlin and 
MR states. Now we address the RR series of states [20| . 
labeled by an integer A: > (the other label M corre- 
sponds to additional Laughlin factors, and does not en- 
ter these considerations). For all k, the CFT other than 
the charge part is the parafermion theory. The field 
ip involved in representing the particles in the conformal 
blocks is the first parafermion current, ?/; = "01, of con- 
formal weight — 1 — 1/ k. These underlying CFTs are 
unitary. Moreover, the perturbation by i/'iV'i is relevant 
for all k. However, we note that for fc > 2, the field ip is 
not equal to its own dual (or "antiparticle" ) 4'i (the ad- 
joint taken in the two-dimensional quantum field theory), 
which instead is ip\ — ij^-i = ipk-i- If the perturbation 

were V'lV'i + V'lV'i: then unitarity of the theory would be 
maintained by the perturbation (as it is in the Abelian 
examples just discussed). But for the perturbation we 
have, the perturbed theory is not unitary, except for the 
cases k = 1 (which is trivial), and k — 2 (the MR state). 

For the unitary perturbation i/'iV'i + of the 

parafermion theory, it is generally believed that the sys- 
tem flows to a massive phase, at least when the perturba- 
tion is weak and positive [6g|. If this were the case here, 
we would immediately conclude that the adiabatic statis- 
tics of quasiholes is the same as the monodromy of the 
conformal blocks. But for the non-unitary perturbation 
by ipi'ipi of the unitary underlying CFT, we may have 
some concern (though the coefficient is positive in the 
convention of Ref. [68|)). It is possible for the flow still 
to arrive at the same IR massive fixed point, with non- 
unitarity only manifesting itself in the behavior of the 
exponentially-small corrections, and if such is the case 
then all is well, and we have trial functions for a topo- 
logical phase. In support of this we can only point out 

that the omission of f/'iV'i from the perturbation may not 
make that much difference, as in any case k insertions of 
Tpi can fuse to the identity in an operator product, so that 
the other terms are not needed for "charge neutrality" , 
unlike the example of the mmn states above. In addi- 
tion, we point out that for the charge sector (and thus 
even for the Laughlin states) , the perturbation is not uni- 
tary, because the positive charges are the point particles, 
and the negative ones form a fixed uniform background. 
Nonetheless, for < 70, the flow is to a screening 
phase, which behaves the same way as for the screening 
phase in the two-component plasma which is a unitary 
perturbed CFT, even though the conditions for reaching 
the screening phase are different in the two cases. We 
conclude then that the non-unitarity of the perturbation 
is not necessarily a problem, and the system may well 



flow to the same massive 2D theory as in the case of the 
unitary perturbation, which would lead to the expected 
non- Abelian statistics for all of the RR series. 



3. Blok-Wen series 

In the Blok-Wen series of states [1^ , the particles (of 
spin k/2) are represented in the conformal blocks by the 
operator -0 = 0^/2, where 4>s is the primary field of spin 
s = 0, 1/2, k/2 in SU(2) level k current algebra 
(or WZW theory). (We have suppressed labels for the 
Sz component on these fields.) The perturbation in the 
CFT is the spin-singlet part of 4'k/2't'k/2- Unlike the RR 
series for fc > 2, in this case one has a unitary pertur- 
bation of a unitary theory for all fc > (except for the 
charge part, as usual). However, we also notice that the 
conformal weight of -0 is /i^, = s{s + I)/(fc -I- 2), which 
is fc/4 for s = k/2. Hence the perturbation is relevant 
for fc < 4, marginal for fc = 4, and irrelevant for fc > 4. 
As the perturbation is of order one, we cannot be cer- 
tain of the fate of each case. But suppose that the IR 
fixed point is the same as if the perturbation were weak, 
as discussed above. Then the cases fc < 4 lead to good 
topological phases. But for fc > 4, the return to the un- 
derlying CFT in the IR would imply that the models are 
not topological after all. We note that Blok and Wen 
proved ^35] that the adiabatic statistics is given by the 
monodromy of the conformal blocks, on the assumption 
that there was screening in the spin sector. This also en- 
tails that there be a well-defined spin of s within a "spin 
screening length" of a quasihole of the type labeled with 
spin s (when the quasihole is far from any others). The 
present scenario would imply that the screening assump- 
tion breaks down for fc > 4, and the statistics calculation 
then fails. In effect, the trial states would not possess spin 
s localized near each quasihole, instead this spin would 
be spread over the region containing several quasiholes. 
This would mean that the states are not topological in 
the spin sector, but behave as a gapless phase, or possi- 
bly a critical point. For the case fc = 4, further calcu- 
lations are needed to establish whether the perturbation 
is marginally relevant, marginally irrelevant, or exactly 
marginal (when weak). The latter two possibilities would 
both lead to behavior that is not that of a topological 
phase. 

The notion that screening in the spin sector breaks 
down for fc > 4 is intriguing, but by no means cer- 
tain. The Blok-Wen states are something like an SU(2)- 
invariant version of the mmn states above, with fc in place 
of TO — n, and corresponding conformal weights and filling 
factor (note that fc = I coincides with the n + l,n + l,n 
state). The relation can be made even closer by writ- 
ing the SU(2) level fc theory in terms of a scalar field if' 
and the parafermion CFT. In view of the mmn ex- 
ample which has screening in pseudospin except for very 
large to -I- n, it may well be that the Blok-Wen states are 
in the screening phase for spin, and possess non- Abelian 
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adiabatic statistics. 

There are other series of trial wavefunctions that have 
been proposed for particles with spin, and for which there 
is a special Hamiltonian for which the trial states are the 
zero-energy eigenstates [s^, [t^ . These behave similarly 
to either the RR or Blok-Wen series, and we will not 
discuss them further. 



4. Non-unitary theories and minimal models 

Now we turn to further examples, but start with gen- 
eral remarks about non-unitary underlying CFTs. In 
such cases, the perturbation by tp may still be relevant, 
irrelevant or marginal, when weak. We first consider the 
possibility of a flow to a massive phase. This would lead 
to adiabatic statistics given by the monodromy of the 
conformal blocks of the non- unitary CFT, as we have 
seen. However, there are difficulties with the interpreta- 
tion of this scenario. We have emphasized previously (Tlj 
that for a gapped bulk phase, the boundary is a massless 
field theory that is both unitary and conformal. If one 
has statistics in the bulk that is given by the monodromy 
of the blocks, then the same CFT would be expected at 
the edge [IJ . For a system with trial wavefunctions given 
by conformal blocks of a non-unitary theory, it is not 
clear what the candidate for this edge theory can be. 
However, there is certainly more than one CFT that pro- 
duces the same fusion rules, braiding in monodromy, etc, 
or in other words the same MTC. 

Here we wish to make a point that is more direct, and 
stronger. If the IR fixed point is massive, and so the 
MTC for the topological phase is that of a non-unitary 
UV CFT, then under some conditions this is not compat- 
ible with any unitary MTC (again, the MTC formalism 
applies to QH systems of bosons, and must be modi- 
fied for fermions, so the arguments here are for bosons, 
but the fermion case will be similar; from the wavefunc- 
tion point of view in QH systems, the two are very sim- 
ilar, differing only in the number of Laughlin-Jastrow 
factors). First, we should explain that there is a no- 
tion of a unitary MTC, which is precisely the correct 
concept to describe topological phases in 2 -f 1 dimen- 
sions that arise from quantum mechanical systems with 
a positive-definite inner product on the Hilbert space, 
and a Hermitian Hamiltonian (with respect to this inner 
product) [5J|. Unitarity of the MTC has certain conse- 
quences. For example, in any MTC each quasiparticle 
type has an associated "quantum dimension" . In a uni- 
tary MTC, all these quantum dimensions must be pos- 
itive. Next, for topological phases that originate from 
trial wavefunctions that are conformal blocks in some 
underlying CFT, the quantum dimensions can be cal- 
culated from the CFT. [Strictly speaking, we have only 
so far argued that when the IR fixed point is a mas- 
sive phase, we obtain the same fusion rules, braiding, 
and (probably) twist as in the underlying CFT. There 
are some remaining structures used in defining an MTC, 



which are needed to calculate the quantum dimensions; 
the derivation of these for trial functions given by confor- 
mal blocks is discussed in Appendix]^] In non-unitary 
CFTs, there are usually some negative quantum dimen- 
sions. [It can be shown using modular transformations 
[l3| that if any conformal weights in a rational CFT are 
strictly negative (h < 0), then some quantum dimen- 
sions are also negative. As a partial converse, if there 
are negative quantum dimensions, then there are fields 
other than the identity with h < 0. It does not seem 
to be known whether a non-unitary RCFT (MTC) must 
contain some negative conformal weights (resp., quantum 
dimensions).] Hence, when this occurs, no unitary MTC 
can describe the resulting behavior (see App.lXlfor further 
discussion). We conclude that such trial wavefunctions 
cannot flow to a massive phase. They must flow to some 
massless fleld theory, either back to the original CFT, 
or to some other non-trivial massless theory, or possibly 
exhibit exponentially-growing overlaps, which may indi- 
cate a more serious instability. We should note that for 
non-unitary theories a massless (scale-invariant) theory 
does not have to be conformal [T^] (for an example, see 
Ref. [Zl). 

As an interesting example, we consider the so-called 
Gaffnian trial wavefunctions that have been proposed re- 
cently ^7^. These are associated with one of the BPZ 
minimal CFTs, which in general are denoted M{p,p'), 
parametrized by a pair of coprime positive integers p, p' 
with p > p' [l7|. The Gaffnian functions are conformal 
blocks from the -M(5, 3) minimal model, which is non- 
unitary. The perturbing field ij) has conformal weight 
3/4, so the perturbation is relevant. The fusion rule for 
-0 is V' X = 1, so that these states may be considered to 
be paired states. Due to the relevance of the perturbation 
(when weak), we expect that the RG flow does not lead 
to the original CFT, and one might expect that the IR 
flxed point is massive (the trivial CFT) , which would give 
the same MTC as that of the underlying CFT. But the 
quantum dimensions da of the primary fields in a rational 
CFT are given by the first row (or column) of the unnor- 
malized version, denoted S*, of the modular S matrix; S 
was discussed above for the Majorana theory [which is 
M(4,3)], though in a different basis than here. One has 
S — DS, where D is the global dimension determined by 
D"^ = Sq^q' where a labels the distinct quasiparticle 
types, and we take the positive square root for D. For 
the minimal models, the primary flelds (f)r,s are labeled 
by ordered pairs of integers (r, s), with 1 < r < p' — 1, 
1 < s < p — 1 (the so-called minimal block in the Kac 
table) . There is an equivalence r —^p'^r, s —>■ p — s, and 
each inequivalent primary occurs just once in the theory. 
The conformal weights and central charge are given by 

M 



1 - 6 



pp' 



/^2 



{pr - p'sf -{p - p') 
App' 



(5.10) 
(5.11) 
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and the modular S matrix is 



Sr,! 



_(^_iy+sp+ra sin(7rrpp/p') sm{Tr sap' /p). 



(5.12) 

For the first column, p = a = 1, and the quantum dimen- 
sions are dr.s — DSr,s;i.i- Now for the case of M{5,3), 
we may put ?■ = 1 also, and for the Gaffnian trial states 
ip — (/>i,4- One finds that both di^^ and di^4 are negative, 
while the other two are positive = 1 for the identity 
always). This behavior is then impossible in a unitary 
MTC, and thus in a unitary topological phase. We con- 
clude that, if the RG flow driven by the perturbation 
does not lead back to the original M(5, 3) CFT, then it 
must go to some other non-trivial massless theory, and 
either way the system is in a gapless phase as a 2 + 1 
system (possibly, at a critical point, as discussed in Ref. 
[74| ) — or a worse instability takes place. We note that 
negative conformal weights, and hence negative quantum 
dimensions, occur in all the minimal models other than 
the unitary ones, for which p = p' -\- I. 

While on the subject of minimal models, we may con- 
sider other examples along similar lines as the Gaffnian. 
For p' > 2, the field at position (p' — 1, 1) [or equivalently, 
{l,p— 1)] in the Kac table always has Abelian statistics: 
the fusion rule for it is (/ip'_i_i x (ppi -i^i = 1, so that all 
these examples are paired states in the same general sense 
as the Gaffnian. [Forp' = 2, the field at (1, 1) or (l,p— 1) 
is the identity, so no interesting functions arise; the sim- 
plest non-trivial example is the Yang-Lee theory, M (5, 2) 
[ITj.] Using this field as -0, one obtains a set of trial wave- 
functions from conformal blocks for any minimal model. 
For many non-unitary cases, the perturbation is relevant, 
but a similar argument applies to all these non-unitary 
theories as for the Gaffnian (see App. However, for 
the unitary minimal models M{p' + l,p'), the scaling di- 
mension is hp'-i^i = (p'^ — 3p' + 2) /A, which is relevant 
(< 1) only in the case p' — 3, which is the MR example 
once again. The filling factors of these states are deter- 
mined hy ly"^ = (p'^ — 3p' + 2)/2 plus integers. The next 
simplest example in the unitary sequence is p' = 4, in 
which the minimal model describes the tricritical Ising 
model, which has many interesting features including su- 
perconformal symmetry. The perturbation is by a sub- 
leading "thermal" operator e" [l^l- Again, in general, 
although the perturbation is irrelevant when weak for 
p' > 3, it may be that these cases flow to a massive 
phase, representing a gapped topological phase in 2 -I- 1 
dimensions. In general, we are not aware of Hamiltonians 
for which these trial wavefunctions are zero-energy eigen- 
states, so it is not clear whether this approach leads to 
the construction of the corresponding topological phases. 
However, S. Simon has informed us that there appears 
to be a Hamiltonian that produces the "tricritical Ising" 
state 
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Recently, some large families of functions have been 
proposed as possible trial wavefunctions [t^, [zl]. It is 
not generally clear if these come from any CFT, and if 
they do not it is not clear whether they truly represent a 



topological phase. However, those in Ref. [77j include the 
Gaffnian and possibly some based on other non-unitary 
minimal models, and so it is likely that many of them 
can be ruled out (as far as gapped topological phases are 
concerned) by arguments similar to those given here and 
in App. \M 



5. Irrational theories 

In the preceding section, we considered non-unitary 
CFTs, but only in the context of rational CFTs, so that 
they still correspond to MTCs. We recall that a rational 
CFT contains a finite set of irreducible representations 
of a chiral algebra, and all the representations are fully 
reducible to direct sums of irreducibles fl§\ . In this sec- 
tion, we will consider some known examples in which the 
CFT may or may not be unitary, but is in any case not 
rational. 

The examples we consider here have been studied ear- 
lier. All of them are paired states. In each case there is a 
"special Hamiltonian" for which the trial states given by 
conformal blocks are the exact zero-energy eigenstates, 
and all were found to be at a critical point in the 2-1-1 
point of view. They are (i) the Haldane-Rezayi state, (ii) 
the permanent state, and (iii) the Haffnian state. First 
we will briefly discuss the CFT whose conformal blocks 
are the trial wavefunctions in each case. 

The Haldane-Rezayi (HR) state [t^ is a spin-singlet 
trial wavefunction for spin-1/2 particles. The CFT (other 
than the charge part) is the symplectic fermion theory 
with c = —2, which has the required SI2 symmetry [l|. 
This central charge corresponds to the values p = 2, p' = 
1 in the BPZ analysis. However, for the M(2, 1) minimal 
model, the minimal block of the Kac table is empty, as 
it is for all cases in which p' = 1. The fields used in the 
HR state [l^l are the symplectic fermion of conformal 
weight 1), and the spin field of conformal weight —1/8, 
which do lie in the Kac table, but are outside the (empty) 
minimal block. Consequently, in this non-unitary theory, 
the representations produced in operator products are 
not fully reducible. It was argued in Ref. 4] that the state 
corresponds to the weak- to strong-pairing transition for 
d — id spin-singlet paired states of fermions. 

The permanent state [l| is also a spin-singlet state. 
The CFT is the /3-7 system Q- On the torus, or in the 
presence of sufficiently many quasiholes, the dimension of 
the space of zero-energy states grows without bound as 
the system size increases (this is because unpaired bosons 
can occupy zero modes). The CFT is non- unitary, and 
also irrational because of the "picture-changing opera- 
tors". This system has been argued to correspond to 
the transition point between the spin-polarized Laugh- 
lin state, and a spin-density wave state, with the per- 
manent ground state itself cor resp onding in fact to an 
anti-skyrmion spin texture p6l . l8l| . 

Finally, the Haffnian state is a paired state for spinless 
(or spin polarized) particles in which spinless composite 
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bosons form pairs of angular momentum —2 [8l|, [S^l- In 
this case, the CFT is unitary, and consists once again of 
a chiral scalar </?'. The particles contain = d^' which 
is a U(l) current of conformal weight 1, and is a good 
conformal field, unlike ip' itself. Using Wick's theorem, 
the correlator of -0 in the case of the trial ground state is 
the Haffnian, 

( • • • dv\^N)) = Hf ( (^^T^) ' (5-13) 
where the Haffnian of a symmetric matrix M is defined 

by 

Hf My = S (Afi2Af34 • • • Mw_i,w) , (5.14) 

in which the symmetrizer S sums over all permutations 
that produce distinct pairings (i, j), similarly to the Pfaf- 
fian (from which the name is derived), but without the 
sign of the permutation. The quasihole functions contain 
the twist field r, around which the ■)/' field has a square 
root behavior of ope: 

V(z)t(0) -z-i/V(0) + .... (5.15) 

This field r features in orbifold theories based on Lp' ^ and 
has conformal weight 1/8 [s^. (There are many other 
possible twist fields with other fractional-power ope's 
|8a | , but only this one can be used to make single- valued 
functions of the particle coordinates when the charge sec- 
tor is included.) Like the permanent case, this state 
has highly-degenerate torus ground states and quasihole 
states, and is believed to be at a transition point, to one 
side of which is the Laughlin state [8l| . 

We will attempt to explain what is happening in these 
examples, focusing primarily on the Haffnian, as it uses 
a more familiar unitary CFT. To obtain unambiguously- 
defined systems, we will invoke here the special Hamil- 
tonians mentioned above. Such a Hamiltonian is an op- 
erator on the particles. The detailed form for each case 
will not be important here. The point is that in seeking 
zero-energy eigenstates of such a Hamiltonian, we are 
forced to fix the ope's of the fields V'e that correspond 
to the particles, in a form that of course depends on the 
Hamiltonian. 

For the Haffnian state, the special Hamiltonian implies 
that the CFT should have the ope's of the U(l) current 
'0 — dip': 

il){z)tl){Q) ^l/z^ + .... (5.16) 

If the scalar field p' were compactified with radius R, 
then the scalar theory would contain such fields as 
exp lap' for a discrete set of real values of a (integer 
multiples of one value, related to R [T^). For the orb- 
ifold of this rational theory, there are similar fields, and 
the ope's among if) and r are independent of the radius R. 
Because the Hamiltonian can only determine the ope for 
'0, in counting states we must consider all possible val- 
ues of R and a. As R goes to infinity, the exponentials 



of p' can be expanded in powers of p' , and insertions 
of these appear to produce the conformal blocks that we 
described as bosons occupying zero modes (on the torus 
or in the presence of quasiholes). In the limit, the CFT is 
clearly not rational. However, one also expects that only 
a rational theory can correspond to a topological phase; 
in particular, the number of quasiparticle types should 
be finite, otherwise one suspects that the system will be 
gapless. 

Similar arguments apply to the symplectic fermion and 
/3-7 systems. For symplectic fermions, being fermions, 
the use of different radii does not directly apply. But 
there is still an infinite set of "rational" versions, in the 
sense of possessing an extended chiral algebra that has 
only a finite number of representations (though they still 
fail to be semi-simple) |84| . 

It is interesting also to pursue some of the other argu- 
ments above in the present cases. We notice that for the 
HR and Haffnian cases, the perturbing field if} is weight 
one, and so is marginal when weak (as a perturbation of 
the CFT without the charge part). The same scaling is 
found in the pairing function g of the p-l-ip paired states 
at the critical point [3] (though then g is not meromor- 
phic, and does not correspond to a conformal field the- 
ory). If we consider this perturbation as weak, and write 
p'{z,'z) = p' [z) ^p'{z) in the non-chiral theory, then in 
the action we have Xdp'dp'. This is the same as the un- 
perturbed action, and can be absorbed by rescaling p' . 
Thus there is no effect on the theory of the perturbation 
(though there would be if the field were compactified). 
The perturbation is not only marginal, it is redundant. 
The same argument goes through for the HR theory also. 
We should note however that since the perturbation is 
actually of order one, there may be concerns about the 
correctness of the argument. At any rate, other irrele- 
vant operators will be generated, whereas the rescaling 
argument would suggest that there are none. We cannot 
enter into this further here. If we accept the argument, 
then the IR fixed point is the same CFT with which we 
started the RG flow, and we argued above that this must 
correspond to a gapless phase. This is consistent with the 
quite different arguments that these states correspond to 
critical points, given previously. 

As a related side remark, we return to examples such 
as the Yang-Lee CFT mentioned previously. The ground 
state wavefunctions (which contain — 1 ioi the CFT) 
can be produced by the usual pseudopotential Hamilto- 
nian that produces the Laughlin states, which of course 
are the same functions. One is free to write down the 
functions derived from the Yang-Lee theory, including 
the only non-trivial field r = 0i,2- But as functions of 
the particle coordinates, the trial functions are indepen- 
dent of the positions w of these insertions, which conse- 
quently can carry zero charge in the charge part. Thus 
these trial functions are not linearly independent of the 
usual Laughlin ground and quasihole states, and are in 
fact the Laughlin functions times functions of the r loca- 
tions alone (thus the theory is a direct product). This is 
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why the Yang-Lee theory leads to nothing new. We no- 
tice that the field ipe (with its dual) generates the chiral 
algebra of the Laughlin state [H , which never includes the 
stress tensor of the Yang-Lee theory. In general, we ex- 
pect that the trial functions based on conformal blocks of 
a CFT are linearly independent as functions of the parti- 
cle coordinates only if the particle field generates the full 
chiral algebra of the CFT (including the charge sector) . 

Put another way, a special Hamiltonian determines the 
field Ve which is part of some CFT. The CFT (includ- 
ing the charge part) that one should use to analyze the 
problem is the "minimal" one, that is the one with the 
smallest possible chiral algebra, and that means that ipe 
and its dual should generate the chiral algebra. For the 
Haffnian and other examples above, the resulting chiral 
algebra is however too small to define a rational theory 
(but extensions exist that correspond to a rational the- 
ory, at least in the sense of a finite number of types). 



VI. CONCLUSION 

The main message to emerge from this work is that 
when conformal blocks are used as trial wavefunctions, 
whether in QH or other systems, the properties of the 
topological phase (if any) that they represent are directly 
determined by, and actually equal to, those of the under- 
lying CFT. (This has frequently been assumed, rather 
than demonstrated, in the literature.) That is, the holon- 
omy (such as adiabatic statistics) is equal to the mon- 
odromy (analytic continuation of blocks) if and only if 
the 2D perturbed CFT is in the massive phase. If not, 
then there are signs that the system is gapless as a 2 -|- 1 
phase of matter (though we did not address Hamiltonians 
for our states). 

When the trial wavefunctions represent a topological 
phase, our results almost completely constrain it. The 
loose end is that the twist (the effect of rotating a quasi- 
particle about its center) should also be calculated by 
adiabatic transport, and this has not been done. How- 
ever, it should be possible by adiabatically varying the 
metric of the 2D surface. A first step in this direction 
is the determination of the Hall viscosity, which deter- 
mines an Abelian contribution to this transport that is 
present in all cases. The fact that this Hall viscosity is 
itself related to the density of conformal weight in the 
trial ground states strongly suggests that the conformal 
weight of the quasihole field will also emerge in the twist, 
as expected. 

If the latter calculation is completed, it will remove the 
final possible loophole in our argument that non-unitary 
RCFTs that contain negative quantum dimensions (or 
conformal weights) do not give rise to a topological phase 
(or unitary modular tensor category) when their blocks 
are used as trial wavefunctions; instead they should ap- 
parently be gapless, or suffer a worse instability. 

Another loose end of this work is that the BCS trial 
wavefunctions for states containing vortices have not 



been shown to reduce to the conformal block form at 
large length scales. Instead, the blocks were used di- 
rectly (see Sec. ITV)) . It would be desirable to show this, 
so as to broaden the domain of explicit calculations to 
the wider class of BCS functions. 

Finally, we note that the rather simple form of all these 
results seems to be related to the fact that the underly- 
ing particles are Abelian, as is the perturbing field in the 
CFT point of view. It would be of interest to consider 
states formed of a non-zero density of non- Abelian parti- 
cles (which could be useful in hierarchical constructions 
of QH or anyon superfluid states) , but this is clearly much 
more difficult to analyze. 

Note added: Recent work has tied up the loose end 
mentioned above, by deriving the twist for quasiholes in 
the trial states from their spin, which can be obtained 
by adiabatic transport in the presence of curved space, 
such as on a sphere [s^. This completes the argument 
that conformal blocks from any non-unitary CFT con- 
taining negative quantum dimensions cannot correspond 
to a topological phase. 
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APPENDIX A: REMAINING DATA OF AN MTC, 
AND NON-UNITARY CASES 

In this Appendix we discuss the remaining structures 
needed for the definition of an MTC, other than the fu- 
sion rules, braiding, and twist, which have been discussed 
in the main text. As we do not intend to give a full ex- 
position of MTCs, we cannot be entirely self-contained, 
and must make use of the references for full technical de- 
tails. Accordingly, this Appendix is intended for readers 
more expert about such matters. However, its contents 
are mostly general definitions, which many readers might 
be willing to take on trust. We do address the Gaffnian 
example again at the end. 

The main goal for the definitions is to be able to calcu- 
late the quantum dimension associated with each quasi- 
particle type, from trial wavefunctions given by confor- 
mal blocks of some CFT. This requires that the full struc- 
ture of a ribbon category [53, 54, 55] be available. (A 
MTC is a ribbon category in which the S matrix, which 
can be defined in a ribbon category, is invertible.) 

Informally, the quantum dimension for a type of quasi- 
particle is defined as the amplitude for a process in which 
the quasiparticle, of type a say, and its antiparticle are 
created from the vacuum, separated, and then annihi- 
lated again. The motion of quasiparticles can be assumed 
to be adiabatic. The definition is not vacuous because 
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the creation and annihilation processes used have to be 
normahzed to satisfy other topological criteria. We note 
immediately that in a unitary (i.e. quantum-mechanical) 
theory (in 2 + 1 dimensions), the creation and annihila- 
tion processes are adjoints to each other, and the quan- 
tum dimension is then the modulus-square of the creation 
amplitude. Hence it must be positive in the ribbon cate- 
gory describing a quantum-mechanical topological phase 
(or unitary TQFT). The latter corresponds to a unitary 
ribbon (and modular tensor) category in the sense of Tu- 
raev [54| . 

In the language of tensor categories, the distinct quasi- 
particle types a correspond to isomorphism classes of 
simple objects Va- The trivial or identity quasiparticle 
type corresponds to Vb = 1. More general objects 
can be formed by direct sums and by tensor products, 
written F ® and V (E)W respectively (the notation is 
based on that for vector spaces, but for us these objects 
do not have to be thought of as such). The fusion rules 
correspond to the decomposition of a tensor product of 
simples as (isomorphic to) a direct sum: 

Va®V0 = ^N2f,V„ (Al) 

7 

where we write NV for V(BVQ) ■ ■ ■ with N terms V in the 
sum. There are isomorphisms li^V = V,V^l = V,ioi 
all V. Although the objects may not be vector spaces, 
the morphisms (or maps) between objects, say from V to 
W, form a complex vector space Hom(F, W). 

The tensor product is in general not strictly associa- 
tive, but there are "natural" isomorphisms F between 
different orders of brackets: 

Fuy.w ■iU(E)V)(g)W ^U(g>{V(g>W). (A2) 

These may be built from the F "matrices" for products 
of simple objects (replacing U, V, W, by Va, Va, V^), 
FafSj. The tensor product on either side decomposes as 
the same sum of simple terms, and after choosing bases, 
Faf3j becomes a matrix in this space |55| . which has di- 
mension J2e,s ^ap^n = ^^-y^L'- t^crc are no 
morphisms between distinct V^s, the F matrices break 
into block diagonal matrices F^^^ of size N^^N^^ = 

T.e'^p-y^L'- Fa/B-y = J2s ^^p-y ^s a tensor product 
is isomorphic to a direct sum, say Va 'S)V/3 = (BtN^^V^, 

the F^ij^ have blocks that can be labeled F^^^^ . The 
decomposition to blocks Ff^^^ is well-defined up to iso- 
morphisms, and these maps become actual matrices (or 
simply numbers if N^^Nf^ and N^^N^^, equal 1) only 
after fully specifying a basis. The F matrices must sat- 
isfy some consistency conditions (pentagon equations) 

[m, M,^]. 

Braiding is strictly defined only for two objects not 
separated by any parentheses: cv,w '■ V(^W W(^V for 
any V, W. The braiding isomorphisms c may be built up 
from simple components, and (choosing bases) give rise 
to matrices c^g on the sum of terms isomorphic to V^ in 



Va'^Vp. In the case of a RCFT, the eigenvalues of these 
can be found from the ope of (j)a with The braiding 
and F matrices must satisfy the hexagon equations. For 
a ribbon category arising from a CFT, the pentagon and 
hexagon equations are satisfied automatically. 

The formula for the number of simple terms, which 
when the ribbon category comes from a RCFT is the di- 
mension of the space of conformal blocks for correlators 
in the plane of corresponding fields, illustrates how in an 
iterated tensor product of simples, of any length, each 
lexicographically valid way of inserting brackets corre- 
sponds to a labeling of states (the labeling is only a par- 
tial labeling if some iVj^ > !)• For each such choice there 
is a basis of the space of conformal blocks. The F ma- 
trices then describe the change between two such bases. 
From the physical point of view, these basis changes are 
"passive" transformations, whereas the braid and twist 
operations are "active". For a CFT, the F matrices can 
be read off from the behavior of the conformal blocks. To 
consider this, it is useful to place all fields (or quasiholes) 
(f>ai on the X axis at positions wi = Xi. If three adjacent 
such fields are types a, /3, 7 as above, with a at a; = 0, 7 
at X = 1, and /S at x, < x < 1, then as a; ^ one would 
use the ope and analyze the block as on the left-hand side 
of the definition of F. As a; 1, one analyzes as on the 
right-hand side. The conformal blocks then determine 
the F matrices. The (w-independent) linear transforma- 
tion between bases remain valid even when the wis are 
not asymptotically close together, which is useful as we 
want them further apart than the screening length. In 
addition, they remain unchanged in the presence of added 
fields ipe, which are in the chiral alg ebra of the full CFT 
(including the charge sector) P, [Zg. For us, the confor- 
mal blocks are the trial states, and we use one of these 
bases, as convenient, and transform between them using 
the F matrices which are determined by our choice of the 
trial functions themselves (elsewhere in this paper, we did 
not need to specify our bases explicitly). For the braid 
group representation, all matrices must be expressed in 
a single basis, and F moves are used to transform this to 
one appropriate for an application of any desired braid 
generator c or c^^. 

The other structure needed in a ribbon category is the 
set of duality maps. These correspond to the creation and 
annihilation processes mentioned above. For each object 

V there is (an isomorphism class of) dual object(s) V* , 
and morphisms ey ■ V* (E) V ^ 1 (evaluation), iy : 1 ^ 

V ^ V* (co-evaluation). The order of factors in these 
tensor products should be carefully noted. These maps 
are required to satisfy relations (as maps V —i- V): 

idy (E) ev o Fvyy o iy (g) idy = idy, (A3) 

where idy is the identity morphism on V, and o is com- 
position of maps, which act from the left. A similar rela- 
tion with V and V* interchanged, and the order of tensor 
factors reversed, is also a requirement. When viewed as 
creation and annihilation processes in spacetime, these 
relations allow us to straighten out worldlines that some- 
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times reverse course and run backwards (by a creation 
or annihilation event). Thus for worldhnes with a zigzag 
of this form, the zig and the zag cancel (however, the 
similar identity for the mirror image zigzag is not a con- 
sequence of these identities). For simple objects, we will 
sometimes write V* as Vq. . (Here and below, for simplic- 
ity we do not explicitl y w rite the maps ly '■ 1 ® V V 
and ry : V ^ 1 ^ V [53[, which strictly should appear 
together with their inverses on the left hand side of this 
relation.) 

Finally, we recall the twist 9 which is a map By : V 
V for each object V. For simple objects Va, we can write 
9va, as a number, which is e'^^^'^" when a ribbon category 
is obtained from a RCFT. The twist must obey 

dy^w — cw,vcy^w(9y ® 9w), (A4) 

and also 9^ — idi — 1, By — {9y)* (the latter uses the 
dual map, which is the transpose matrix — not complex 
conjugation!). 

The quantum dimension for any object V can now be 
formally defined as a composite map from 1 to 1, and 
hence can be represented as a complex number: 

qdimy eyCyy,{9y ® idy. )iy, (A5) 

which obeys qdimy^^y = qdimy + qdim^y and 
qdimy^^ = (qdimy)(qdim^y). For simple objects Va we 
define da = qdimy^ (and note that do = !)• The quan- 
tum dimension can be defined more simply, using maps 
e'y , defined by e'y^ = ey^Ca,a' {9y^ (g)idyj ) -.Va^V* ^ 1, 
and similarly i'y. Then one has qdimy = e'yiy = eyi'y. 
As mentioned above, in a unitary ribbon category, e'y 
(ey) is the adjoint of iy (resp., iy), and so c?q > for all 
a. 

Now we address these structures for our conformal 
blocks viewed as trial wavef unctions, in the case when 
the flow in the 2D theory is to a massive IR fixed point. 
Then we have argued in the main text that the overlap 
matrix is proportional to a sum of identity matrices, in a 
basis in which the braid group acts by unitary transfor- 
mations (for this, any one of the bases above is suitable), 
with proportionality constants that are independent of 
quasihole positions w, provided they are well separated. 
Again, we will consider here wi on the real axis. The 
F matrices allow us to transform the basis in this state- 
ment. We can assume that a basis exists in which the F 
matrices are unitary, as otherwise there will be inconsis- 
tencies. For the ey maps, we will adiabatically transport 
quasiholes corresponding to Va and V* together, but not 
closer than the screening length. By definition, from far 
away they can be viewed as a single quasihole of types oc- 
curring in the fusion rule for Va with V*^ and this must 
include the identity 1, just once. (We must digress to 
point out that the fusion of the negatively-charged quasi- 
holes in the QH effect will produce a negatively-charged 
result. But the destruction of the underlying particle 
produces a quasihole of exactly this type, and is viewed 
as equivalent to the identity as it has trivial braiding 



with all quasiparticle types.) There is also a trial state 
with the pair replaced by the quasiparticle type 1, which 
can also be normalized. The map ey^ is a map between 
these, and can be defined without further consideration 
of any short-scale physics of the states. It is tempting 
to say that the map ey^ is the map between these two 
normalized states, times some phase factor. But this is 
not generally correct, and the reason is that the i and 
e maps must be normalized by relation (|A3[) above. As 
we will see, the F matrix element is something like 1/da, 
while e and i are like ^/d^. 

In terms of conformal blocks viewed as trial wavefunc- 
tions, if a quasiparticle type a* that is suitable to be the 
dual of type a exists for all a, then it can be shown that 
Hom(l, Va (Xi V*) and Rom{V* (g) Va, 1), to which iy^ and 
ey^ belong, are one-dimensional. Then the normaliza- 
tion condition allows the definition of iy^ and ey^ maps 
to be varied only by multiplication by scalars that are 
inverses of each other for each a, which cannot change 
the quantum dimensions. This is true for the CFT itself, 
also. But the quantum dimensions can be evaluated in 
the CFT, and if any is negative, then there is no way to 
avoid a confiict with the quantum mechanical definition 
which leads to positive quantum dimensions. 

Here all structures other than the duality maps were 
already fixed, either from properties of the conformal 
blocks used, or by adiabatic transport that agrees with 
the monodromy of the blocks. However, there was some 
uncertainty about 9 (which is independent of the choice of 
-Fs, c, i, e), because the adiabatic calculation of the twist 
by varying the metric has not been performed. Perhaps 
the system can be rendered consistent by the presence of 
different values 9'y , with the other structures which we 
have either defined or calculated (c, F, i, e) unchanged. 
A different choice for the values of 9 on simple objects can 
sometimes change the values of the quantum dimensions, 
making some positive, as we will see in a moment. But 
this cannot change the value of the F matrices, or the 
normalization condition for i and e. We will see that, in 
many cases, no consistent choice exists for the twist that 
makes all quantum dimensions positive (for the given Fs, 
c, i, e). 

We already know one set of values of 9y^ that is con- 
sistent, which come from the CFT. A different choice 9' 
must still satisfy relation (|A4p . so that 

(^Vc,(»VfJ (^Vc(SVfi = {(^v^/SyJ{9y^/9y^). (A6) 

Moreover, 9'^/9x = 1. This implies that the map 
(j)a — > {Ov^lQva)4'a is an automorphism of the ring de- 
fined by the fusion rules (the Grothendieck ring). This 
and 9'y, = 9'y^ then imply further that {9'y^/9yJ'^ = 1, 
so 9'y /9y^ ~ ±1 for all a. Now as the definition of 
quantum dimension used 9, if we keep c, z, e fixed we 
may be able to reverse the signs of some of the quantum 
dimensions. 

As an example, for the Gaffnian state, where the CFT 
apart from the charge sector is the M(5, 3) minimal 
model, the fusion rules are the same as for SU(2) level 3, 
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and [using the notation of SU(2) spin s, a = s = 0, 1/2, 
. . . , 3/2 for the fields], the only automorphism is given 
by multiplication by (—1)^*. This leaves di (di^s in the 
notation in the main text) still negative. More generally, 
for the class of states mentioned in Sec. IVC4[ in which 
the CFT is a BPZ minimal model, the twist of the field 
(/>r,s can be changed by multiplication by (—1)' "^, or by 
(— or both. Using the formula (|5.12p for the mod- 
ular S matrix of the minimal models, for the non-unitary 
cases p ^ p' + 1 there are some negative quantum dimen- 
sions, and one can show that not all of the negative ones 
become positive under any of these changes in the twist 
(indeed, additional ones are generated). Thus the fam- 
ily of trial states constructed from the minimal models, 
mentioned in Sec. IV C 4[ fails to produce valid (i.e. uni- 
tary) topological phases except in the case of the unitary 
minimal models M(jj' + l,p'). 

We also record here a basis-independent relation in- 
volving the part of the F matrix that enters the normal- 



ization condition, together with braiding and twist ma- 
trix elements, and the quantum dimensions. The result, 
which is close to Moore and Seiberg [3| (and references 
therein), eqs. (7.12) and (C.15), is 

(Co,q) i'^a,a') ^aa'afiO — (^a/da- (A7) 

All the factors in this expression are simply numbers. 
(We note that c°»„c°q,, = 9~'^ and cJ^^qCo.q = 1 in any 
basis, so up to some choices of basis the result says that 
Paa'afio ~ l/^a-) H wc Can change the sign of 9a as 
discussed above, then da also changes sign, while the 
left hand side is invariant. The relation is obtained by 
evaluating 

Ca.o o (idy„ ® e'yj O (idy^ (8) „-) ° ^Sa-a ° «) idy^), 

(A8) 

which is a map from 1 (g) to itself, in two ways, using 
the relation (|A3|) . 
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